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ABSTRACT 


The Interplanetary Program to Optimize Simulated Trajectories (IPOST) is 
intended to support many analysis phases, from early interplanetary 
feasibility studies through spacecraft development and operations. The 
IPOST output provides information for sizing and understanding mission 
impacts related to propulsion, guidance, communications, sensor/actuators, 
payload, and other dynamic and geometric environments. 

IPOST models three degree of freedom trajectory events, such as 
launch/ascent, orbital coast, propulsive maneuvering (impulsive and finite 
bum), gravity assist, and atmospheric entry. Trajectory propagation is 
performed using a choice of Cowell, Encke, Multiconic, Onestep, or Conic 
methods. The user identifies a desired sequence of trajectory events, and 
selects which parameters are independent (controls) and dependent (targets), 
as well as other constraints and the cost function. 

Targeting and optimization is performed using the Stanford NPSOL (Non- 
linear Programming Stanford Optimization Laboratory) algorithm. IPOST . 
structure allows sub-problems within a master optimization problem to aid in 
the general constrained parameter optimization solution. An alternate 
optimization method uses implicit simulation and collocation techniques. 

IPOST has been developed by Martin Marietta under contract NASI- 18147 to 
NASA/Langley. IPOST runs on a SUN and a Silicon Graphics computer. 
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FOREWORD 


This report describing the formulation of the Interplanetary Program to 
Optimize Sim ulated Trajectories (IPOST) is provided in accordance with. 
Statement of Work Part 4.9 of NASA Contract NAS1-18147. The report is 
presented as follows: 

Volume I - Interplanetary Program to Optimize Simulated 
Trajectories User's Guide 

Volume II - Interplanetary Program to Optimize Simulated 
Trajectories Analytic Manual 

Volume III - Interplanetary Program to Optimize Simulated 
Trajectories Programmer’s Manual 

Volume IV - Interplanetary Program to Optimize Simulated 
Trajectories - Sample Cases 

This work was conducted under the direction of Mr. Richard W. Powell of the 
Space Systems Division, National Aeronautics and Space Administration, 
Langley Research Center. 

A number of people contributed to the development of IPOST and to this 
report. The first issue (Reference 1) was delivered in March, 1990, with the 
support of Garry Brauer, Sandy Fitzgerald, Phil Hong, Perry Kent, Mac 
Milleur, Dave Olson, Fred Petersen, and Beth Swickard. The current 
revision and expansion of IPOST and this report were developed by Garry 
Brauer, Phil Hong, Perry Kent, Dave Olson, Larry Rockwell, Brian Sutter 
«nd Candy Vallado with acknowledgement to Dick Powell, Scott Striepe and 
Prasun Desai at Langley Research Center for their inputs, suggestions and 
support. 
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1.0 INTRODUCTION 

The Interplanetary Progr am to Optimize Simulated Trajectories (IPOST) is 
intended to support many analysis phases, from early interplanetary feasibility 
studies through spacecraft development and operations. The IPOST output 
provides information for sizing and underst an d ing mission impacts related to 
propulsion, guidance, co mmuni cations, sensor/actuators, payload, and other 
dynamic and geometric environments. 

Much of the overall architecture for IPOST has been derived from the Program to 
Optimize Simulated Trajectories (POST) (Reference 2). Indeed certain POST 
parameters and capabilities have been incorporated into IPOST to aid in POST- 
IPOST user compatibility. IPOST has extended trajectory capabilities to target 
planets and other celestial bodies with intermediate and velocity correction 
maneuvers. IPOST capabilities and limitations are summarized in Table 1-1. 


FEATURE 

CAPABILITY 

Optimization method 

Explicit (Master/subproblems), Implicit (collocation) 

Optimization algorithm 

NPSOL 

Optimization 

parameter* 

AV magnitude, mass, time, — 

Maximum controls 

25 (Master), 45 (subproblems), 1700 (collocation) 

Control parameters* 

Values of event criteria, AV, arrival conditions, thrust, . . 

Maximum targets 

25 (Master), 45 (subproblems), 1700 (collocation) 


Time, position, velocity, orbital conditions, . . . 

Targeting method 

NPSOL, Newton-Raphson, special Onestep 

Sensitivity matrix 

Finite differencing, analytic for special interplanetary 
targeting 

Maximum events 

100 

Event criteria* 

Time, distance, speed, closest approach, . . . 

Event activities 

Info, impulsive AV, launch, orbit insertion, mass jettison 

Maximum maneuvers/ [ 
subproblems 

15 


Conic, Onestep, Multiconic, Encke, Cowell, implicit 

Planetary bodies 

Sun, nine planets, Earth's moon, any user-defined bodies 

Ephemeris 

Analytic, precision (JPL) 

Trajectory perturbations 

Central body, perturbing bodies, radiation pressure, J2, 
aerodynamics, thrust 

Input/Output frames 

Ecliptic or planet equator. Mean 1950 or Mean 2000 

| * User selectable 1 


Table 1-1. IPOST Features/Capabilities 

IPOST, along with members of its family, such as POST and 1PREP, can analyze 
and support almost every activity associated with space exploration. 
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IPOST is event driven. That is, the user defines a sequence of events which are 
executed in the simulation process. The events can be triggered by different 
criteria, such as absolute or relative time, distance from a body, or propellant 
consumption. At the event times, various activities can be initiated or ter mina ted, 
such as employing a different thrust steering law, changing trajectory propagators 
or propagation step size, performing an impulsive delta velocity maneuver or 
jettisoning a probe or stage. 

The time period between two contiguous events is called a phase. Trajectory 
propagation takes place in each phase. Five types of propagators are a vailabl e 
(listed in order of increasing accuracy and decreasing computational speed): Conic, 
Onestep, Multiconic, Encke, Cowell. Propagator selection depends upon user needs, 
such as simple fast simulations for parametric feasibility analysis, or precision 
detailed trajectories to support subsystem design. 

IPOST can run a single trajectory simulation or it can run multiple simulations. 

For multiple simulations, one can run a parametric scan and/or an optimization 
mode. The search mode will vary one parameter, such as planetary arrival time, 
over a specified interval and increment size, and perform a simulation (or 
optimization) for each search parameter value. 

The optimization mode will optimize a user cost/objective function, such as 
m a ximum mass that can be placed in a desired orbit, subject to user-specified 
constraints. The constraint variables, such as periapsis altitude or orbital 
inclination, are called dependent variables or target parameters. The parameters 
which are free to vary, such as maneuver delta velocity (AV), are called independent 
variables or control parameters. As part of, or instead of, optimization, trajectory 
targeting can be performed. In this case, there is no cost function and the IPOST 
problem reduces to finding a set of control parameter values that meet specified 
target parameter conditions. 

Generalized targeting and optimization uses the Stanford NPSOL algorithm. For 
certain types of problems, a trajectory decomposition method is available. There is 
a master optimization process which requires that the trajectory be divided into legs 
or sub-problems. Each subproblem is an optimization problem in itself, wmtniwing 
controls, constraints and an (optional) objective function. A special application of 
decomposition is the Interplanetary Targeting and Optimization Option (ITOO). 
This technique uses analytical partials generated during nominal trajectory 
propagation to determine minimum AV (or mass) trajectories, usually for gravity 
assist (swingby) missions. 

In addition to the classic method of explicit optimization, there exists an option to 
perform implicit optimization us in g the collocation method. In this case, each phase 
is divided into independent segments which are allowed to vary subject to 
intersegment continuity and the equations of motion. Optimization using 
collocation is less sensitive to faulty initial guesses, but requires much greater CP 
time t h a n explicit optimization to achieve the same level of accuracy. 
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IPOST input is via three namelists: $TOP, $TRAJ and $TAB. $TOP contains a 
description of the targeting and optimization problem. It must be input first. 
$TRAJ mntjaing data that describes each mission event/phase. It mu6t follow $TOP 
and there must be one $TRAJ for each event. $TAB is used to input tabular data 
anpV> as thrust vs. time or drag coefficient vs. mach number and angle of attack. 
Input and output units are metric. 
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[SSION ANALYSIS TOOLS 


The general area of mission analysis indudes both atmospheric and exo- 
atmospheric flight. Although many software tools have been developed, only a few 
have proved to be useful for practical engineering studies. The POST famil y has 
evolved over the last 20 years to provide government and industry users with the 
capabilities to analyze all regimes of flight. 

2.1 THE POST FAMTT.V 


POST was developed in 1970 as a Space Shuttle Trajectory Optimization Program. 
Since that time, the program has been significantly improved with additional 
capabilities in die areas of vehicle modelling, trajectory simulation, and targeting 
and optimization. As the need has grown to analyze diverse missions, POST has 
spawned a number of related programs. In 1988, IPOST development was initiated 
in order to analyze future interplanetary missions. The complete POST family of 
tools consists of IPREP (Interplanetary PREProcessor), LPREP (Lunar 
PRE Processor), IPOST, TESTE (Testbed), GMAP (General Mission Analysis 
Program), NPOST and 6DPOST. (see Figure 2-1) 

POST is capable of simulating and optimizing trajectories for aerod ynami c vehicles 
operating in the vicinity of a single planetary body. POST can optimize a variety of 
trajectory problems, such as launch and ascent to orbit or escape, entry and landing 
of a probe, optimal aircraft cruise profiles, and vehicle guidance and control. 
Previous applications of POST have included launch vehicle and trajectory design 
for Voyager (Earth liftoff to hyperbolic escape). Mars ascent vehicle design (MRSR - 
Mars Rover and Sample Return), Manned Earth entry (return from Mars) 
simulation, and single stage to orbit (e.g. NASP - National AeroSpace Plane). 

NPOST expanded on POST by using the NPSOL algorithm rather than PGA 
(Projected Gradient Algorithm) to provide an alternative technique which proved to 
be more robust. NPOST is currently used to optimize Titan II and Titan IV ascent 
trajectories. 


6DPOST added rotational degrees of freedom and guidance feedback to POST, along 
with an improved input processor. An example of a 6DPOST application is Titan IV 
guidance coefficient and autopilot design for flight profiles requiring wind relief. 

GMAP, developed in the mid- 1970's, performs orbital mission optimization. 
Emphasis has been directed toward the evaluation of Earth synchronous missions, 
sun synchronous missions, drag influenced low altitude missions, and Shuttle 
applications. GMAP is applicable to a wide variety of mission analysis problems 
such as maneuvering and stationkeeping. Maneuver models include impulsive and 
fimteburn with various guidance and targeting schemes. GMAP contains a passive 
orbiting target vehide to allow for simplified rendezvous analysis. There is no 
interplanetaiy capability beyond that of Ejecting a given payload to the desired 
outgoing asymptote. Examples of GMAP applications indude multi-vehide 
rendezvous simulation and multiple impulse LEO to GEO orbit transfer (minimum 
fuel). 
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TESTB is a software testbed containing simple 2 degree of freedom simulation 
models of launch vehicles and interplanetary missions. It allows for rapid 
prototyping and characterization of trajectory optimization concepts. Before 
implementation of new algorithms into highly complex programs such as IPOST 
and POST, TESTB can be used for proof of concept demonstration. TESTB has been 
used to compare NPSOL with PGA , as well as with stochastic optimization 
methods, and for collocation development with trajectory applications. 

IPOST was developed from the POST architecture in direct response to 
interplanetary space exploration needs. Many of the trajectory simulation models 
were changed, while retaining the basic constrained parameter optimization 
capabilities of POST. As IPOST development progressed, it was dear that the non- 
linearities of interplanetary flight required feasible initial guesses for the 
optimization control parameters. Thus were bom IPREP and LPREP. 

IPKEP and LPREP provide rapid grid-searches across launch and arrival windows. 
These scans of mission opportunities result in selection of a minimum energy 
solution. For speed of execution both programs use simple propagators (no 
integration) and minimal perturbing forces. IPREP/LPREP can be used as 
standalone mission opportunity characterization tools and/or to provide initial 
conditions for IPOST. 

IPREP/LPREP model planetary phases of launch, arrival, flyby and orbit. 

Planetary constraints indude event dates, flight time and trajectory/orbit 
conditions. Mission options indude ballistic partial and multiple revolutions (Type 
I, II, HI, etc.), low thrust (e.g. SEP and NEP), high thrust (impulsive), and gravity 
a88 *8t. Applications of IPREP/LPREP have induded Voyager mission verification, 
solar probe mission design, and space station departure with lunar free return 
mission design. 

The POST family has been developed to operate as a set of coherent mission 
analysis tools. For example, recent interest in manned space exploration has 
identified la ndin g on Mars as a top national priority. One representative mission is 
illustrated in Figure 2-2. 
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Figure 2-2. Mission to Mars 

Application of IPREP, IPO ST and POST to the representative Mars mission is 
shown in Figure 2 - 3. A number of alternative mission options can be examined 
including nuclear thermal propulsion, Venus gravity assist for outbound or return 
legs, abort trajectories and modes, and propulsive vs. aerobraking orbit capture. 
IPREP is initially used to define Earth launch windows and escape energy, Mars 
arrival time and orbit insertion energy. Mars stay time and escape energy from 
Mars, and Earth return windows and orbit capture energy. This information is 
used in IPOST to define precise trajectory and performance data for each major 
mission leg: Earth to Mars, Mars oibit, Mars to Earth. POST is used to simulate 
and design the atmospheric flight portion: Earth liftoff to escape. Mars entry and 
landing, Mars surface ascent to orbit, and Earth entxy to touchdown. Together, all 
three programs provide a synergistic simulation and optimization capability which 
complements the overall system design and analysis. Input and output for each 
program has been designed to allow smooth connectivity across programs. 
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Manned Mars Mission 

• Opposition Gass 

• Launch in 2005 

• Return in 2008 


IPREP 


i 


IPOST 


i 


POST 


• Examine dates and energy requirements for Earth 
departure and Mars arrival, Mars departure and 
Earth return 

• Evaluate nuclear propulsion, Venus swingby and 
deep space maneuver options 

• Select date, C3, etc. for Earth to Mars and 
Mars to Earth with Venus swingby return 


• Optimize: Earth-Mars- Mars-Earth mission 
(maximize return to Earth mass) 

• Controls: Earth and Mars C3 departures 
Earth VINF return. Mars VTNF approach. 
Nuclear propulsion: thrust magnitudes, 
directions, segments (on/off times), power 

• Constraints: Mars stay time. Mars orbit and 
mass to Mars surface 

• Determine geometry, dynamic conditions and 
habitation conditions 


• Optimize: Earth LV GLOW (gross liftoff 
weight). Mars ascent vehicle mass 

• Constrain max q, g's , escape energy 


Figure 2-3. Interplanetary Mission Analysis Example Problem 




2.2 COORDINATE SYSTEMS 


There are many types of coordinate systems used in mission analysis applications. 
IPOST provides a number of systems to allow the user a maximum amount of 
analysis insight and flexibility. 


2.2.1 INERTIAL ECLIPTIC SYSTEM 


This system is used during heliocentric (sun-centered) interplanetary flight, 
although it is sometimes used as a fixed reference frame for the entire mission. 


Center of system 

Primary plane 
X - axis 


the sun, planetary bodies, or 
planetary satellites 
Earth ecliptic 

Earth vernal equinox direction 



Figure 2-4. The inertial ecliptic coordinate system. 

(The center can also be at any of the planetary bodies and 
satellites) 
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2.2.2 INERTIAL PLANETARY EQUATOR SYSTEM 

This system is used during flight near a planetary body. 


Center 

Primary Plane 
X-axis 


planetary body 
planetary equator 

rotation of planetary vernal equinox 
direction through right ascension and 
dedination angles of the planetary 
pole vector. 



Figure 2-5. The planet equatorial coordinate system. 


10 


& 


2.2.3 


This system is used for vehicles whose longitudinal axis or whose thrust axis is 
along the velocity vector. 

Center - S/C center of mass 

Primary Plane ■ orbital plane of S/C 

Primary (X) - cross product of S/C velocity U = V x 

W direction and S/C angular 
momentum direction 



— ► — » 
Irxvl 


U = V x W 


Figure 2-6. The UVW Coordinate System. 
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2.2.4 RTN SYSTEM 


This system is used for vehides whose longitudinal axis is along the local vertical, 
such as gravity gradient stabilized s/c. 


Center 

Primary Plane 
X-axis 


S/C center of mass 
orbital plane of S/C 
radius vector direction 



R = 


Irl 


N = 


rx v 


Irx vl 

— » -» -* 

T = NiK 


Figure 2-7. The RTN Coordinate System. 
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2.2.5 ] 
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BODY FRAME COORDINA TE SYSTEM 


The body frame is used in conjunction with other reference frames to orient the 
vehicle in celestial space, and to identify locations and orientations of vehicle 
components, such as the primary thrust vector and antenna boresights. 


Center - S/C center of mass 

x b axis - from center of mass through nose of S/C, 

z b axis - from center of mass through bottom of S/C, 

orthogonal to x 

y b axis - x b x z b 

P - pitch 

R - roll 

Y - yaw 



Figure 2-8. The body coordinate system. 
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2.2.6 



The B-plane coordinate frame is used for hyperbolic approach to, and departure 
from, a celestial body. It is often the most numerically stable system for describing 
planetary approach/departure conditions. 

HYPERB OLIC PATH 
OF SPA ^RAFT 



B = Impact Parameter (Vector from Planet Center to Aiming Point) 


0 = Orientation of B Relative to T 


S = Parallel to Incoming Asymptote 
— ► 

T = Parallel to Reference Plane (Ecliptic Unless Otherwise Specified) 
R =S xT 

Figure 2-9. B-Plane Coordinate System 
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2.2.7 CONE-CLOCK 


The cone-dock system is used for determining the orientation of vehide sensors and 
actuators. One application is to transform IPREP and QTOP thrust data into 
IPOST usable thrust acceleration profiles. 

z 



plant _L to tht oollptle 


<15 


Figure 2 - 10. Cone/Clock Coordinate System 
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3.0 INTERPLANETARY MISSION ANALYSTS 


Mission analysis is used to support all aspects of space exploration and applications 
missions (Figure 3-1). It covers both breadth and depth. In terms of breadth, 
mission analysis provides quantification, and optimization if needed, of desired 
performance measures for any major mission option. These options include 
planetary body encounters, means of propulsion, payload configurations, technology 
status, and mission goals. In terms of depth, mission analysis provides data used to 
support all levels of system/vehicle/ subsystem design and assessment. 

Mission analysis provides information related to geometry, dynamics, event 
times/conditions, and mass/energy relationships. With this information, a great 
number of activities are supported. In addition to those activities listed in Table 3- 
1, mission analysis can have strong interactions with structural design and stress 
analysis, human factors and life support, including habitation and bio-engineering, 
payload design, including science and instrumentation analysis, and life cycle 
cost/risk assessment. 


Mission Performance 

Assessing measures of effectiveness; Performance 
optimization; Trading off of mass-energy-time 


Optimization of propulsion performance; mode vs. 
phase, thrust vs. time, specific impulse, stage jettison, 
propellant and system mass 

Communications 

Analysis and optimization of coverage, link margins, 
positioning/pointing, blockage 

Power 

Assessing types of power sources and consumption 
patterns; Determining vehicle and solar array pointing 


Error analyses; Development of guidance and 
navigation requirements and capabilities; 
Design/analysis of sensor control systems 

Thermal 

Determining peak and cumulative thermal loads for 
active and passive systems 


Developing mission profiles and event sequences; 
Determining command and telemetry data and rates; 
Sizing of software and data processing 


Table 3-1. Representative Mission Analysis Applications 




The breadth and depth of mission analysis depends upon the development phase. 
Pre-phase A studies involve early feasibility and concept definition. Usually the 
starting point is a desire to depart Earth within some large time window and to 
arrive at a desired target body, e.g., Saturn or a libration point. Only a few mission 
objectives are specified, such as robotic vs. human, or surface probe vs. orbiter. 
Mission constraints, for example launch vehicle type or interplanetary propulsion 
mode, are often ill defined. Intermediate mission objectives, such as flyby of 
another object of interest, are also open for investigation as mission enhancements 
or to improve mission viability. For mission concept studies simple software tools 
are used because large numbers of concept variations must be assessed and 
compared. This is the realm of tools like IPREP and LPREP. The output of pre- 
Phase A studies is a preliminary definition of mission opportunities and constraints, 
with identification of potential difficulties, or "tall poles". 

Phase A studies examine mission feasibility, including preliminary definition of 
subsystem characteristics and performance. At this level, vehicle design is still very 
flexible, or "rubbeiy". System definition must address solutions to the tall poles 
identified in pre-Phase A studies. The need is for rapid, but reasonably accurate, 
vehicle trajectory information. Simulation models must be flexible and contain 
sufficient detail to provide reliable comparison of options. Robust performance 
optimization is needed because of widely varying initial conditions and 
environments. Phase A output is a set of feasible mission designs and performance 
comparisons, including preliminary subsystem design, usually with a recommended 
baseline. IPOST is the tool of choice, using simplified model options. 

Phase B studies define the mission and system in detail. Each mission phase from 
liftoff through mission completion is simulated using models which accurately 
reflect all relevant mission activities and environments. In order to provide 
information on vehicle and subsystem performance, specific phases often are 
optimized, such as launch vehicle ascent, as well as multiple phases, such as Earth 
departure to Mars landing. Both IPOST and POST are the tools of choice depending 
upon what mission phases are to be examined. The output of Phase B studies is a 
complete mission and system design from which hardware and software can be 
developed and tested. 

Phase C/D refers to conducting the operational and post-flight activities. Here the 
emphasis of mission analysis is on data base generation and validation, planning 
activities, such as maneuver design and implementation, anomaly resolution, and 
simulation support. IPOST (and POST) can be used to play "what if' games, such as 
contingency planning, troubleshooting, or mission redirection. 

A general flow of mission analysis activity using IPOST is shown in Figure 3-2. 

This flow can be modified and used for all development phases. It demonstrates (a) 
early concept evaluation using IPREP/LPREP, with possible verification using 
IPOST, (b) end-to-end simulation/optimization needed for feasibility studies, and (c) 
the more detailed design/assessment used for system definition, production, 
integration and operations. IPOST output is used by the mission analyst to support 
all elements of the system (such as Figure 3 - 1) in a highly interactive and iterative 
fashion. 
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The following sections discuss the major elements of mission analysis. These 
include early concept comparison characterized by mission opportunity definition, 
vehicle and trtyectory simulation, mission performance optimization, and how to set 
up IPOST to support the various mission analysis needs. 
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IPREP 

LPREP 


DEFINE PURPOSE 





Figure 3 - 2 IPOST ACTIVITY FLOW 









3.1 MISSION OP PORTUNITIES 

A definition of mission opportunity would be the time to launch a S/C such that the 
mission meets requirements as well as optimizes certain mission characteristics, 
such as total AV, or initial mass. In other words, what are the mission conditions 
that get the most performance for the least cost Occasionally, certain mission 
constraints prevent reaching an absolute optimum. Opportunities can still be found 
which meet mission constraints and represent a local optimum. There are also 
situations, due to the constraints imposed, where no feasible mission exists. 
Constraints include such things as boost capability vs. payload required, delta-v 
capability vs. delta-v required, launch times allowed due to planetary alignment, 
etc. There are a multitude of constraints which can be imposed as well as a 
multitude of variables to be optimized. 

Simple trajectory models, such as conics and patched conics, are used to do quick 
mission opportunity studies. Grid-search type optimization can be used to look at 
large time spans, in which many missions are run and contour plots can be created 
over time to view changes in the cost function relative to departure and arrival 
dates. 

To do pre-Phase A interplanetary mission opportunity studies, IPREP was created. 
Control of the mission design is left up to the user, such as, where the mission will 
fly, how to get there (high, low thrust, hybrid), what are the main concerns to 
optimize ( m i nimiz e delta-v, maximize payload), and regions where the search for 
mission opportunities should occur. Once the basic mission is designed, the user 
may use IPREP to find a "best” mission, the top number of "best" missions, or a 
large set of data displaying encounter time and each mission cost so graphical 
mission studies may be performed. After the user has identified desirable mission 
performance regions, EPOST may be used to generate more accurate trsyectoiy 
analyses and improve on this data. 

IPREP uses a simple conic propagator and Lambert method (Reference 4) for 
ballistic segments, and Chebytop polynomial representations (Reference 5) for low 
thrust segments, of interplanetary trajectories. Gravity assist is modelled as an 
instantaneous velocity change. The grid search is performed on a user input time 
window for any body encountered during the mission. The best mission is 
calculated by the minimum cost, which is a user-selected weighting of initial 
and each maneuver AV , arrival velocity and departure velocity. IPREP can be used 
to generate the initial conditions for IPOST for interplanetary missions. By using 
user is able to generate more accurate initial conditions, thus saving 
IPOST computer and user analysis time. 

To perform pre-Phase A lunar mission opportunity studies, LPREP has been 
created. As an offshoot of IPREP , it m a in t ain s the capability to search large regions 
of time quickly to perform mission opportunity analyses, and to generate initial 
conditions for IPOST. 
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Since the earth-moon problem is a three body problem rather than the more typical 
two body problem of interplanetary formulations, LPREP uses patched conics and 
actual targeting with a Newton-Raphson technique. Lunar mission conditions are 
generated to the same fidelity as interplanetary conditions. As in IPREP the grid 
search is performed on encounter time Gaunch and arrival) and optimization is 
done on total mission A V or with a generalized cost function. 

A complete guide to usage of IPREP and LPREP may be found in Appendix B. 


21 



3 2 TRAJECTOR Y SIMULATION 

A trajectory simulation is described by characteristics of the vehicle and by how the 
vehide behaves over time (Table 3-2). It is important to define vehide related 
characteristics because they represent necessary conditions that must be met, as 
well as the limits of trajectory fidelity, or accuracy. 


VEHICLE 

MODEL 

Propulsion Type 

Nudear 

Solar 

Chemical 

Aerodynamics 

Lift (constant coeffitient or 

Drag tabular) 

Mass 

Variable (finite and impulsive thrust) 
Jettison 

Event Activities 

Launch 
Orbit insertion 
Midcourse correction 
Jettison 

Planet/body or other changes 

TRAJECTORY 

MODEL 

Event Criteria 

User selectable trigger parameter 
e.g., dosest approach 

Event Categories 

Primary 

Secondary 

Roving 

Propagation Methods 

Explidt (Conic, Onestep, Multiconic, 
Encke, Cowell, and time step size) 
Implidt (Hermite polynomials over 
time segments) 

Environmental Forces 

Central body aspheririty 

Other body gravitational disturbances 

Atmosphere 

Solar Pressure 


Table 3 -2 . Simulation Characteristics 


Vdhtde characteristics in a trajectory sense are associated with forces acting on the 
vehide. Thus, it is necessary to describe propulsive modes (if any), aerodynamic 
properties (if any), and changes which affect mass properties, such as staging or 
reconfiguration. 

Trajectories are organized by time ordered events. The events can represent 
in sta n ta n eous activities, such as initiation of a pitch rate or rocket ignition, and 
functional activities which are approximated as an instantaneous process, such as 
launch to escape or an orbit insertion maneuver. Events can also be user directed 
changes for simulation purposes, such as changing propagation methods or 
activating output displays. A phase is the time period between consecutive events. 
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At any point in time, the vehicle can be described by a seven element state, e.g., 3 
cartesian position components, 3 cartesian velocity components, and mass. The 
vehicle state must always obey the equations of motion which operate between 
contiguous events, that is, during phases. 

In general, there are two methods of computing a vehicle state between events, 
explicit and implicit. The explicit method requires an initial state, at completion of 
the starting event. The equations of motion are integrated, either analytically, 
numerically or both, such that the state at the beginning of the destination event is 
computed. States at intermediate times are computed by interrupting the trajectory 
propagation at the desired time(s). Historically, explicit methods are the most 
common process for predicting future states. 

The implicit method requires vehicle states on each end of a time interval, either a 
phase or a phase segment. Hie start and end states of a phase or phase segment 
are used to generate coefficients of a vector function, such as a third order 
polynomial for each state component. Intermediate states are obtained by 
evaluating the vector function (interpolation). Implicit methods are often used 
durin g post-processing of results generated by explicit methods. Certain 
optimization methods, such as collocation, make use of implicit simulation 
representations. 

3.2.1 EVENTS AND PHASES 

Events and phases for a mission simulation must be organized with the specific 
application in mind. For feasibility assessments, many rapid simulation runs may 
be needed, either through a parametric scanning process or using 
targeting/optimization for a number of cases. The simulation should have as few 
events/phases as possible to speed up the run time. Typically, impulsive AV 
maneuvers and instantaneous activity approximations are used. Examples of the 
latter include launch to escape (LAUNCH) and orbit insertion from a hyperbolic 
approach (ORBINS). 

For quick simulations, the criteria, (CRLTK), for triggering events should be mission 
time, as opposed to a non-time criteria parameter, such as closest approach, because 
the simulation process must compute time-to-go in an iterative fashion for each 
propagation step. 

As the accuracy of mission design increases, more detailed models and a finer 
structure of events/phases are needed. Formerly impulsive maneuvers may become 
finite burns with an event to initiate thrusting and an event to terminate thrusting. 
For even greater accuracy, the thrust period may be divided into thrusting and 
coasting segments. Thrusting segments may be divided further into segments 
which alter the thrust profile, such as changing the throttle levels or thrust vector 
orientation controls. 

There is also the option to introduce different force models (see 3.2.2) in different 
propagation algorithms (see 3.2.3) and to alter step sizes depending upon mission 
phases. Usually there is also a need to change coordinate frames and/or planetary 
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bodies for ease of user understanding or for propagation accuracy. All of these 
simulation alternatives are accomplished through c h a n ges in the event activity 
specification. 

In addition to specifying event criteria and activities, IPOST allows the user to 
specify the event category. The primary event category is the most widely used. 
It is always triggered whenever the respective event criteria is encountered. 
Primary events must occur in consecutive order. A roving event can occur any 
time after its initiation when its criteria are met. 

3.2.2 VEHICLE FORCES 

The forces applied to the S/C during flight have a significant effect on the actual 
trajectory. Different mission phases and desired simulation accuracies dictate 
which forces are important. Studies can be done by turning forces on and off to aid 
the decision-making process of which forces are necessary to achieve the desired 
trajectory accuracy. Table 3 - 3 lists the vehicle forces implemented in IPOST 
including what trajectory propagators canbe^use<L 
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ORIGINAL PA3E to 
OF POOR QUALITY 


VEHICLE FORCES 

IPOST RELATED 
VARIABLES 

FORCE 

MAGNITUDE* 

Central Body Gravity 

IPBODY, GMU 

F = - r 2 

Propulsion 

• Low Thrust 

• Nuclear or solar electric 

• High Thrust 

- Impulsive/event 

- General, tabular 

IFORCE (5), PRPDT, SPI 
THRUST, WPROP, WJETT 

_ 2 7) P ( r, t) 

glSP 

p SCO 

F = F (m, t) - AePa 

Gravitational Perturbations 

• Nonsphericity of central 
body 

* Noncentral bodies 

IFORCE (2), GJ2 
IFORCE (1), RSOI, TSOI 

F = F(J 2 ,p,?) 

-if 

i=l 

Radiation Pressure 

IFORCE (3), SCSFA “ 


Aerodynamics 
* Drag and Lift 

IFORCE (4), CD, CL 
IATMOS, ATMS, SCDRAG, 
ALT, PRES, TPT, ITABSZ 

F=|pv2C L A 

* Atmospheric Density 

- Exponential Model 

- Pressure/Temp. vs. alt. 


F = v 2 C D A 

P = PO e - 
p = p (p, T) 

| * See Analytic Manual for detailed force model descriptions 


Table 3-3. Vehicle Forces 


A force which is always active is the central body gravitational acceleration. During 
interplanetary flight the vehicle passes from one central body influence to another. 
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A major vehicle force is propulsion. During early concept feasibility, propulsive 
thrust is often modelled as an instantaneous velocity change, or impulsive AV 
(F = oo). The simplicity of AV allows for rapid simulations which are needed in the 
parametric investigations of concept definition. As the vehicle and mission design 
progresses, higher fidelity propulsion models are needed, in particular finite thrust. 
Nuclear and solar electric propulsion provide veiy low thrust, on the order of .Olg’s 
or less. Low thrust propulsion must be activated for long mission times to 
accumulate sufficient orbit energy changes. High thrust propulsion typically 
operates over short durations. These include chemical systems, such as pressure 
regulated and blowdown. A moderate thrust propulsion technology is currently 
being advocated, namely nuclear thermal systems. All of these propulsion systems 
can be modelled in IPOST including user specified models, such as through tabular 
specification. 

The most often utilized non-central body forces in interplanetary space flight are the 
gravitational forces of other planetary bodies. These are usually perturbations until 
the S/C approaches the sphere of influence (SOI). In this case the primary and 
secondary body forces are the same order of magnitude and the secondary forces can 
no longer be considered a mere perturbation, but also a major force. In 
interplanetary flight, the forces of other planets can have an effect on the shape of 
the trajectory. When in orbit about Earth, particularly for GEO S/C, the Sun and 
moon can have a significant effect on the orbit. This is the same for other planets 
and their satellites. While, from a mission opportunity standpoint, these effects 
may be small and less important, they need to be considered for highly accurate 
trajectories to be obtained. 

When in orbit about planetary bodies, aspheridty of the gravitational field can also 
have a significant effect on the S/C trsyectory. IPOST models the predominant 
aspherical term, J2. 

Solar radiation exerts small forces on a vehicle with trajectory effects that 
accumulate over time. Radiation pressure effects are amplified for long mission 
durations (years), inner planet missions, and high area to mass ratio vehicles such 
as those having very large solar arrays. 

For interplanetary missions, aerodynamic forces have limited, but important 
applications. Drag due to the atmosphere of the central body can be a significant 
force if the S/C altitude is low enough. This is the case for low altitude orbits, and 
for ascent and entry trajectories* Another aerodynamic consideration, involving 
both lift and drag, is for aerobraking techniques which can be used for orbit capture 
and orbit c h a n ges. An example of the latter is an aeroassist planetary encounter 
based upon the hypersonic wave-rider concept. 

IPOST models the atmospheric interaction of lift and drag. Simple and tabular 
models for lift and drag coefficients are available. Exponential and 
pressure/temperature based atmospheric density models are also available. 
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Not all forces are active in each mission phase. In constructing the simulation event 
sequence, it is important to select or deselect the appropriate forces and force 
models. This is not only computationally efficient, but also minimizes errors caused 
by inadvertent force perturbations. 
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3.2.3 EXPLICIT PROPAGATION 


The method of explicit propagation is specified as part of the event activities. Once 
selected, it will be used for the simulation until changed by a subsequent event. The 
type of study the user wants to perform will usually define the propagation method 
needed to achieve the level of accuracy required. There are 5 propagation methods 
in IPO ST (Table 3-4) with different levels of capability as well as accuracy. The 
perturbing forces which are applicable for each propagator were described in Table 
3-3. 


Propagation Model 

Assumptions/Method 



Two-body motion: central body and vehicle 

^ 



mmn 

MULTICONIC 


ESsmU 

ENCKE ; 

Perturbed two-body motion with step size, 
numerical integration or nm«11 perturbations 

EaaBl 

COWELL " 

Numerical integration of equations of motion 
using any or all forces 


where: 

C = Central body gravity N = Nonsphericity of central body LT = Low thrust 
SB = Secondary body RP = Radiation Pressure HT = High Thrust 

L = Lift D » Drag 


Table 3-4. Explicit Trajectory Propagation Models 

The simplest propagation method is CONIC. CONIC propagation solves the two- 
body problem. Given an initial state, it will provide the final state at the new time. 
This propagator is fast, but does not allow for any forces other than central body 
gravity. When only two-body accuracy is desired (quick mission analysis with low 
accuracy), the conic propagator would be chosen. 
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On a higher fidelity level than the CONIC propagator, is the ONESTEP propagator. 
ONESTEP takes advantage of the quickness of the CONIC propagator, but is more 
accurate in that it allows for secondary body forces to be taken into account 
(Reference 6). No other perturbative forces are allowed; but it is a step up in 
accuracy from the CONIC propagator alone by providing an approximate solution to 
the three-body problem. ONESTEP is useful for propagation from low planetary 
orbit to outside the planet’s SOI into interplanetary space or vice versa. It does no 
integration but computes the trajectory in a single step, hence the name ONESTEP . 

MULTICONIC is more accurate than the previous two. The MULTICONIC 
propagator accounts for secondary body effects, and improves on ONESTEP's 
accuracy by allowing multiple steps and the addition of perturbative forces. It 
operates by performing a primary and secondary body conic trajectory each 
integration step, and adds in the perturbative accelerations halfway through the 
step. 

In using the ONESTEP and MULTICONIC propagators, the secondary body must 
be in orbit about the primary body. If this is not true, the mission setup will be 
inconsistent with the algorithm assumptions and the trajectory calculated will be 
useless. For example, with either of these two propagators, whether the S/C is 
within the SOI of a planet or is escaping the planet toward another body (other than 
a moon), the planet should be the secondary body with the Sun being the primary 
body. 

Numerical trajectory integration is performed with the COWELL propagator. The 
COWELL method directly integrates the equations of motion, including all 
perturbating forces. Even if the non-central body forces are on the same order of 
magnitude as the central-body forces, the COWELL propagator will deliver highly 
accurate trajectories. Although the COWELL propagator is the most accurate, it is 
also the slowest because it requires small integration steps. 

Another trajectory integration technique is the ENCKE propagator. ENCKE 
performs explicit integration like COWELL, but integrates only the perturbations 
and adds them to two-body CONIC motion. The assumptions behind this method 
are that the perturbations are very small when compared to the central body force 
and therefore change much more slowly. Thus the integration step size can be 
much larger than that for COWELL, up to 10 times larger during interplanetary 
flight. Whenever the perturbations are small flow thrust, solar radiation pressure) 
ENCKE would be a good choice. Whenever the perturbation forces become large 
though (e.g., continuous high thrust or near the planetary SOI) the assumptions are 
no longer valid and a switch to COWELL should be performed. 

It is not wise to mig propagator types within a simulation. The overall trajectory 
accuracy is only as good as its weakest link, or propagator in this case. One possible 
exception is to match similar type propagators, such as CONIC and ONESTEP or 
ENCKE and COWELL, for different phases. One must also be careful to be 
consistent in the step size adjustments. 
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3.2.4 


•rtl 


PLICIT SIMULATION 


Implicit simulation computes the vehicle state between two junction points, or 
nodes. This trajectory segment can be a portion of a simulation phase or the 
entire phase. At nodes between events, the vehicle state must be file same for each 
side, that is, just prior to and just after the node. At nodes that correspond to 
events, the states on each side may be different. This discontinuity would occur if 
instantaneous event activities were specified, such as mass jettison or an impulsive 
AV maneuver. 


The vehicle state is represented between nodes as Hermite, or cubic, polynomials. 
Each state component corresponds to a single Hermite polynomial. To compute the 
polynomial coefficients requires position, velocity, and acceleration (from evaluating 
the equations of motion) at each node. Figure 3-3 describes the methodology and 
compares it with explicit propagation. There is a crude correspondence between the 
integration step size of explicit propagation and the segment interval of implicit 
simulation. The latter is usually larger than the former for the same level of 
accuracy. 
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• Assume x\ and X^ 2 are known 


* Determine cubic polynomial for each component ofX ( x ), 1 1 £ x < 1 2 , 

For example, x = ao + aj.t + a 2 t 2 + a 3 t 3 

where a Q , a j , a 2 , a g are computed from boundary conditions .^2 

• Evaluate x(t) = ao + ajt + a 2 t 2 + a 3 t 3 

y(t) = bo + bit +b 2 t 2 + b3t 3 

Z (t) = fo + fit + f 2 t 2 + f3t 3 

m(t) = g 0 + g 1 t + g 2 t 2 + g 3 t 3 

Figure 3*3. Explicit and Implicit Simulation Methods 
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To initialize implicit simulation (and optimization) requires some a-priori 
knowledge, however poor. Initial guesses usually come from an explicit simulation. 
The easiest method is to use a mission opportunity tool, such as IPREP, to compute 
states at major events. IPOST allows two methods of 

determining undefined nodal states: Linear interpolation between defined states 
and Hermite interpolation from a previously defined state. Another initialization 
method is to use states from a previous implicit simulation/optimization run, such 
as when increasing from N segments per phase to M segments per phase (M > N). 

In general, it is better to start with as few segments per phase as possible, 
consistent with the initial guess accuracy, and then to build more accurate solutions 
as the simulation/optimization process progresses, as illustrated in Figure 3-4. 
Furthermore, the user needs to be aware of the dynamics involved in each phase. 
Some phases involve much more trajectory motion than others, which means there 
should be more Hermite polynomials, or segments, for the active phases then for the 
relatively quiescent or linear phases. These considerations will help control the 
region of stability and help prevent optimization divergence. 
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Figure 3-4. Implicit Simulation Accuracy 



3.3 PERFORMA NCE OPTIMIZATION 

Mission optimization (Table 3-5) is a critical part of the system design process. 
Quite often the function to be optimized is either total system mass or masB 
delivered to the final destination. However, the objective function can also be flight 
time, electrical power, or any system parameter which is important to be optimized. 
So m eti m es the objective function is called the system cost, in which case it is 
always minimized (a negative cost means the objective function is maximized). 


Optimization 

Techniques 

Maximize F ( u,y) 
Subject to u (y) 

ymin £y £ymax 


Master Problem and Subproblems 
Collocation 

None (Simulation or targeting only) 

Optimization 

Variables 

F = F( u, y) 


Final mass 
Sum of A Vs 

Initial mass in Low Earth Orbit 
Flight Time 

Control 

Variables 

(independent) 

u = u 1 ,u 2 

upper and lower 
bounds, weighting 


Initial maM 

Earth departure V. magnitude 
Thrust start time 
AV direction 
Encounter conditions 

Constraint 

Variables 

(dependent) 

y =yi.y2 

upper and lower 
bounds, weighting 

MSJLiiTf ra il 1 

Altitude 

Closest approach time 
BDT and BDR 
Inclination 
Velocity 


Table 3-5. Optimization Characteristics 

In order to optimize a function there must be some allowable degrees of freedom. 
Free parameters are called control or independent variables. Examples include 
initial mass and thrust management parameters (throttle setting and direction). 
Controls are usually given bounding values that reflect physical limitations, such as 
0 to 1 for throttle setting. 
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Another important optimization consideration is the presence of constraints. 
Constrained or dependent variables are parameters that must be met within 
bounds in order to have a feasible mission. Examples are minimum altitude (closest 
approach to a planet surface), rendezvous velocity, and peak dynamic pressure or 
acceleration. 

The optimization process must adjust the control parameters to jointly meet 
constraints and minimize cost. Optimization methods should be accurate, robust, 
and fast. Accuracy means find ing a global optimum solution, as opposed to a local 
optimum. Robust means findin g the optimal solution in the presence of widely 
varying conditions such as poor initial guesses of the control values (often the case) 
and non-linearities in cost function, constraints, trajectories, etc. Fast means 
rapidity in compute time and in user interaction/understanding. These goals are 
seldom achieved, particularly for complex missions. 

Among the methods of optimization that a mission analyst must choose are problem 
formulations and mathematical solution techniques. In IPOST the optimization 
problem can be formulated in an explicit and an implicit form. The solution 
technique uses NPSOL or PGA algorithms. 

Occasionally, no optimization is desired or needed. In this case, the analyst is 
interested in generating single pass trajectories or simulations, such as for 
sensitivity analyses. Alternately, the analyst may be looking for a mission that does 
the best job of meeting constraints (targeting only), and optimization is either 
premature or unnecessary. 

3.3.1 OPTIMIZATION TECHNIQUES 

There are two aspects to solving trajectory performance optimization problems: 
problem formulation and numerical solution. Problem formulation refers to the 
structure and organization of the optimization problem and its connection with the 
simulation process. This includes identification of control parameters, constraints 
and the objective function. The numerical solution refers to the algorithm which is 
applied in order to generate an optimal solution which meets all constraints. 

Classical optimization formulations are straightforward. The control parameters 
correspond directly to mission relevant parameters, such as thrust magnitude, that 
a mission designer can select within limits. Constraints are mission relatable 
variables which are associated with mission goals. The objective or cost function is 
the most dominant parameter that characterizes mission success or viability. This 
classical formulation is difficult to solve, except for very simple simulations, due to 
non-linearities in the kinematics and dynamics associated with the time-ordered 
simulation. Because of these instabilities, IPOST uses decomposition techniques. 

Decomposition breaks up the optimization problem into smaller, more tractable, 
sub-problems which are then unified or integrated by a master problem to achieve 
the overall solution. Decomposition (Reference 7) is applied to both explicit and 
implicit simulations. 
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Explicit optimization is similar to the classical approach in that all parameters are 
mission relatable quantities. Decomposition is applied by breaking up the 
simulation into major sections, or subproblems. Each subproblem represents a 
targeting problem, that is, meeting only constraints. The general subproblem is 
nonsquare, allowing the number of subproblem controls to be different than the 
number of subproblem constraints or targets. Each subproblem is solved in 
sequence with the beginning of one subproblem contingent upon the results of the 
previous subproblem. When the sequence of subproblems is completed, this 
represents a single simulation or mission pass. 

The master problem has its own set of controls and constraints and is responsible 
for optimizing the objective function. The non-linearities of optimization are 
m i nimi zed by (a) reducing the dimensionality (number of controls and constraints) 
for the subproblems and the master problem, as opposed to one single master 
problem, as in the classical approach, and (b) reducing the simulation time span 
through subproblems. This decomposition method has been successfully used for 
very non-linear simulations, for example JPL's PLATO/MOSES were used to design 
the Voyager and Galileo missions (Reference 6). Section 3.3.2 discusses explicit 
optimization in more detail, and Section 3.4.2 discusses EPOST related parameters. 


Implicit optimization also applies decomposition principles. In this case the 
simulation is divided into multiple segments per phase. The segment junctions are 
called nodes, with variable vehicle states on each side of the node. The controls are 
now the mission relatable free parameters plus all the vehicle states for all the 
nodes. The constraints are the classic mission constraints plus mathematically 
imposed restrictions to ensure continuity across nodes and to ensure that the 
equations of motion are met at the nodes and at the center of each segment. 

Clearly, the dimensionality of implicit optimization, or collocation, is quite large. In 
this case, the large dimensionality is offset by a reduction in non-linearity. Section 
3.3.3 discusses implicit optimization in more detail, and Section 3.4.2 discusses 
IPOST related parameters. 

The numerical solution of parameter optimization problems is as much art as 
science. IPOST uses primarily the NPSOL algorithm, although some versions allow 
the option to use PGA. 

The projected gradient algorithm (PGA) is an iterative technique designed to solve a 
general class of nonlinear progr amming problems. PGA employs cost-function and 
constraint-gradient information to replace the multidimensional optimization 
problem by an equivalent sequence of one-dimensional searches. In this manner, it 
solves a difficult multidimensional problem by solving a sequence of simpler 
problems. In general, at the initiation of the iteration sequence, PGA is primarily a 
constraint-satisfaction algorithm. As the iteration process proceeds, the emphasis 
changes from constraint satisfaction to cost-function reduction. PGA is a 
combination of Rosen's projection method for nonlinear programming and Davidon's 
variable metric method for unconstrained optimization. 

The NPSOL algorithm (Reference 8) uses major and minor iterations to solve the 
targeting and optimization problem. In the major iterations NPSOL seeks a 

significant decrease in the merit function along a direction of search p*. NPSOL 
defines the merit function as: 


F(5*) =X [li * (ci (i?)-si)]+|x [ri* (ri (i?) - si) 2 ] 
i z i 

where F(u*) is the objective function 

li is the vector of lagrange multipliers 
d is the i-th constraint function 

si is the set of slack variables used to handle inactive constraints and 
ri is the penalty vector for constraint violations. 

Therefore, during NPSOL operation, the problem becomes better targeted and more 
optimized in a simultaneous fashion. 
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In the minor iteration NPSOL seeks the search direction P for the major iteration 
by minimizing the quadratic programming problem: 


Q = p > +|p > TH^ ) 

subject to a set of constraints, where is the gradient of F(u*), and H is the quasi- 
Newton approximation to the Hessian of the merit function. 

Both PGA and NPSOL are gradient based algorithms. They require derivatives of 
the objective function and of the constraints with respect to the control variables. 
The former vector of derivatives is called the objective gradient and the latter 
matrix of derivatives is called the Jacobian, or sensitivity, matrix. In addition, 
NPSOL constructs a second derivative matrix or Hessian. 

One of the keys to optimization success is a well conditioned Jacobian matrix. 
Usually, this is formed by finite differencing. In IPOST, the finite difference control 
perturbations can be input by the user or computed by NPSOL. If NPSOL computes 
the perturbation size for each control parameter, the appropriate perturbation is 
that control value which produces the most accurate partial derivative for the 
control-constraint (or objective) combination, which is just above the numerical 
noise threshold. It can use either forward or central differencing techniques. Quite 
often about a third of the run is spent computing an accurate Jacobian because 
many simulation passes, or function evaluations, are needed to produce the finite 
difference perturbation. NPSOL also has the option of "verifying" an input 
perturbation size by constructing and comparing a finite differenced Jacobian. 

Another key to optimization guesses are "normalized" conditions. The ideal problem 
would have all values of the control parameters within one order of magnitude of 
each other. Similarly, all values of the constraints would be dose to each other, 
all non-zero elements of the Jacobian matrix would lie wi thin a one order of 
magnitude range. The no rm a liz ed conditions can be approached to some extent by 
adjusting the control parameter weights (WVU) and the constraint parameter 
weights (WVNLC). Section 3.4.2 discusses IPOST wei ghting techniques «nH 
options. 


3.3.2 


ELICIT OPTIMIZATION • MASTER/ST T BPROBT.BMS 


To make the targeting/optimization for missions more stable, IPOST allows 
decomposition of the trajectory into a master problem and subproblems . If the 
trajectory has several legs (planetary encounters, impulsive corrections, etc.), each 
leg can become a subproblem and targeted individually. 


For example, the user could choose a mission to launch from Earth park orbit, flyby 
Venus, do a deep space maneuver on the way to Mars and then enter park orbit at 
Mars, while m a xim i zin g the arrival payload (Figure 3 - 5). This problem could be 
set up as one entire master problem, the classical optimization approach, with 
approximately 10 control values, 6 constraints, and an optimization variable. 
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Figure 3 * 5. Explicit Optimization Decomposition Example 



Alternately, the problem can be decomposed with targeting performed on the 
subproblem level. The Earth to Venus trajectory would be solved with three 
controls and three targets, and the Venus to Mars trajectory would be solved by 
controlling the deep-space maneuver to target Mars. Above these two subproblems, 
on the master level, the user would control the encounter times as well as the flyby 
conditions at Venus to find the best conditions to maximize the arrival payload. 

A special application of this master/subproblem description is the use of ITOO 
(Interplanetary Targeting and Optimization Option) for grav ity assist missions, e.g., 
Voyager. ITOO uses analytic Jacobian partials and the ISTEP propagator. It is 
much faster than conventional methods. 

Which is the better approach, to do three simple problems or one complex problem? 
The three simpler problems are much easier to follow, but of more importance is the 
improved stability of the overall optimization process. For very simple missions, all 
of the targeting and optimization can be done on the master level; but when the 
missions become more complex, decomposition is needed to make the problem 
tractable. 


3.3.3 
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PLICIT OPTIMIZATION . 


COLLOCATION 


The technique of collocation (References 9 and 10 ) requires setting up the trajectory 
as an implicit simulation (see 3.2.4). Each simulation phase is divided into one or 
more segments. Nodes occur at segment boundaries, inclu ding events. In 
collocation, the states on each side of each node, that is, the pre- and post-node 
states, are allowed to vary. These free states become additional control variables. 

Additional constraints must also be introduced because the states are not totally 
free. For example, the state on each side of an event node is related by the event 
activity. Thus, an additional constraint is any state discontinuity arising from the 
activity, such as an impulsive AV. 


Another set of constraints is introduced because the state on each side of a node 
which is internal to a simulation phase, or between events, must be the same. In 
some collocation formulations, the zero state differences for these internal nodes 
are introduced as constraints. However in IPOST, the state on each side of an 
internal node is explicitly set equal to the state on the other side. This has the 
effect of reducing the number of controls and constraints such that only a single 
state at each internal node become control variables. 


One other IPOST refinement eliminates the state differences of the first event node 
by explicitly setting the post-event state equal to the pre-event state plus any state 
discontinuity associated with activities in the first event. In theory, all double 
states can be reduced to single states by this procedure, but this remains a future 
modification. 


There are other constraints introduced by collocation. These co ns traints involve 
forcing the Hermite polynomials to meet the equations of motion at the mid-point (tc 

) of each segment. The constraints, called defects, (cf) are formulated as: 
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derivatives 


(f = X* Q?, tc) - X*hP = o 


x = [ x , *y , *z , x , y , z , m ], the seven element vehicle state 


where tc ) = equations of motion evaluated at t*. using the state values (x) 
interpolated from the Hermite polynomials. 

= analytic derivations of Hermite polynomials, e.g., 

if y = b o + b 1 1 + b 2 1 2 + b3t 3 

then yjjp = bi + 2b2 + 3b3t 2 


The net result of collocation is that for any reasonable simulation/optimization 
problem there are hundreds more control and constraint parameters than in an 
explicit simulation/optimization(see Figure 3-6). In the IPOST formulation, an 
oqunl n umb er of controls and constraints are added via collocation. For example, if 
a classical problem has 10 controls and 6 constraints, and the simulation has 5 
ptiftBAn with 4 segments per phase, then the total number of controls and constraints 
are : 


Number of controls = 10 + nc * 7 = 192 
Number of constraints = 6 + nc * 7 = 188 
where nc = (number of segments per phase + 1) (number of phases) + 1 
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controls: x = state 

■* 

constraints: E = event state (discontinuity 
3 = state derivative defect 



internal internal event 

node node node 


event internal internal internal event 

node node node node node 


igure 3 - 6. Controls and Constraints Introduced by Collocation 


Although the dimensionality of the optimization problem has been greatly 
increased, there are several compensating features to collocation which ultimately 
expand its region of convergence. The first benefit is that each segment is relatively 
independent of the others because of the free states at each node. Thus, a change 
early in the mission is not amplified into extremely large, and unpredictable 
changes later in the mission, as in explicit simulation/optimization. 

Another collocation benefit is that most derivatives are analytic, and therefore 
rapidly, and accurately, computed. The few derivatives that must be done by finite 
differencing are of the "local" variety. That is, there is very little time mapping, at 
most over a single segment, which makes the finite differenced partials relatively 
stable and well-behaved. 
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Finally, the total Jacobian, although quite large, i6 a sparsely populated banded 
matrix. There are many techniques for efficient generation and manipulation of 
sparse ma trices. In IPOST, the zero elements of the Jacobian are identified and 
pre-set, t-h |1R eliminating some of the finite differencing in NPSOL for Jacobian 
elements that should be zero. 
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3.4 IPOST PROGRAM SETU P AND OPERATIONS 

In formulating a simulation/optimization problem for an interplanetary mission, it 
is usually best to construct the simulation flow first, and then op timi zation This is 
because simulation is easier to visualize, and optimization relies on the simulation 
activities, sequences, and results. It is then a straightforward mapping into IPOST 
usage. 

Each simulation event (activity, criteria and category) corresponds to a nameli st 
($TRAJ) input block. Whatever is established in the event description will be 
applied in the following time interval, or phase, until the next simulation event is 
encountered. After the simulation is defined in terms of $TRAJ blocks, the 
targeting/optimization namelist ($TOP) can be constructed. 

$TOP will identify such options as master/subproblemB vs. collocation, control vs. 
independent parameters (type, initial value, bounds, weights), constraint or 
dependent parameters (type, desired value, bounds, weights), and objective function 
identification. 


If the mission is very complex, involving many events/activities and/or many control 
and constraint parameters, it is best to start as simple as possible, and build the 
problem solution in stages. This serves three purposes: (a) easier diagnosis of 
input and simulation emirs, (b) developing a better feel for the complexity and 
sensitivities of simulation - optimization coupling, (c) b uil ding an orderly 

PJ°£T? 88ion of parameter values, such as bounds and weights, to maintain solution 
stability. 


For explicit simulation/optimization this means working the major mission phases 
nrst, with simple and few subproblems. For implicit simulation/optimization it also 
means starting with a portion of the mission, as well as using only a few segments 
per pb^se initially. Another option is to use a mission opportunity solution, such as 
irom IPREP, as a starting point for a simple collocation problem. Then transition to 
the more complex explicit optimization formulation. Although tempting, it is 
usually disastrous to start immediately with the full blown complex problem unless 
one has successfully completed a very similar problem. 


™J«iP* is via aamefists: $TOP, $TRAJ and $TAB. $TOP must be input 
first. $TRAJ must follow $TOP , and there must be one $TRAJ for each event. 
$TAB is used to input tabular data for thrust, aerodynamic coefficients and 
spacecraft mass. There may be one or more $TAB, following each $TRAJ. 
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Input and output units are metric. Unless otherwise specified, the following units 


Time 

days 

Mass 

kg 

Position 

1cm 

Velocity 

km/sec 

Area 

met 2 

Volume 

met 3 

Density 

kgtaet 3 

Force 

Newtons ( = kg-met/sec 2 ) 

Pressure 

Newtons/met 2 

Power 

watts ( = kg-met^sec 3 ) 


3.4.1 TRAJECTORY SIMULATI ON AND STRAJ 

IPOST trajectory simulation is based on event cycling. An event is a point in the 
trajectory where certain user input criteria (CRTTR) are met. The event CRlTK can 
be time-based or non-time-based. The user inputs an EVENT sequence number 
(e.g., EVENT = 5), a CRTTR variable (e.g., CRTTR = TDURF or time of duration 
since the last event) and an event criteria VALUE (e.g., VALUE = 20.d0 or twenty 
days since the last event). This event will not begin until the variable in CRTTR has 
achieved the state of VALUE. 

Care must be taken in choosing variables and values of CRTTR which are not time- 
based. If non-time-based CRTTRs are chosen, such as when a specified radius is 
reached, the user must be sure that the trajectory will pass through the specified 
value. Otherwise the event will never be reached and an event cycling error will 
occur. 

Once an event has been reached, any number of simulation changes can be made via 
$TRAJ inputs. Table 3-6 describes some commonly used inputs. These inputs 
include changin g the primary and secondary bodies, switching out propagators or 
trajectory acceleration forces, requesting impulsive maneuvers, or changing the 
input, output frames. For a more complete list of trajectory simulation inputs, see 
Section 5.2. 

Input state parameters are used only if the simulation state is to be initialized, e.g., 
launch, or reinitialized, e.g., return to Earth after a planetary stay. 

The maneuver mode (MANTYP) identifies certain instantaneous activities. These 
activity approximations include LAUNCH (impulsive injection into a hyperbolic 
trajectory), IMPULS (impulsive midcourse maneuver AV) and ORBINS (impulsive 
orbit insertion from a hyperbolic trajectory). 
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For explicit simulation, the propagation mode (IPROP) and step size (DT) are 
important to specify. For implicit simulation, and collocation, the user specifies the 
number of segments per phase with the variable NSEGPH. 


event 

data 


fEVENT 
CRITR 
< VALUE 
IDBODY 
IJDFRAMO 


Event/phase number 
Criteria for initiation 
Criteria value 
Planetary body 
Fame of reference 


flNPUTX 

input J X ( ) 

state ] IEPOCH 

Idateo 


Mode of state input/update 
State vector 
Type of epoch time 
Calender date 


propagation 


flPROP 
J IPBODY( ) 
j IFORCE 

Idt 


Propagation mode 
Primary and secondary bodies 
Perturbing forces 
Propagation step size 


maneuver/ 

guidance 


fMANTYP 

SPI 

SCMASS 
< THL 
WPROP 
PRPDAT 
vPITCH & YAW 


Maneuver mode 
Specific impulse 
Vehicle mass 
Throttle level capabilities 
Propellant maas 

Propulsion system characteristics 
Vehicle pointing coefficients 


input/ 

output 


[NAMLST 

PINC 

'PRINC 


Flag for next namelist input 
Printout interval 
PROFIL output interval 


Table 3 - 6. COMMONLY USED $TRAJ INPUTS 

of R * gs in I f ( ? R , CE ide ntify which forces are active for the current phase. 
A summary offeree related parameters was shown in Table 3 -3. For example 
finite bum segments are activated by IFORCE(5) which specifies the type of 
propulsive thrust (chemical/blowdown, electric propulsion, etc.) To deactivate 
thrust, set IFORCE(5) = 0 which is the a-priori default. WPROP specifies 
propellant mass, which is depleted as the vehicle thrusts and consumes propellant. 

Depending upon the type of thrust, other parameters must be specified in PRPDAT 

nSnSJSvS l0 ^ t 1 hrust 8 ° lar °J nuclear propulsion, require PRPDAT(9) through 
FKPDAT(15) which set values for thrust efficiency, housekeeping power, etc. Peak 
thrust is determined by the propulsion type through PRPDAT data. The thrust 
vector is assumed to be aligned along the vehicle x (roll) axis. 
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When thrusting, the user must specify throttle level (THL parameters), and 
direction, such as vehicle pitch (PITCH parameters) and yaw (YAW parameters). 
The throttle, pitch and yaw parameters can be represented by quadratic functions 
operating over the thrusting phase. For example, THLO is thrust bias (initial thrust 
at beginning of phase), THL1 is the thrust rate of increase, and THL2 is the thrust 
acceleration coefficient. These coefficients can also be used as optimization control 
parameters along with thrust on and off times. 

Spacecraft mass is a critical mission parameter and is often the objective function to 
be optimized. SCMASS in $TRAJ defines the mass. As propellant is consumed 
during thrusting phases, or mass is jettisoned, SCMASS diminishes accordingly. 
Two other mass related parameters are WPROP and WJETT. 

The propellant mass, WPROP, can be specified at any event. Thereafter, WPROP 
will diminish (as does SCMASS) whenever thrusting occurs. A variable initial 
WPROP can be specified by using the parameter, FRACI, where WPROP = FRACI * 
SCMASS. Alternately, a mass fraction table (FMASIT) can be input which then will 
override FRACI. 

Quite often a mission will jettison mass, such as a probe or spent fuel tanks or a 
complete stage. The form of jettison mass is dictated by the flag, JETT. 

For JETT = 0, the value of WJETT input at any event is ejected. 

For JETT = 1, and for negative values of JETT, a calculated value of mass is 

jettisoned at the current event. 

For JETT = 2, the jettison mass (defined internally as WJETTM) is the input FRAC 
value, or WJETTM = FRAC. Alternately, FRAC can be replaced by values defined 
in the table FMASST. WJETTM is calculated and saved until another event 
specifies JETT = 1. JETT = -2 will result in the calculation and jettison of mass at 
the current event. 

For JETT = 3, the calculated and saved jettison mass is 
WJETTM = WPROP * ( 1/FRACI - 1) 

Alternately, FRACI can be replaced by values defined by the table FMASIT. JETT 
= - 3 will result in the calculation and jettison of mass at the current event. Note 
that for JETT * 3, the parameter FRACI (and its related table FMASIT) is used 
differently, in particular, to initialize WPROP as described earlier. 

Atmospheric flight requires specification of the lift and drag coefficients, CL and 
CD, respectively, and the atmosphere model (IATMOS). The user must identify 
atmospheric characteristics. For example, the exponential density model (IATMOS 
= 1) requires base density (ATMS (1)), scale height (ATMS(6)), and minimum 
altitude (ATMS(7)). In all cases ALTATM(I) represents the maximum altitude of 
the sensible atmosphere for the Ith celestial body. 
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3.4.2 PERFORMA NCE OPTIMIZATION AND STOP 


The $TOP namelist contains inputs which describe the targeting and optimization 
process. Table 3 - 7 illustrates commonly used inputs. 


control or 

independent 

parameters 


f INDVR ( ) 
INDPHO 
U() 

INDXI ( ) 
INDPLBO 
INDPUBO 
WVUO 


parameter name 
phase 

initial value 
active parameter 
min value 
max value 
weighting 


target or 
constraint or 
dependent 
parameters 


f DEFVR ( ) 
DEPPH ( ) 
INDXDO 
DEPTLO 
DEFVLBO 
DEPVUB ( ) 
UWNLC( ) 


parameter name 
phase 

active parameter 
tolerance/weighting 
min value 
max value 
weighting 


Optimization 
Master problem 


rOPTVAR 

OPTPH 

Ifesn 

) SRCHM 
OPT 

MVDOTOP 


Optimization variable 
Optimization phase 
Last event of simulation/optimization 
Method of optimization 
Targeting/optimization mode 
Maximum iterations 


TABLE 3 - 7. COMMONLY USED $TOP INPUTS 


Control parameters (independent variables) are defined in namelist $TOP by the 
array INDVR. Table 3 - 8 identifies allowable control parameters. In theory, any 
real variable in namelist $TRAJ can be used. An example of usage would be 
INDVR(3) = 6HPITCH1, which specifies that the third control parameter is the body 
pitch rate. 

CHARACTER DEFINITION 


ALTIT 

C3 

CRITR 

DEC 

DEPSVLij 

DVX 

DVY 

DVZ 

INC 

PITCHO 

PITCH1 

PITCH2 

RA 

RAPOAP 

ROLLO 

ROLL1 

ROLL2 

SCMASS 

THLO 

THL1 

THL2 

VINFXO 

VINFY1 

VINFZ2 

vx 

VY 

VZ 

WJETT 

WPROP 

X 

Y 

YAWO 

YAW1 

YAW2 

Z 


Park orbit altitude 

Escape energy 

Event criteria parameter 

Escape asymptote declination 

Subproblem dependent parameter ij 

Delta velocity x-component 

Delta velocity y component 

Delta velocity z component 

Orbit inclination 

Initial pitch angle 

Pitch angle rate 

Pitch angle acceleration 

Escape asymptote right ascension 

Radius of apoapsis 

Initial roll angle 

Roll angle rate 

Roll angle acceleration 

Spacecraft mass 

Initial thrust throttle level 

Throttle level rate 

Throttle level acceleration 

Escape asymptote velocity X component 

Escape asymptote velocity Y component 

Escape asymptote velocity Z component 

Velocity X component 

Velocity Y component 

Velocity Z component 

Jettison mass 

Propellant mass 

Position X component 

Position Y component 

Initial yaw angle 

Yaw angle rate 

Yaw angle acceleration 

Position Z component 


TABLE 3 - 8. CHARACTER NAMES FOR CONTROL PARAMETERS 
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In Table 3-8, the parameter CRITR is used to represent another parameter. When 
INDVR = 5HCRITR, then the control parameter is the event criteria parameter 
that triggers the appropriate phase (INDPH). Event criteria parameters are 
defined in namelist $TRAJ and in Table 3-9. 


CHARACTER 


DEE 


■tiwmc 


ALTIT 

ANGAT1 

ANGAT2 

ANLONG 

ARGP 

BDRO 

BDRI 

BDTO 

BDTI 

BTHETA 

C3 

DEC 

DVMAG * 
DVSUM * 

ECCEN 

FPA 

HYPTA 

INC 

LAT 

LONG 

LONGP 

MEAAN 

PERIOD 

PITCH 

RA 

RADIUS 

RAPOAP 

ROLL 

RPERI 

SCMASS 

SMA 

SMACH 

SPEED 

TDURP 

TFP 

THL 


Altitude 

Angle of attack (in-plane) 

Angle of attack (out-of-plane) 

Longitude of ascending node 
Argument of periapsis 
B dot R (outgoing) 

B dot R (incoming) 

B dot T (outgoing) 

B dot T (incoming) 

B-plane theta angle 
Launch energy 

Declination of launch asymptote 
Delta velocity magnitude 
Summation of all DVMAGs over the 
course of the trajectory 
Eccentricity 

Relative flight path angle 

Hyperbolic turn angle between asymptotes 

Inclination 

Latitude 

Longitude 

Longitude of periapsis 

Mean anomaly 

Orbital period 

Pitch angle for the S/C 

Right ascension of launch asymptote 

Radius magnitude 

Radius of apoapsis 

Roll angle for the S/C 

Radius of periapsis 

Vehicle mass 

Semi-major axis of orbit 

Mach number 

Velocity magnitude 

Time since last event 

Time from periapsis 

Throttle level 


TABLE 3 - 9. CHARACTER NAMES FOR TARGETS/OPTIMIZATION/EVENT 

CRITERION/TABULAR VARIABLES 
*Not available for CRITR. 


TIME 

Time in days 

TIMRFi 

Time of activity i (initiated by DTIMR(i)) 

TPERI * 

Time of closest approach 

TRUAN 

True anomaly 

VINFX 

V-infinity X component 

VINFY 

V-infinity Y component 

VINFZ 

V-infinity Z component 

VMAG 

Velocity magnitude 

VX 

Inertial cartesian X velocity 

VY 

Inertial cartesian Y velocity 

VZ 

Inertial cartesian Z velocity 

WPROP 

Propellant mass 

X 

Inertial cartesian X position 

Y 

Inertial cartesian Y position 

YAW 

Yaw angle for the S/C 

Z 

Inertial cartesian Z position 


TABLE 3 - 9. CHARACTER NAMES FOR TARGETS/OPTIMIZATION/EVENT 

CRITERION/TABULAR VARIABLES (Continued) 

♦Not available for CR1TR. 

Hie initial values of control or independent parameters are specified in the U array. 
It is always best to specify values that are closest to the optimal values. For 
missions with some degree of uncertainty, a mission opportunity tool, such as 
IPREP, should be used to generate initial guesses. For more defined missions, 
results of previous simulations and optimal solutions should be used. 

Sometimes it is convenient to input all likely controls by specifying INDVR, INDPH, 
U,..., and then, for the problem at hand, selecting a subset for active controls. In 
this way, the user can select the most well-behaved controls for the problem or 
mission phase. The array INDXI allows this flexibility. 

INDPLB and INDPUB are lower and upper bounds for each control parameter. If 
the bounds are too loose, there is the potential for the solution to wander and 
possibly converge within a mathematically feasible, but physically undersireable, 
region. If the bounds are too tight, there may not be enough freedom to find the 
true global optimization. It is often best to keep the bounds reasonably small, and 
then relax those bounds that are reached by the solution. This requires a certain 
amount of user interaction and awareness. 

Another important control related parameter is the weighting (WVU). "Proper" 
weighting will improve the speed and robustness of the solution process. Improper 
weighting often leads to divergence. 
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Dependent, or target or constraints, parameters are described by a comparable set 
of variables to the control parameters. One key difference is that controls require 
an initial guess (U), whereas constraints are specified by their bounds (DEPVLB 
and DEPVUB). Another difference is that constraints permit usage of a normalized 
tolerance (DEPTL) which affects constraint error metrics. Is is recommended that 
the default value of DEPTL be used except where significant experience dictates 
otherwise. 

T he op timization process is characterized by OPTVAR (parameter to be optimized), 
OPTPH (phase of optimization parameter), FESN (final event number), SRCHM 
(optimization formulation type), OPT (maximize or minimize optimization 
parameter), and MXITOP (maximum number of iterations for performing selected 
optimization activities). 

IPOST allows great flexibility in selecting variables. 

Table 3-9 describes the allowable variables for target parameters (DEPVR), 
optimization parameter (OPTVAR), event criteria (CRITR), and table input 
independent parameters). Variables denoted by an asterisk are exceptions to 
CRITR availability. 

The care and feeding of NPSOL is very important to the success of finding an 
optimal solution. One important aspect is the scaling or weighting of the control 
and constraint variables, WVU and WVNLC, respectively. Bad scaling of the 
problem can exacerbate the ill-conditioning and non-linearity of the problem. For 
independent variables experience shows that the best weighting is the magnitude of 
the variable. Since weighting is used in IPOST as a denominator, NPSOL will then 
operate on variables that are all roughly near unity. For dependent variables the 
weighting should be chosen so that the range of magnitudes in the sensitivity 
(Jacobian) matrix is minimized. For collocation, IPOST provides a set of default 
weightings for the supplementary controls (node states) and supplementary 
constraints (defects and event discontinuities). These defaults are described in 
more detail in Section 5.1 within the definition of WVNLC and WVU. 

One important option is the ability to rescale independent and dependent variables 
Initially, the user should input values of WVNLC and WVU. Good first values of 
WVNLC and WVU are the magnitudes of the dependent and independent variables, 
respectively. For the venturesome user, initial WVNLC and WVU values can be 
overridden by input flag (IRSCL). A robust option is to automatically renormalize 
the Jacobian after it is calculated by rescaling the dependent parameters (IRSCL = 

2 and MXTTOP(2) = 0). For example, MXITOP = 0,10, would mean that IPOST 
would first evaluate the trajectory and its jacobian, rescale the weightings for the 
problem, and run the problem for 10 iterations. In the NPSOL iteration output 
check COND T for improvement after rescaling. See variable IRSCL for rescaling 
options. 
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3.4.2. 1 


PLICIT OPTIMIZATION ■ M ASTER/SUBPROBLEMS 


For explicit optimization, the problem can be decomposed or partitioned into a 
number of subproblems which are tied together by a master problem. The 
subproblems are defined such that the beginning of one is the end of the previous 
subproblem. Each subproblem represents a targeting and/or optimization problem. 

The parameters which define the master problem are the same as those in the 
previous section. For subproblem definition, there is an analogous set of parameters 
as illustrated in Table 3 - 10, with some important differences. The independent 
parameter active ID (INDSXI) specifies the subproblem for that control. For 
examp le, INDSXI(5) = 2 means that the fifth parameter corresponds to a control 
parameter for the second subproblem. 


Variable 

Master 

Subproblem 

Independent Parameters 



Name 

INDVR () 

INDSVRO 

Phase 

INDPHO 

INDDSPH ( ) 

Initial value 

U() 

USUBO 

Value lower boundary 

INDPLB ( ) 

INDSLB ( ) 

Value upper boundary 

INDPUBO 

INDSUBO 

Active ID 

INDXI ( ) 

INDXSIO 

Weighting 

wvuo 

wvuso 

Perturbation size 

PERT ( ) 

PERTSB ( ) 

Dependent Parameters 



Name 

DEPVRO 

DEPSVR ( ) 

Phase 

DEPPH ( ) 

DEPSPH ( ) 

Value upper bound 

DEFVUBO 

DEPSUB ( ) 

Value lower bound 

DEPVLB ( ) 

DEPSLB ( ) 

Value target 


DEPSVL ( ) 

Active ID 

INDXDO 

INDXSDO 

Tolerance 

DEPTLO 

DEPSTL ( ) 

Weighting variables 

WVNLCO 

WVSNLC ( ) 

Weighting errors 


WGTSO 

Objective Parameters 



Name 

OPTVARO 

OPTSVRO 

Objective magnitude 

OPTO 

OPTSO 

Phase 

OPTPHO 

OPTSPH ( ) 

Weighting 

WOPTO 

WOPTS ( ) 

Problem solution method 

SRCHM 

MODELT ( ) = SUBOPT 

Number of iternations 

MXITOP 

MXITAR ( ) 


Table 3 - 10. Master and Subproblem Parameters 
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There is &n option to use & subproblem target variable as a master problem control. 
For example, if INDVR = 8HDEPSVL02, then the master level control parameter is 
the L-th subproblem second dependent parameter. As with subproblem controls, 
the index INDSXSD specifies the subproblem for a given dependent parameter or 
target. For example, INDSRSD (3) = 1 means that the third dependent parameter 
is a target parameter for the first subproblem. 

When using decomposition, it is important to note that each master level simulation 
pass requires that all the subproblems be targeted/optimized in sequence. This 
means that finite difference construction of the Jacobian matrix and objective 
gradient requires each perturbed pass to successfully perform subproblem 
targeting/optimization. This can be a very time consuming process for complex 
master-subproblem formulations. 

Dependent parameter weighting for subproblems is different than master problem 
weighting. For subproblems, target weighting depends upon the solution method 
For Newton-Raphson (MODELT = ’NRAPH’), which is the recommended and 
default method, the target weights are the target values (DEPSVR), unless the 
target value is zero (then the weight is unity). For specific targeting (MODELT = 
TARG1S'), target weighting is specified by WVSNLC. 

There are other properties associated with the Newton-Raphson method of solving 
subproblems. It is not necessary to have an equal number of controls as target 
parameters. IPOST solves the non-square targeting condition with a Singular 
Valued Decomposition algorithm. Whereas a master problem dependent parameter 
value 18 specified by an upper and lower bound, a subproblem dependent parameter 
is specified by a single value (DEPSVL). There is also a master-subproblem 
difference in usage of the dependent parameter tolerance. Whereas the master 
problem tolerance (DEPTL) is a normalized numerical tolerance that is not typically 
changed from its default value, the subproblem target tolerance (DEPSTL) for 
Newton-Raphson is an important user input, and serves as a bound on the desired 
target value. The value of DEPSTL affects both the speed and accuracy of 
subproblem convexgence. An important guideline is to set the DEPSTL values for 
earlier subproblems tighter than for later subproblems to reduce numerical error 
particularly at the master problem level. Inconsistent DEPSTL’s will lead to 
erroneous solutions or divergence. 

The use of analytic partials (ISTM = ‘ANAL’ in $TOP) will speed up run time and 
improve convergence stability in most cases. Analytic partials are available only for 
the Come, Onestep, and Multiconic propagators, and available only for 
Interplanetary Targeting and Optimization Option (SRCHM = ‘ITOO’) or for Master 
problem only or for Sub-problems only. 
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For subproblem optimization each subproblem is constructed with parameters 
roughly equivalent to the method used in the master level. Subproblem 
optimization is selected by setting MODELT = 6HSUBOPT for the applicable 
problem. 

Subproblem optimization usually requires much more time than Newton-Raphson 
targeting, and there is no guarantee that the optimum subproblem is consistent 
with the overall master problem optimum. 

If the user intends to use the NPSOL41 algorithm to target the subproblem, but not 
optimize any variable simultaneously, the OPTSUP = 6HTARGET option should be 
used. This creates an objective which contains the magnitude of the errors of the 
constraint functions divided by the constraint target error weights, WGTS. In this 
way, the objective and constraints will produce a direction of search that is 
consistent for both sets and result in faster and more robust convergence. 

All guidelines in selecting variable weights for the master problem should be used in 
the selection of subproblem variable weights. Similarly, all the variables normally 
contained within the NPINPUT file for master level optimization should be 
reproduced in a dedicated NPINPUT01, NPINPUT02, etc. file for each subproblem. 
The use of these variables is unchanged from the master level, and can be tailored 
for each subproblem. 

The user may select a master problem control variable as the subproblem dependent 
variable in much the same manner as discussed earlier in subproblem targeting. 
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3.4J2.2 


101 


PLICIT OPTIMIZATION ■ COT.T.OCATTON 


The collocation method requires breaking up each phase into one or more segments. 
Extra degrees of freedom (controls) are introduced by adding the vehicle state vector 
at each segment boundary (node). The extra controls are mitigated by introducing 
extra constraints in the form of pre and post node state connectivity and by reducing 
the defects to zero (that is, forcing the state derivatives at segment centers to match 
the equations of motion). 


For implicit optimization, or collocation, SRCHM must be set to 'COLLOC'. The 
independent (INDVR) and dependent (DEPVR) related variables must be 
augmented to include all nodal state and defect related parameters. The following 
checklist may be helpful in setting up a targeting and optimization problem using 
collocation. 


First, set up the simulation with event data in $TRAJ. This defines how the 
implicit simulation is initialized, and how each phase is modelled. This discussion 
includes only those parameters directly related to collocation, and assumes all other 
mission data is already established and transformed into input. Some IPOST 
features are restricted by collocation, such as not allowing non-time CRITR’s or 
LAUNCH or ORBINS (MANTYP) maneuvers. 

1. For each $TRAJ, which corresponds to the be ginning of a phase, set 
NSEGPH equal to the number of segments per phase. It is usually best to 
start with as few segments as possible. It is also wise to make sure 
NSEGPH is consistent across the simulation by having those phases with 
high trajectory curvature use more segments than those phases with amall 
curvature. One should avoid having zero length phases, such as not setting 
TDURP = 0. 

2. For each $TRAJ, set NSGPHO, which is the number of input segment 
states for that phase. The default (NSGPHO = 0) means all internal node 
states are interpolated based upon the initial and final phase states. Setting 
2 < NSGPHO £ NSEGPH means that the user must input the state values for 
those internal nodes specified by NSGPHO (see item 6 for the order of the 
state input), and the post-event state. 

3. There are two options which allow overriding the default control «nH 
constraint weightings. These are not recommended except for the 
experienced user. NSGPWD and NSGPWI are parameters in each $TRAJ 
which allow the user to override the constraint weightings (WVNLC) and 
control weightings (WVU) for this phase. The default values of NSGPWD and 
NSGPWI automatically set them both equal to NSGPHO. See item 8 for the 
order of weightings input. 

Once the simulation has been set up, the optimization parameters in $TOP must be 
included. In addition to the standard parameters, such as OPTVAR and OPTPH, 
the following collocation specific parameters should be input. 

4. SRCHM = COLLOC' 
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5. After the mission control parameters have been established, the 
collocation related control parameters (U and WVU) must be input. The 
order in which control parameters are input is also the order in which all 
independent parameters are displayed as output. In general, the order of 
input for U (initial control values) and WVU (control parameter weights) is: 
all mission controls, all pre-event states (7 components each in local body, 
IPBODY (1), cartesian form) starting with the second event and ending with 
the final event, all post-event states and internal node states in time order. 
Table 3-11 illustrates the input order for a mission that has four mission 
controls, 3 phases (4 events), and 2-2-3 segments per phase. Note that the 
number of input U's and WVU's may be affected by NSGPHO and NSGPWI, 
respectively, that are in the $TRAJ inputs. When using NSGPWI to override 
default weights, a zero value of WVU will reinstitute the default value. The 
seven components of the state are the three components of position, in local 
body, IPBODY (1), ecliptic cartesian form, three components of velocity, in 
local body ecliptic cartesian form, and the spacecraft mass, SCMASS, in 
kilograms. The user can override the default and specify which body, with 
icbody, instead of local body, IPBODY(l), is the reference body for the 
collocation U's. 

6. It is not necessary to input other collocation control related parameters. 
For example, the collocation control parameter upper bounds (INDPUB) and 
lower bounds (INDPLB) are automatically set to +10E10 and -10E10, 
respectively. For some problems, these bounding values may be too loose. 

7. After the mission constraint parameters have been established, the 
collocation related constraint weights (WVNLC) must be input. In general, 
the order of input for WVNLC is: all mission constraints, all event state 
discontinuities, all segment defect states. This order of WVNLC parameter 
input is also the output order of all dependent parameters in IPOST display. 
Again, Table 3-11 shows an example of the input order. The number of 
input WVNLC's may be affected by NSGPWD which is in the $TRAJ input. 

8. It is not necessary to input other collocation related parameters. For 
example, the collocation constraint parameter upper bounds (DEPVUB) and 
lower bounds (DEPVLB) are all set to zero because these constraints must be 
met exactly. 

9. The user may select the method of state interpolation for those collocation 
states which are undefined (NSGPHO < NSEGPH) by inputting COINTP. 
The default is linear interpolation. 

10. As with explicit optimization, an important option is to rescale the 
constraints and/or weights using IRSCL. If rescaling is desired, it is 
recommended that IRSCL = 2 because this method of rescaling appears the 
most effective. 
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With collocation, the solution is very sensitive to control and constraint weightings, 
and with control parameter bounds. This applies to the mission parameters, as weil 
as those parameters introduced by the collocation formulation. A rule of thumb is to 
begin as simple as possible, understand the sensitivities and quirks, and then build 
up to the desired solutions. As the problem/solution evolve, so must all tolerances, 
accuracy parameters, etc. in a coherent fashion. When to rescale weights, which 
control/constraint combinations to select, and event/segment adjustments will also 
evolve. 
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Phase 1 ^ ^ — Phase 2 ^ ^ Phase 3 ► 

Seg A I Seg B I Seg A ^ Seg B I Seg A Iseg B I Seg C 


Controls: Earth AV, Venus fly. 
Constraints: Mars orbit R and 

P 


Earth 
(Event 1) 


Venus 
(Event 2) 


DSM 
(Event 3) 


Mars 
(Event 4) 


CONTROL PARAMETERS 


U, WVU, INDPLB, INDPUB, 
INDPTL 


Parameter Parameter Phase 
Name Number 




PARAMETERS: 


WVNLC, DEFVLB, DEPVUB, DEPTL 



Parameter Parameters) 
Number 



Earth DVY 


Earth DVZ 


Venus 



5-11 

pre Event 2 
state 

12-18 

pre Event 3 
state 

19-25 

pre Event 4 
state 

26-32 

post Event 1 
state 

33-39 

post Seg. 1A 
state 

40-46 

post Event 2 
state 

47-53 

post Seg. 2A 
state 

54-60 

post Event 3 
state 

61-67 

post Seg. 3A 
state 

68-74 

post Seg. 3B 
state 

75-81 

post Event 4 
state 



Mars i 


Event 1 state 
discont. 


Event 2 state 
discont. 


Event 3 state 
discont. 


Event 4 state 
discont. 


Defect Segment 1A 


Defect Segment IB 


Defect Segment 2A 2 


Defect Segment 2B 2 


Defect Segment 3A 3 


Defect Segment 3B 3 


Defect Segment 3C 3 


Phase 

Number 


3 


3 


0/1 




Table 3 - 11. Example Collocation Input/Output Parameter Order 
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3 - 4 -2- 3 MULTIPLE PROBLEM R UN CAPABILITY 

IPOST has the capability to do several problem runs consecutively. After the first 
problem, the user may change some of the $TOP parameters and then run a second 
problem. The user requests a second problem by inputting NAMLST = TOP', at the 
end of the first problem. IPOST will run the second problem if NPSOL terminates 
with INFORM = 0, 1, 4 or 6 or if SRCHM = ’NONE'j otherwise, the first problem is 
considered a failure and the job is terminated. The user may change any $TOP 
variable, except the following: 

Depph Indph IndrH 

Depsph Indsph Indii 

Depsvr Indsvr 

Depvr Indvr 

The user can run a third or a fourth problem, after the second problem, subject to 
the same restrictions. 
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3.4.3 PROCESSING 

As with all programs of this size and complexity, there will be blatant and subtle 
discrepancies. The more complex the problems become, the more difficult they are 
to solve, and there becomes more room for error in problem setup, simulation, 
targeting, optimization, etc. 

There are generally two types of errors which occur, external and internal. External 
errors are user errors which cause program malfunctions. These are generally 
keystroke errors, model selection errors, improper data vs. initialization, etc. These 
are usuall y found quickly and easily. Input error diagnostics are in Section 6.4, 
Error Messages. 

Internal errors are generally more difficult to find and to correct. Many times these 
can also be caused by input, and modelling, but in a more subtle way. An example 
of this type of problem would be that the subproblems will not successfully target. 

It could be that the parameters chosen as controls, or their allowable range of 
values, are inadequate to reach the specified targets. If the user has chosen forward 
differencing to calculate the sensitivities, possibly central differencing should be 
tried. Possibly, the size of the perturbations used in calculating the sensitivities 
was too large or too small. Perhaps control or constraint weightings were not 
chosen properly so that the Jacobian matrix elements vary by many orders of 
magnitude. As a general rule of thumb, weights should be chosen to normalize the 
control parameters, Jacobian matrix, and constraint parameters, in that order. 

If the program has run successfully, and a solution has been given, but NPSOL has 
delivered that infamous message "EXIT NPSOL - Current point can not be 
improved upon," another question comes up. How can the user find a global 
solution rather than a possible local minimum? How can this solution be improved? 
Again, there are many possibilities and the user must use some common sense and 
knowledge in choosing the proper path to take. Checking the output and seeing 
what the program tried to do is a good starting point. Were there any controls 
constantly being pushed to one of their bounds? Perhaps they need to be loosened 
up some, or tightened up. Perhaps the weightings or perturbations are limiting the 
step sizes too much. Is the program taking giant leaps after each iteration? Again, 
perhaps the weightings or perturbations are allowing steps that are too large due to 
noise or non-linearities in the sensitivities. Check the objective gradient and 
Jacobian matrix. Do the elements seem to be near the same order of magnitude, or 
are there large differences? All of these things and more can be modified to improve 
solutions. 

Past use of the NPSOL algorithm has shown that the problem is frequently badly 
scaled and that NPSOL terminates with an INFORM=6, 'Current point could not be 
improved upon' message. The NPSOL user's guide states that this can happen if: 

1. Overly stringent accuracy has been requested, i.e.. Optimality 

Tolerance is too small. 

2. Routines OBJFUN or CONFUN may be incorrect or inaccurate. 
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3. The search direction has become inaccurate because of 
illconditioning in the Hessian. This tends to be reflected in large 
values of ItQP in the output. In this circumstance it may be 
worthwhile to rerun NPSOL from the final point. 

4. The matrix of constraints in the working set is ill-conditioned, as indicated 
by an extremely large value of Cond T. This matrix is used in the QP minor 
iterations to find the search direction. In thin circ umstan ce it may be 
worthwhile to run NPSOL with a relaxed value of Feasibility Tolerance. 

In addition to these recommended measures, a capability to rescale the targeting 
and optimization problem and resubmit to NPSOL has been incorporated into 
IPOST. This improves IPOST's ability to solve the problem. 

More user flexibility is provided by a file called NPINPUT. The directives in thia 
file (Table 3 - 12) allow direct access to some of the NPSOL options. These 
directives should only be altered by the experienced user. It is recommended that 
the NPSOL user's guide be a supplemental text. 
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Cold Start 


The user does not provide NPSOL with die initial working 
set. 


Derivative Level i 


Specifies whether IPOST provides NPSOL with analytic 
derivatives. Always input i = 0 except when using ITOO. 


The user can specify the perturbation interval used to 
estimate the gradients and sensitivities. If the user does not 
specify a value V, then NPSOL will estimate the finite 
difference interval. For a class of similar problems specifying 
Difference Interval will save execution time. 


Difference 
Interval r 


Feasibility” 
Tolerance r 


V defines the maximum acceptable absolute violations in 
linear and nonlinear constraints at a "feasible" point. Using 
this directive sets both the Linear Feasibility Tolerance and 
the Nonlinear Feasibility Tolerance. 


Hessian No 


'Hessian No’ will help NPSOL run faster 
'Hessian Yes' is recommended for later warm starts 


Hessian Yes 


Major Iteration 
Limit i 


The mmrimnm number of major iterations allowed before 
NPSOL termination. 


Controls the amount of mcgor iteration printout NPSOL ~ 

produces. 

= 0 None 

= 1 Final solution only 

= 5 One line each major iteration. No final solution 
printout. 

> 9 Final solution and one line for each iteration 

> 19 The objective function, eudidean norm of the nonlinear 
constraint violations, values of the scaled nonlinear 
constraints and values of the scaled control variables each 
iteration. 


Msgor Print Level 


Nonlinear 
Feasibility 
Tolerance r 


'r' defines the maximum acceptable absolute violations in the 
nonlinear constraints at a "feasible" point. 


Optimality 
Tolerance r 
Warm Start 


'r' is the number of correct figures desired in the objective 
function at the solution. 

The user provides NPSOL with the initial working set. 


Table 3 - 12. NPSOL Directives in NPINPUT File 
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4.0 SAMPLE CASES 

The following sample cases illustrates a mission application of IPOST. This is by no 
means intended to cover all capabilities nor is it realistic in every detail, but it does 
provide a meaningful example for constructing and unders tanding a typical mission 
application. The sample case represents a Voyager 2 mission. 

4.1 APPLICATIONS 

Before the Voyager case is discusses, the context within which IPOST is applied 
should be explained. 

Examples of mission applications are shown in Table 4-1. They illustrate some of 
the diverse IPOST model capabilities, including lunar, interplanetary, orbital, 
thrusting (impulsive, low, high), and gravity assist. 

The typical IPOST application is usually in the form of run sequences which 
compare various mission options. These analysis threads build on each other, 
culminating in a reference mission which is used to support detailed system design 
and analysis. Table 4 - 2 illustrates representative mission threads which exercise 
key IPOST capabilities. 

For example, in the Comet Rendezvous thread, each case is actually several runs of 
IPOST to generate parametric data such that mission decisions and refinements can 
be made. The sequence of multiple runs per case feeds each succeeding case with 
each case becoming more realistic in terms of model fidelity, and encompassing 
more system objectives and constraints. 

The first step in the Comet Rendezvous thread is to define mission requirements 
(necessary conditions) and goals (desired conditions), as well as any known 
constraints, such as technology status. In this thread, a mission requirement would 
be to successfully rendezvous with a specific comet in a specific time frame. A 
mission goal might be to collect data on planets or bodies that are enco un tered 
during the interplanetary trajectory from Earth to the comet. A mission constraint 
mq?ht be the availability of a Cesium ion thrust engine powered by solar arrays, 
which provide limits on power/thrust levels and on specific impulse. 


The first case is an approximate impulsive AV solution using a Venus gravity assist. 
This examines coarse energy requirements, benefits of gravity assist, and optimum 
mission opportunities (launch-arrival dates with payload/launch mass sensitivities). 
The mission may be analyzed as separate phases, e.g., Earth to Venus, Venus to 
comet approach, rendezvous and stationkeeping. 

The second case models a low thrust mission using a single thrust segment with 
variable steering and variable throttle. Implicit optimization (collocation) is used. 
This recognizes the coarseness of the initial guess, and provides rapid solution 
searches. A determination is made whether available technology is sufficient to 
provide the required payload mass at comet rendezvous. 
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As mission knowledge evolves, the third case introduces multiple coast/thrust 
segments. These added degrees of freedom provide more flexible, and more 
realistic, mission solutions. The optimization method can be implicit or explicit 
(with master-subproblems), depending upon how many, and what level, of mission 
decisions need to be made. This would include interplanetary and dose encounter 
geometries, flight times, subsystem performance, etc. 

The final case provides an end-to-end predsion optimized reference trajectory for 
system analyses and subsystem design support. Sdence and mission objectives can 
be assessed with a high degree of confidence. 
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TEST MISSION PLANETARY 

CASE DURATION BODIES ACTIVITIES SIMULATION OPTIMIZATION 
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• Lunar Mission Thread 

LI Earth departure to lunar orbit with patched conic 

L2 Space Station to moon with free return 

L3 Space Station with finite bum escape to libration point 

• Voyager II Thread 

VI EJS portion with ITOO (analytic partials) 

V2 EJSUN with finite differences partials 

* Comet Rendezvous with Solar Electric Propulsion Thread 

Cl Approximate impulsive AV (DSM) with Venus gravity assist. 
C2 Single thrust segment, simple collocation 
C3 Multiple thrust segments, complex collocation 
C4 Explicit optimization comparison 

* Human Mission to Mars Thread 

Ml EM launch/arrival date search, simple aerocapture 
M2 MVE return leg optimization. Earth orbit capture 
M3 Mar6 surface ascent to orbit (MAV design) 

M4 EMMVE round-trip optimization 



Bodies 


E-Moon 


SS-Moon-E 


Propagators) Forces 



E-V-C 


VI 


V2 


Cl 


C2 


C3 


C4 


Ml E-M 


M2 


M3 


M4 


A = Aerodynamics 
J2 = Zonal (J2) 
HT = High thrust 


1-ste 


Multiconic, 

Cowell 


1-step, Conic 


1-ste 



Optimization 

method 


M 


M 


MS(2) 


MS(3) + ITOO 


SB, J2 


DB.HT 


SB 


B 


SB 


DB, LT 


MS(2) 


collocation 


DB, LT, SP collocation 


DB, LT, SP M 


Search 


MS(2) 



1-step, Encke I SB, J2 


LT = Low thrust DB = Disturbing body 
SP = Solar pressure M = Master only 
SB = Secondary body MS = Master + Sub(s) 


Table 4-2. Example Mission Threads 


The Voyager 2 case illustrates a master subproblem formulation and multiple 
planetary encounters. Only the Earth-Jupiter-Satum phase of Voyager is 
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performed. The simulation has 8 events, starting in Earth park orbit and ending 
with a Saturn flyby. Total impulsive delta V is minimized in the master problem. 
The two subproblems target each of two legs, Earth to Jupiter and Jupiter to 
Saturn. 

In setting up an IPOST problem, the trajectory simulation is defined first, as 
opposed to the optimization process, because it describes the primary mission. 

The first $TRAJ namelist is event 5. This utilizes parameters such as S/C mass and 
propulsion characteristics. The initial date of July 31, 1977 precedes the actual 
Voyager launch date. The S/C is placed in a circular orbit about Earth, and 
trajectory propagation will use conic, or two-body, equations of motion. 

The second event (#10) is triggered by a flight time of 20 days. The 1STEP 
propagator is activated in connection with activation of the LAUNCH mode. Hence, 
Voyager orbits the Earth for 20 days and then is impulsively injected onto an escape 
hyperbola. Using flight time as a control parameter would allow variation of initial 
launch date. For 1STEP, the primary body is the Sun and Earth is the secondary 
body. 

The third event (#15) is initiated after a flight time of 20 days from the launch 
event. At this point, the S/C is well outside of the Earth's sphere of influence. The 
secondary body is now defined as Jupiter for 1STEP propagation. The reference 
body for trajectory calculations is also set to Jupiter. 

At event 20, the triggering criteria is mission time. The intended value is Jupiter 
closest approach time which is not specified explicitly, but indirectly through the 
optimization process. A conic propagator is used with Jupiter as the primary body. 

Oneday later, at event 23, an impulsive trajectory correction maneuver is executed. 
1STEP is reactivated as the propagator as the vehicle flies away from Jupiter. The 
combination of gravity assist and midcourse correction will set up the traiectorv for 
the flight towards Saturn. 

The fifth event (#25) occurs 300 days after the midcourse maneuver. The secondary 
body for 1STEP is switched to Saturn, as is the reference body. A general rule of 
thumb is that planetary sphere's of influence for 1STEP (and Multiconic) are about 
20 days for small inner planets (Mercury to Mars) and 200 days for large outer 
planets (Jupiter to Neptune). 

The next event (#30) is Saturn's closest approach, on August 25, 1981. It occurs 
1467 days after Earth park orbit escape. 


A final event (#90) is used as a "dummy" event, and coincides with the previous 
event. This event is needed because final event computations are done only on the 
"minus" side of the event. The namelist specification of "NONE" means no more 
input follows; that is, this is the end of IPOST problem specification. 

We now return to the optimization definition, or $TOP namelist. The NPSOL 
algorithm is used with finite difference perturbations bong calculated internally by 
NPSOL. This is the normal mode, as opposed to the user specifying perturbations. 

A maximum of 10 iterations is allowed for master problem optimization and 600 
iterations for each subproblem targeting process. A subproblem targeting iteration 
is typically much shorter than a master problem iteration because the latter must 
solve all subproblems successfully. 

The mas ter problem controls are Jupiter closest approach time (defined as the 
criterion of event 20), B dot T and B dot R of the Jupiter approach, or incoming, 
hyperbola (defined as the 2nd and 3rd dependent variables of subproblem 1). In 
addition to the initial guesses for control values, the upper and lower bounds, and 
weighting values are important inputs. These define the performance manifold and 
often mpan the difference between problem convergence and divergence. The 
objective f unc tion is the s um of all delta V magnitudes, which in this case 
corresponds to the launch/escape maneuver at Earth and the Jupiter midcourse 
correction. 

The two subproblems are defined next. The first subproblem ends at event 20 
(Jupiter closest approach) and the second subproblem ends at event 30 (Saturn 
closest approach). IPOST automatically assumes that subproblems are non- 
overlapping. Both subproblems use a Newton-Raphson targeting technique, as 
opposed to subproblem optimization with NPSOL. 

The control variables for subproblem 1 are the V-infinity vector of the Earth 
departure hyperbola. Controls for subproblem 2 are the AV components after 
Jupiter flyby. As in the master problem, important inputs are the control initial 
guesses (USUB), bounds (INDSLB and INDSUB), and weightings (WGTS). 

The constraint or target variables for subproblem 1 are time from periapsis, B dot T 
and B dot R at Jupiter. For subproblem 2 the constraints are time from periapsis, B 
dot T and closest approach distance at Saturn. The use of B dot T as a control 
parameter (with loose bounds) at Saturn affects oxbit inclination as well as which 
side of Saturn the S/C flies by. B dot T and B dot R are often used because of their 
stability in the targeting and optimization process. The choice of constraint 
parameters for subproblem 2 reflect termination of the mission at Saturn flyby. For 
the actual Voyager mission, which continues on to Uranus and Neptune, the 
constraint parameter types at Saturn would have been identical to those at Jupiter. 

Because IPRINT was not input, the default value of IPRINT = 0 is used. This will 
result in only summary information of the master problem, plus the final trajectory, 
being output. Except for well tested production runs, it is recommended that more 
detailed levels of print be exercised. 
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The first page of output summarizes IPOST input, including initial conditions, event 
summary, targeting and optimization definition, master-subproblem structure, find 
NPSOL options. 

The second page of output completes NPSOL parameter definition and then prints 
an iteration summary for each major (master level) iteration. This includes the 
objective function value (sum of delta V) and other conditions of the optimization 
process. Of some interest are the number of objective function evaluations (NFUN) 
and the condition number of the Jacobian (COND T). This problem does not have 
any nonlinear constraints, only control bounds, but if it did, the value of the 
nonlinear constraint norm is displayed, as well as the constraint values. 

As the iterations progress, certain key parameters should be monitored. These are 
the objective function value (which should decrease), the condition numbers of the 
Hessian and Jacobian (which should remain small), and the value of the nonlinear 
constraint norm (which should decrease). Also of interest are the convergence 
indicators at the right side of the summary. When all these flags are "T", then 
successful optimization has been achieved (according to NPSOL). One cautionary 
note is that reliance on a few indicators, such as the convergence flags, can be 
misleading. It is important to examine all measures, including the final solution. 

The iteration summary continues until page 4 of the output, where a iwavimiim 
iteration limit has been reached (INFORM = 4). NPSOL is exited, and the final 
trajectory is displayed. Conditions at each event are output, including the minus 
and plus sides of each impulsive maneuver. At the end of each subproblem, a final 
iteration su mmar y is displayed. In this case, Subproblem 1 convergence is 
summarized between trajectory blocks for event 15 and event 20 (Jupiter closest 
approach), and Subproblem 2 is summarized just before event 30 (Saturn closest 
approach). 


The last sdution has a total AV of 7.27 km/s, reduced from the initial solution of 
7.67 km/s. Jupiter closest approach time is 711.9 days from launch, and B dot T B 
dot R were selected as 1898454 km, 130291 km. 


Finally, a master problem summary output is given which includes the number of 
iterations and CPU time. A valid solution may or may not meet all convergence 
criteria. Only the user, with adequate engineering experience, can make that 
judgment. 


69 



-THU 


Utt'# TT* r/P^AU/ 

:p r 


70 




o 

c 


xJ 

u 

€> 

a 

o VO 

XJ I 

* 3 -o ”0 
*o «a * o o • 
1-4 o> * on«H 


i> i E <h a ^ 

i C 41 ~ O «fl 

: — 4 xi jc c EU 4 JH 

i u ® aBiJ o 

i au «) « « a x x x> 

► — I —4 T 4 VM — I c £ EM-1 




O * CM 

* 


X 

* * * k H • 

—4 


> 

o • • • o * 

XJ 


TJ 

T3 O O O O 4-4 

•o 



• HHHOO 

A 



O 1 o • 

* 


* 


* 


* 

* • • • 44 

* 


o 

i-4ooo com 

a 


N 

IT »— 4 i— 4 f— I r-4 

4*4 


M-i 

CM 1 - * 

4J 


c 

VO % * *i-4 VO 

* 


-H 

»-4 • • • O TJ 



> 

•00000} 

* 



tnHHHOM) 

* 



1 O CM 

—4 



***** O' 

u 

% 

o 

f-4 • • • CM 

TJ 


>4 

O' O O O *«H 

XI 

xz 

4*4 

m CM CM CM 1-4 • 

* 

a 

* C 

lO 1 o * 

* * 

« 

<N -4 

v or- 

CM * 

u 

* > 


*— 4 * 

c 

CM * 

n w • • • o O' 

CM O 


* «, 

cm o o o • 

*TJ co 

* 

CM * 

* * CM CM CM CO 

CM XI * 


o m i *o 

* ro 

X •> o » 

* X 

r- i-4 co 

* * 

ao o o 

4-4 VM 

* Cl * * * O i - 4 

4—4 * * 

« ro TJ TJ 

k COM) • • * O • 

% ao 

u *o o % 

H<HH(MOOOO * 

r-4 4-4 CM 

CO * • o 

* > * • CM CM CM O «H 

* x> * 

* (\ HHH 

*-4 * 

CO CM 1 • O 

4—4 * CO 


XkH 

-H U JZ -H 3 

x > a a ax> 

TJ TJ TJ tTJ Up 

c c c c c e > 

■H *H — I 3 —4 — 4 CL 9 


V 

N 


i—l 

« 

•H 


X) 

H 

£ 

Li 

o 

4i C O 

-4 

« JZ 

U 

1—1 « Ll 

c 

> a u c 

a 

4> l> Vm 3 


4J 4J 4J a 40 

X) 

C 

£ 

a a a o *j 

D 

o a o o a o 


o o o s *> 

« 

£ « 4J xi C U 

a a 

o 

u o 

C 

c 

c 

c 


■h ux: X5XJ - 

« > a p 

X » •> -D « » « 4 
TJTJTJ 3TJTJ p 

c c c n c c > 

•H — 4 —4 3 —4 —4 % 


PRECEDING PAGE BLANK NOT FILMED 



TJ T3 T5 O T3 
• • • O CM 

O O' H H • 

O O' o 
o as o 
o O' o 
o O' o 
cm «-h cm 

vo vo vo r vo 
*0*0*0 0*0 
o o o o o 


vo 

o 

eg 


VO VO VO «— I VO 
75*0*0 • *0 
00 OD CD 00 

VO VO VO VO 

cm cm cm cm 

OV O' O' O' 

CM CM CM CM 


*0*0*0 ■ *o 

O' O' O' O' 

^ 

00 00 00 00 

o o o o 

00 00 00 00 


* * *.«n 

H fH H O 

o o o o 
o o o o 
o o o o 


I I I 


*o 

4J O 

V * ->f-l 
0»«-g m « 

* H 

«. m 

*o * 

C r- 
0) * 

CM 

C 

« 

i 

o 

* 

4J % 

a u 


Wi 1 CM «. 

r» 

-h c 

vo in • 

m * *o 

o r* 

u o 

*0*0 0 

4JHOH 

% O O' 

« * * 0 

• • 00 

o 1 CM • 

in ^ *h 

« m m * 

rH CM 

III! 

INI 

III! 

( 1 1 


> 


» © 

o * 

«H • 

c <n 

0 vo 
u in 
* vo 

1 I 


**-« a 

u « 

E -» 4J 


X 

u 

c 

3 

«e 


1 *** ** * 


vo vo « 

ID ID H 

vo vo i 


x x ^ o 


i 4 -5 

X 

E >.>i 

r *j 

X 

X 


>i 

a 


a 

>H 30H 

V 

> a «a 

V o 

40 * 0*0 a 

« ■» 

O 

o 

lg c 

a*o 

s£x 


40 

m m ot o» c 

« 

« W 0) 

■n C O O 

u O 0 0 

4 3 

3 

■n tg 

C 4J 3 

o o 

4J 0 

tg 

o 

o* a a a • 

T-l 

*j 4j a 

* «j au 

*HjQX) M 

E M-h 

a 

<0 40 

0 *H H 

M X 

c c 

U 

a 

««»««> 

x 

a a a o 

M > ft) « 

*c a*o a 

o x a 

c 

M « 

> *g 40 

a a 

<0 

a 

40 

•o *o *o *o * 

o 

o o o * 

4J « •*« *0 

-»H -H 

« 4J 0 

-H X 

O 

* O > 


E -H 

u 

u • 


a u o 


ID- V> 

a o 


72 


a 

u 

3 

**0 * tn 


ft) U > -H 


I 


change idbody 


jupiter tea 
7 / 9 / 1979 . 


«— I tO CD 

•d w 

EfOH 

o 

*> 

*o 


E 4J u 
•B C 0) 
^ 4) C 
•0 > CD 
a O O' 








*■ 






% % % 




* 


* 


r-t •. 








* 

* 

M-l U 


a ^ a 


a 


U4 * 


a 

* 

U 


M 0 V 


M *. 


M ■ 


M 

© 

EC*. 

>t 

p au «* 


3 • * 


Er- 

*o 

3 

c 

%-H o O 

XI 

% T5 »E Din 


«>T3 O SO 


0D 

c 

**o *> 

o 

o ii u * 

o> 

m 4J • *-H *. 

>* 

m iio ». «. 


O iiT 

© 

O ii ■ 

c 

cn * *. m 

c 

(N* H*. * O 

TJ 

cm *■ m o vo 


«n *. *h 


<Ti * O 

* 




O 




© 



• i i l 

* 

1 1 1 1 I 1 

X) 


<d 

l l l 

x: 

ii a a 

a 


•» 


*© 


O »H 


ij 








4J 00 






T3 




ON 


n 




• 


« 


C *H 






M 

P* >1 

O' 

>i >. 

M ^ 




4-1 

jj M Q|*0 

© 

4J u © S O.'O 

c 

4J M © *0 *T3 

3 m 

jJ M © 

m 

ii M v 

« 

C -U O O 

T*» K 

C ii 7 U O O 

-r-> id 

CiiPOC 

f-* 4J eg 

Cii 3 


Cii 3 

H 

O^H Mi) 

id o 

c MX) 

«d x: 

© — 1 <-H Xl .Q 

id * ^ 

© *H r-H 

© X 

d *H H 

E 

> m a* a 

m a 

> m « « a a 

M o 

> m «d a*o 

M MOD 

> M d 

M iJ 

> m m 

<d 

® o **H 

4J 

• o > £ «h -h 

4J 

® o > 

jJ 

4) O > 

4J 

© u > 

c 


vy w 


v> v> 


vy w 


w v> 




a u 

u 

a o 

u 

a o o 

O 

a o 

o 



73 


requirements for problem 



0 

O' 

<0 

a 


m 

cn 


u 

o 

3 

4 = 

H 


U U 
0 *H 
0 44 

CL 
O *M 
O rH 

o 0 
o « 


E x: 

44 


Li 

JC 


0 
» o 
»h c 

m « 

» Li 

0 


m 

m c 


o o o o 

+ h- -► + 

O O Q Q 
O O O O 
O O O O 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 


+ I 

a a 

CD m 
O CN 
CN m 

m *h 

rH m 
O <N 
m cn 
cn to 


>000 to rH 


>000 
>000 
' + + + 
>000 

> CN O O 

> CN O O 
(WOO 

• woo 
10000 

• moo 

• moo 

• CN O O 

• moo 

• o\ o o 

r- o o 

r*~ o o 


o o 
♦ -l 
0 a 

in m 
cn vd 


M 

O 

*uaaoaaaa 

O'-H s (|y 4 9) II M 
<0 C iJ U C 44 44 44 4 J 

Ouonnomnnn 

O Oh h O «h rH rH «H 


43 0 

U rH 

V c I) I) P « Ml 

. coccaccc 
>*O 0 OOEOOO 
C «H C C *h c c c 




m m 

*h r- 

VP CN O 
W rH *H 

cn o li 

H«iJ 

U I 


m o m o Our- rH r 

0000-H0044 

-L + + +rHH-+-r 4 

QQQQOQOC 
OOOO VO H 3 

OOOO WflD 
OOOO «0 VP 4 J 
OOOO DOW 3 
OOOOHiTOD CL 
OOOO Ll VP VP .u 
OOOOiJHM 9 

0000 crn o 
m 0 m r 


«H o 
+ -► 
o o 
o o 
o o 
O o 
O o 
000 
3 0 0 
— < o o 
<000 
>00 
o o 
o o 


HHO 

000 
+ + + 
000 
000 
000 

OOO 

000 

OOO 

OOO 

OOO 

OOO 

OOO 

OOO 

OOO 


cn m o 

OOO 
+ + + 
000 

OOO 

OOO 

OOO 

OOO 

OOO 

OOO 

OOO 

OOO 

0 r- o 
O 00 o 
OTO 


■H 

E 

0 «r • £ 

W O VP O 

0 0* OV «. 


t 

4J CN • 0 

in 0 m 0 

0 m cn 0 


3 

3 • • Li 

♦ 

• • • 

■H • • ^4 

■O 

0 

a • «m 

\p O VP O H OV CN u 

C 


C • 0 



0 44 

0 

44 

-H 44 *0 



JZ 0 


3 

0 0C 



I 

0 

a 

>* 4J *0 0 




4i 

c 

0 0 

c 


0 1 >« 


•H 

O 0*0 Li >,44 

0 


44 Ll 

*0 


44 O 0 TJ 0 -H 


0 0 0 


4 4 

0 a c *0 0 44 

n 


44 44 g 

c 

0 

000 C43 0 

0 

O 

0 -HE 

0 

O 

-m -H 0 

4 4 


C 3 

*H 

a 

0HHHHH 

Ll 

c 

rH 3 0 

44 

>H 

LI 0 3 0 0 0 

0 

0 

0 

3 


44 -H *n 0 rH -H 

O 

0 

-H 44 44 

U 

* 

4J 44 4J 



44 PC. 

0 


* -H -H -H 



-h a 0 

X 

* 

* c c c 



c c > 

0 

* 

* -H rH «H 



•H -H 0 


H(S|(N<MH|C|HO 


u a«M a a«M a 

0'EE3E33E3 
t 7 *-H H *0 «rTO TJ «H xi 

u 

44 


OOOOOOOO 

OOOOOOOO 

OOOOOOOO 

momommoo 

rH *-» CN CN CN m cn 


1 


44 

c 

0 

> 

0 


VP o 

o o 
-f + 

w w 
r* 0 
CN o 
O o> o 
CN CN O 


2 


m o 
cn «h m 
u 

2 

Li 


E *-h 









0 0 
+ + 





3 O 









M U 





0 0 









in 0 





> a 









00 0 





*0 C 






O O 


OOO 



m 








>lH 


cn a> 0 



CN 


1 






Lh 

•H 


rH rH 





Ll £ 






C 

44 


• ■ 


Ph 


44 

0 



0 


•H 

V 


O O 


> 


*o 

0 -C 

Ll 



Li 


> 

jQ 




•O 


43 

♦ 44 44 

0 



0 










♦ 0-H 

44 



44 


CN 



m ^ 





* E Ll 

0 



0 





0 0 





0 0O 
44 Ll O' 

E 

0 



E 

0 





l 1 
U U 





p 0rH 

Li 



Ll 


0 



O O 





a a 0 

0 



0 


H 



O O 





c 

a 



a 


O O 


OOO 



m 


-H c C 



cn m 



0 X rH 


CN O O 



CN 


0 0 

a 


0 0 

a 


E 'Ll 



«h m 




C -H -H 

0 


rH H 

0 


0 c 

a 




X 


a 

O 44 44 

*0 


Ll > > 

•0 


»H «H 

'Ll 


0 0 


> 


'Ll 

•H 0 0 

c 


44 0 0 



43 > 

44 


1 


*0 


44 

44 N M 



*h aa 

44 


O 








0 -H *H 


0 

Lt 0 0 

0 

0 

M *H 




CN 




NEE 

H 

t 

0 *0*0 

O' 

£ 

a 







-H *H -H 

0 

0 


14 

*} £} 0 44 


0 

0 


0 

44 


E 44 44 

Ll 

C 


0 

c 

3 -H C 

0 44 

0 


rH 

C 

0 

*H a a 

44 



44 


0 Ll O 0 44 

c 

0 c 

Li 

O 

0 

44 

44 0 0 

C 





0 Ll > 

0 

0 

0 0 

0 

Ll 

> 

0 

a 

O 





L 44 0 

O' > 

3 Li 

A 

44 

0 

O' 

0 E E 

O 



E 


O E C 

Ll 

V H 4} 

1= 

c 


Ll 

>*0 0 




0 


3 O 4J 

0 


0 rH 

3 

0 

44 

0 

t7*rH H 

■ 





L» C O 0 44 44 

> O 

C 

0 

0 

4 4 

C JO 43 

JQ 



4a 


0 


0 

44 





•H O O 

O 



0 


Zee 

E 


E 

E 

E 


E 

44 Ll M 

u 


OOO 

Ll 


E 0 0 

0 


0 JZ 

0 

0 


0 

0 a a 

a 


OOO 

a 


3 'H i-H 

,-4 


rH 44 rH rH 


rH 

O' 



OOO 



C 43 43 

43 


43 -H 43 43 


43 

U Li Li 

N E 

• ■ a 

Ll E 

O O 

O 


O » 

O 

O 


O 

0 0 0 

0 

3 0 0 0 

0 

DH Ll Ll 

Ll 


Li 

Ll 

u 


Li 

44 44 44 

44 

C CN CN CN 

44 

C 

0 a a 

a 


a 

a a 


a 

♦ 0 0 

0 4 4 


0 

44 

44 X>43 

43 


43 

43 43 


43 

* 5 S 

0 

> 


0 

> 

0 3 3 

3 


3 

3 

3 


3 

* EE 

E 

0 


E 

0 44 0 0 

0 


0 

0 

0 


0 


74 


subproblem values -0.10000E-03 0.20613E+06 0.20000E+06 

with tolerance 0.50000E-04 0.10000E+03 0.10000E+03 



VERIFY LEVEL 0 

DERIV LEVEL 0 

DIFFERENCE INTERVAL l.E-4 


CN 


o 

H 

O 

H 


m 

03 

w 



o 

O 



o 

o\ 







H 

ev 



• • 




• • 




• • 




• • 

c 



* • 

o 



• • 

*mI 



• t> 

V) 



• o 

-H 



• c 

o 


o 

• *o 

0) 

H 

C 

• U 

kl 

01 

41 

4J 0> 

cu 

> 

3 

U *-H 


V 

tr 

<0 O 

c 

H 

o> 

4-1 4-1 

o 


wi 

<n 

^-4 

>, 

*M 

X 

4-1 

VM 


Q « 

0 

-H 

t) 

^ <0 

c 

U 

> 

O H 

3 

01 

*0 

o o 

bu 

> 

CO 


> 

c 

o 

CJ 

H 

T> 

C 

o 

u 


N 

X 

T> 

c 

o 

u 


r~ vx> 

( i 

w u 

o o • • O 

n HDH j: 


03 

u 

§ 

o: 

w 

m3 

o 

fx 


H 

M 

x 

M 

m3 


>* 

€-* 03 

*“» u 

m 3 z 
* »H < 
3 ffl cC 


o 

w 

CO 


Jl.. 

V) UfflU 
Z > M 03 
O UWh 
H J<0 

h w h 

< H Cb 

os z >« e-* 

w h At h a 

HOCSh < 
HfoWJZH 

z < in 

t£ oc m £ m 

O O J h W D 

*r> z H w m3 

< < O Ou W O 

X X Z O X U 


I 

4-1 

£ 


0) 

4J 




noo 

VO P- rH o 

o o o 

to 



O CN CN 

o o o 

f-H 

o 



1 •+ + 

i I t 


1 



uuu 

WWW 


w 



VO o o 

o o o 


H 



<7> O O 

o o o 











in h h 

r- 1 cn ov 


lO 



i * • 

, , , , 

. . « 

4—4 



I . a 

toil 

4 • H 

4) 



• a • 


• • -rl 

> 



a a a 

O 4-1 4 • 

• • *4-4 

kl 



. o . 

U -H O • 


0) 

cr» 

i 

>iN 0 

c*^ c * 

• • c 

4-1 

m 

R 

iJ-H N 

<0 *<H <0 • 

H H 0 

c 

O' 

■ 

*H « -H 

k4.fi UH 

o> oi 

-H 

«-H 

1 

«— i n 

01 -1 0) 01 

> > *0 



a 

-H *T3 

mh n mh > 

0) 4> 01 



a 

43 c a 

0 4 0 91 

H H ^ 

o 

> 

a 

•H 3 0) 

4J 0) jj mH 

<0 

*w 

o 

a 

11 Oil 

VM 

4-1 4-1 « 

*4-1 

z 

a 

*0 43 w 

>4 X 01 

C C 

•H 


a 

V 

4J U 0 > 

•h -h a 

•o 


a 

o a 

«H <0 kl -H 

k4 U43 



a 

4J 4J 

h a 4 -u 

o. a 

«— 1 

m 

a 

M -H 

« c a o 

o 

* 

o 

a 

10 C C 

£•»<«> 

kl kl 4J 

u 


a 

01 ->-4 -X 

-HH a -H 

o o 

4-1 

*r 

a 

C *w *m 

4J C C M 

•nCZ 

c 


a 

•H C C 

ao-H a 

<0 *H Z3 

0) 

c 

a 

JHH 

OZ JO 

X X DC 

u 


TJ 


0 « 

-H • 

n « 

ki i 

v i 

> > 
a 
i 

1 a 

i i 

i i 
« 
i 

-3 I 

O * 

10 I 

04 fl 

z « 


Z Z 


O ON 

o 

ID 

CN 


H 

a. 

o 



ono 

ocorom 

o o o 




o 

H 

m m 

o 



+ 

1 


1 



03 

03 


03 



O 

m 


o 



o 

kO 


o 

z z 



• 


• 

M M 


CN 

CO 


i—4 







T3 



• dum. 

• • • 

• 

• 



• ki b c 

4-> 4-1 * 

• 

V 



•><0*0 

--4 —4 • 

• 

*0 c 



4J > > -H 

E £ 0) 

• 




c « 

•H —-4 H 

i—4 

0) 


n • • 

•H C > *H 

HH-H 

<0 

T3 * 


4J • • 

*0 *0 -H o 

*4-1 

> 

0) _ 


c * • 

kl -H 4J 4k 

m « 

ki 

*0 E 


-r4 . • 

4-i 43 O ki 

C c E 

0) 

-4 0) 


<0 • • 

« o 0) CL 

o o o 

4-1 

> 


u • • 

c o **"* 

•H-H U 

c 

O 43 


4-1 * * 

O *0 43 4> 

4J 4-> VW 

-H 

ki o 


n • • 

O O C 

*0 «0 


a k4 

« 1 

C • 4-> 

-*-4 

k| k4 *0 

0) 

a 

U 1 

o •» -4 

ki U U JC 

01 4) 4) 

u 

0) 

a i 

o t> E 

e «o * o 

4-1 4-1 T3 

c 

O 0) 

jj i 

H ■H 

01 0) 01 <0 

•H*H a 

0> 

<0 > 

4! * 

ki 43 

C C C £ 

o 

kl 

0u— 1 

1 1 

<0 <0 

«H «H "H ' 

ki kiH 

1) 

*) O 

<0 1 

«-h a 

H H >— 4 

0 o 

*4-4 

4C « 

U 1 

C M « 

C C C CO 

T-» C Z 

*4-1 

ki 

<0 1 

<H Oi> 

0 O o 04 

<0 ^4 3 

— — 1 

O O 

04 1 

m3 > 10 

Z Z Z 03 

X X DC 

Q 

Z fx 


75 


The user sets 0 out of 3 objective gradient elements. 

Each iteration, 3 gradient elements will be estimated numerically. 
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initial conditions for phase 
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state relative to idbody: jupiter 

x - 0 . 433346698E+06 y - 0.555192933E+06 i — 0.910720487E+05 radius - 0.710156794E+06 scmass - 0 . 213473875E+06 
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esn * 90.000 fesn- 90.000 

time- 2.44484250D+06 
normal termination 
cpu - 40.300 second* 
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5.1 NAMELIST STOP 

The following is a description of input parameters for Namelist $TOP . 


VARIABLE STORED 

NAME TYPE VALUE DESCRIPTION 

ATARWT REAL Contingency targeting adjustments 

parameters. Applicable only for 
MODELT(i) = TARG1S’. 

2. D0 ATARWT(1) = fraction of target error 

that causes DR adjustment 

7. D0 ATARWT(2) = maximum number of target 

error increases (INCR) allowed 

3. D0 ATARWT(3) = DR scale factor for 

target error increase 

1.D0 ATARWT(4) = fraction of DVMAG that 

causes DR adjustment 

1.D1 ATARWT(5) = fraction of DR that 

causes DR adjustment 

8. D0 ATARWT(6) = maximum number of DR 

increases (INCPR) allowed 

5.D0 ATARWK7) = DR scale factor for DR 

increase (big target error) 

0.7D0 ATARWT(8) = DR scale factor for DR 

increase (small target error) 

2*0.D0 T.B.D. 

DEPPH(i) INT. 25*900 The event at which dependent variable (target) i 

is to be satisfied. DEPPH is of type esn=INT. 

DEPSLB REAL 500*-10D10 The lower bound of variable DEPSVR. 

DEPSPH(i) REAL 45*999.D0 The event at which subproblem or dependent 

or variable (target) i is to be INT. satisfied. 
DEPSPH is of type esn=INT. 

DEPSTL(i) REAL 45*1.D0 The desired accuracy level within which 

DEPSVR(i) is considered to be satisfied. 

DEPSUP REAL 500*10D10 The upper bound of variable DEPSVR. 

DEPSVL(i) REAL 45*0 .DO The desired value of the subproblem 

dependent variable DEPSVR(i). 

DEPSVR(i) CHAR. 45* ' ' The Character name of subproblem 

dependent variable (target) i. 

(See Table 3 - 9). 
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VARIABLE STORED 

NAME TYPE VALUE 

DEPTL(i) REAL 25*1.D0 

DEPVLBG) REAL 25* 

-10.D10 

DEPVR(i) CHAR. 25* * ' 

DEFVUBG) REAL 25* 

10.D10 

FESN INT. 999 
IDEB INT. 0 


IEPHEM INT. 0 
IFDEG INT. 25*0 

INDPH(i) INT. 25*0 
INDPLBG) REAL 25*-10.D10 
INDPUB(i) REAL 25*10.D10 


DESCRIPTION 

The desired accuracy level within 
which DEPVR(i) is considered to be 
satisfied. 

The lower bound on the dependent 
variable DEPVRCi). 

The Character name of dependent 
variable i.(See Table 3 - 9). 

The upper bound on the dependent 
variable DEPVR(i). 

The final event sequence number of 
the master problem. FESN is of type 
esn=INT. 

Trial step summary printout flag (i.e, 
print 'normal termination' at the end of 
each trajectory, call NOMOUT at the end 
of CNFUNC and print gradients) 

= 0, do not print trial summaries. 

= 1, print trial step summaries. 

Flag to specify type of planetary 
ephemeri8 to be used. 

= 0, analytic. 

= 1, precision. 

A flag to allow angles with a range 
of 0 to 360 degrees to be used as 
target variables. 

= 0, no adjustment is made 
= 1, if ABS(E(i)) is greater than 
180.0, set EG) = EG) plus or 
minus 360.0. 

The phase at which INDVRG) is to be 
initiated. INDPH is of type esn=INT. 

The lower bound on the independent 
variable INDVRG). 

The upper bound on the independent 
variable INDVRG). 
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VARIABLE 

NAME 

TYPE 

STORED 

VALUE 

DESCRIPTION 

INDSLB(i) 

REAL 

500* -1.D38 

The lower bound on the subproblem 
independent variable INDSVR(i) 

INDSPH(i) 

REAL 

or 

INT. 

45*999.D0 

The phase at which INDSVR(i) is to 
be applied. INDSPH is of type 
esn=INT. 

INDSUB(i) 

REAL 

500*1.D38 

The upper bound on the subproblem 
independent variable INDSVR(i). 

INDSVR(i) 

CHAR. 

45*” 

The Character name of the ith sub- 
problem independent variable used 
for the targeting/optimization 
process.CSee Table 3 - 8). 

INDVR(i) 

CHAR. 

25*” 

The Character name of the ith 
independent variable used for the 
targeting/optimization process. 
(See Table 3 - 8). 

INDXD(j) 

REAL 

or 

INT. 

1,2,3, ...25 

An array containing the indices of 
the dependent variables to be used 
for targeting/optimization. The 
indices must be input in INDXD(j) in 
ascending order; e.g. INDXD = 1, 5, 

INDXRj) 

INT. 

1,2,3,...25 

An array containing the indices of 
the independent variables to be used 
for targeting/optimization. The 
indices must be input in INDXI(j) 
ascending order; e.g. INDXI = 3, 5, 

9, is acceptable. 

INDXSD(i) 

INT. 

45*0 

Mapping of subproblem dependent 
variables to individual subproblems. 

INDXSI(i) 

INT. 

45*0 

Mapping of subproblem independent 
variables to individual subproblems. 

INTT 

INT. 

0 

Flag to return to Master problem 


initialization if warning errors 
occur in the input data: 


= 0, do not return to master 
initialization. 

= 1, return to master 
initialization. 
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VARIABLE 


STORED 


NAME 

TYPE 

VALUE 

DESCRIPTION 

INOUT 

INT. 

1 

Flag for displaying processed inputs. 
= 0, do not display processed input. 

= 1, display processed input. 

IPRINT 

INT. 

0 

Global printout flag (each value 
will produce print from all smaller 


values): 

= 0, Print diagnostic and master 
problem summary output: 

(IDEB=0, IPROs-1, LISTCN=-1, 
MSGLVL=1, PRNTPD=.FALSE.). 

= 1, Print nominal summary output: 
(IDEB=0, IPRO=0, LISTIN=-1, 
MSGLVLsslO, PRNTPD=. FALSE.). 

= 2, Print more summary information: 
(roEB=l, IPRO=l, LISTIN=-1, 
MSGLVL=15, PRNTPD=. FALSE.). 

= 3, Print subproblem nominal 
summary output: 

(IDEB=1, IPRO=2, LISTIN=-1, 
MSGLVL=15, PRNTPD=. FALSE.). 

— 4, Print massive debug output: 
(IDEB=1, IPRO=5, LISTIN=1, 
MSGLVL=99, PRNTPD=.TRUE.). 

= 5, Same as 4 with more debug output 
-NOTE- input of the individual 
flags, IDEB, IPRO, LISTIN, MSGLVL 
or PRNTPD will override that facet 
of the IPRINT flag . 

INT. 0 A flag to control the printing of 

trajectories during targeting/optimization. 
=-l. Print the final trajectory. 

= 0, Print the first and final trajectory. 

= 1, Print all nominal trajectories. 

= 2, Print all trajectories. 
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VARIABLE STORED 

NAME TYPE VALUE DESCRIPTION 

IRSCL INT. 2 Control weighting (WVU) and constraint 

weighting (WVNLC) rescaling option. Used 
with MXTTOP (I), I > 1. Rescaling based upon 
elements of Jacobian matrix and Objective 
gradient. 

= 0, rescale WVNLC based upon largest 
constraint elements. 

= 1, rescale WVNLC and WVU based 

upon largest constraint and control elements. 
= 2, rescale WVNLC based upon 

renormalization of Jacobian elements. 


ISMDEP INT. 0 

ISMIND ( ) INT. 25*0 

ISTATE INT. 60*0 


ISTM CHAR. 'ANAL' 


ISUB(i) INT. 15*0 


An index identifying the master 
level dependent variable DEPVR(i) to be 
used in the parameter search. 

Index of master level independent 
parameters that are reported/recorded after the 
search (ISMDEP * 0). 

NPSOL argument "ISTATE", indicates the 
status of every constraint with respect to the 
current prediction of the active set. = -2, This 
constraint violates its lower bound by more than 
FEATOL(j) in a QP subproblem. 

= -1, This constraint violates its upper bound by 
more than FEATOL(j) in a QP subproblem. 

= 0, The constraint is not in the predicted active 
set. 

= 1, This inequality constraint is included in 
the predicted active at its lower bound. 

= 2, This inequality constraint is included in the 
predicted active at its upper bound. 

= 3, The constraint is included in the predicted 
active as an equality. 

Flag to define method of generating state 
transition matrices. 

= 'ANAL', analytic matrices. 

= 'CENTRAL', central finite differencing. 

= 'FORWARD', forward finite differencing. 

Flag to control subproblem solution printout. 

= 0, do not print each subproblem solution 
* 0, print each subproblem solution 
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VARIABLE STORED 

NAME TYPE VALUE DESCRIPTION 

ITARGKi) INT. 15*9 Flag identifying the type of target parameters 

for subproblem i in ITOO targeting. See also 
SRCHM = TOO* and MODELT. Overrides 
DEPSVR. 

= 0, inactive. 

= 1, X, Y, Z 
= 2, VX, VY, VZ 
= 3, C3, RA, DEC 
= 4, VINFX, VINFY, VINFZ 
= 5, RPERI, BTHETA, TCA 
= 6, SMA, RPERI, INC 
= 7, DVX, DVY, DYZ 
= 8, RPERI, BTHETA, FPA 
= 9, RPERI, INC, MEAAN 

ITOOWT CHAR ‘SVMASS’ Type of weighting for ITOO (SRCHM = 

‘TTOOO 

= ‘SUMDV2’, sum of squares of AV. 

= ’SUMDV’, sum of AV magnitudes. 

= ’SVMASS’, sum of AV mass changes. 

USTTN INT. -1 INPUTX input display selector: 

=-l, no input display. 

= 0, display tables. 

= 1, display tables and namelist input. 

= 2, display namelist input. 

LNCNLB REAL 10*-10.D10 Linear Constraint lower bounds, BL. 

LNCNTJB REAL 10*10.D10 Linear Constraint upper bounds, BU. 

MODELT(i) CHAR. 15*’NRAPH' Subproblem targeting model selection 

flag, subproblem i. 

* 'NRAPH', N ewton-Raphson method 
= TARG1S', Special Onestep method 
= SUBOP+’, NPSOL41 Op timiza tion 

MODEW INT. 1 Flag to control type of weighting to be used for 

independent variables. 

= 0, user input wei ghti ng. 

= 1, auto control wei ghtin g. 

MSGLVL INT. 10 This variable was superceded by 'Major 

Print Level’ in npinput file. 
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VARIABLE 

NAME 

MXITARCi) 

MXITOP(i) 

NAMLST 

NUNES 

NPAD(l) 

NPAD(2) 

NPAD(3) 

NPI 


STORED 
TYPE VALUE 

INT. 10 
ARRAY 10, 5*-l 


CHAR. TRAJ' 


INT. 

60 

REAL 

9.D0 


REAL 4 .DO 


DESCRIPTION 

Maximum number of subproblem targeting 
iterations. MXITAR is a singleton which applies 
to each subproblem. 

The maximum number of major iterations i =1,6 
to be performed on NPSOL pass (i), before 
rescaling. MXTFOP overrides the value of the 
"Mqjor Iteration Limit" in the file NPINPUT. 

See also IRSCL for rescaling and section 3.4.2. 

The next expected namelist. 

= NONE', no more namelist input. Solve 
problem. 

= TABLE', read namelist TAB next. 

= TOP', read namelist TOP next. 

= TRAJ', read namelist TRAJ next. 

NOTE: TRAJ and TABLE namelists can 
be read only in the first problem. 

The maximum number of lines per page 
of output. 

The desired average of minimum and maximum 
number of digits different between perturbed 
dependent or optimization variable values to be 
achieved by adjusting PERT(i). Increasing the 
value ofNPAIXl) will cause PERT(i) to become 
larger. This option is disabled if NPAD(l)is 
input as zero. 

The number of significant figures different 
between nominal and perturbed dependent or 
optimization variable values below which the 
variable is ignored in selecting minimum and 
maximum required for PERT adjustment. 


REAL 14.44 94D0 The number of significant figures different 

between nominal and perturbed dependent or 
optimization variable values, above which the 
variable is ignored in selecting minimum and 
maximum required for PERT adjustment. 


INT. 3 Number of previous iterations to compare error 

in N e wton-Raphson subproblem targeting 
method before dedaring solution divergent 
(MODELT = NRAPH’). 
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VARIABLE STORED 


NAME 

TYPE 

VALUE 

DESCRIPTION 

NPINPT 

CHAR*30 'NINPUT' 

The name of the NPSOL NPFILE input. 

NPVAR 

INT. 

0 

Number of output parameters for PROFIL 
(PRNC >0). 

NSMINC 

INT. 

0 

Number of search increments (ISMDEP * 0). 

NSMIND 

INT. 

0 

Number of ismind listed controls. 

OPT 

INT. 

0 

A flag to select the type of optimization to 
be performed: 

= -1, Minimize the variable OPTVAR. 

= 0, No optimization. 

= 1, Maximize the variable OPTVAR. 

OPTPH 

REAL 

999 .DO 

The phase at which the variable OPTVAR 
is to be optimized. OPTPH is of type esn=INT. 

OPTS 

INT. 

25*0 

A flag to select the type of subproblem 
optimization to be performed. 

OPTSPH 

REAL 

25*900.D0 

The phase at which the variable OPTSVR 
is to be optimized. 

OPTSVR 

CHAR 

25*0 

The character name of the subproblem 
variable to be optimized. 

OPTVAR 

CHAR 

» t 

The Character name of the variable to be 
optimized, (see Table 3 - 9). 

OUTNAM 

CHAR 

100*’ • 

Output parameter names for PROFIL (PRNC > 
0). 

PDLMAX 

REAL 

2.D0 

Trigger level for a second (central differences) 
perturbed trajectory. The magnitude of the 
difference between NPAD(l) and the average of 
the minimum and mnrimiim number of digits 
different. 

PERT(i) 

REAL 

25*.0001D0 

The perturbation (increment) to be added to the 
independent variable, INDVR(i), whose value is 
currently U(i). Used to determine the 


sensitivity dE(j)/dU(i). 
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VARIABLE STORED 

NAME TYPE VALUE 


DESCRIPTION 


PERTSB(i) REAL 45*.0001D0 The perturbation (increment) to be added to the 

subproblem independent variable, INDSVR(i), 
whose value is currently USUB(i). Used to 
determine the sensitivity dESUB(j)/dUSUB(i). 


PRNTPD 

LOG. 

.FALSE. 

Rftj) 

REAL 

0 . 

RFSOI(i) 

REAL 

15*0.5D0 

SMINCR 

REAL 

0 . 

SPFESNG) 

REAL 

15*999.D0 

SRCHM 

CHAR 

‘NONE’ 


A flag to activate printing of the number of 
significant digits difference due to the 
perturbations in the control parameters when 
generating the sensitivity matrix. 

The Upper-Triangular Cholesky factor of the 
Hessian of the NPSOL problem. This must be 
input if a WARMSTART is used. 

The fraction of sphere of influence beyond which 
XYZ targeting occurs, subproblem i, TARG1S 
subproblem targeting model. (See MODELT). 

Size of each search increment (ISMDEP * 0). 

The final event sequence number for 
subproblem i. SPFESN is of type esn=INT. 

Master problem targeting optimization 
algorithm selector: 

= 'COLLOC', Use collocation algorithm 
= TTOO\ special Interplanetary Option: 
analytical-partials, master and 3x3 
subproblems 

= 'NONE', No master level targeting or 
optimization. 

= 'NPSOL', Use NPSOL at master level. 


SMINCR 

REAL 

1.D0 

TARGMT(i) REAL 

1.D0 

TOLF 

REAL 

1.D0 

TOLU 

REAL 

1.D0 


Size of search increments. 

Target error magnitude threshold for c hang ing 
target weights, subproblem i, TARG1S 
subproblem targeting model. 

Tolerance for convergence of the subproblems, 
tested against the sum of weighted target errors. 

Tolerance for convergence of the subproblems, 
tested against the change in the subproblem 
controls. 
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VARIABLE STORED 

NAME TYPE VALUE 


DESCRIPTION 


U(i) REAL 25*0.D0 The initial values to be used for the variables 

specified by INDVRfi) at the beginning of phase 
INDPH(i). 

Collocation: The user inputs m initial values for 
the m user independent variables, U(l) to U(m). 
Next the user inputs initial values for the event 
nodes (and internal nodes if needed). IPOST 
can interpolate the remaining controls for the 
internal nodes(NSGPHO in $TRAJ). 

USUB(i) REAL 45*0.D0 The initial values to be used for the subproblem 

variables specified by INDSVR(i) at the 
beginning of phase INDSPH(i). 

WGTS REAL 500*1. DO The dependent variable scale factors used when 

optimizing to magnitude of dependent variable 
errors. 

WOPT REAL 1.D0 The weighting for the optimization variable. 

WOPT should be input as approximately one 
over the nominal value of OFTVAR. 

WOPTS REAL 25*1.D0 The weightings for the subproblem optimization 

variable. 

WVLC R EAL 10*1.D0 The linear constraint weighting. 
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VARIABLE STORED 

NAME TYPE VALUE 

WVNLC REAL 25*1.D0 


WVSNLCG) REAL 25*1.D0 


DESCRIPTION 

The dependent variable weighting for DEPVR(i). 
Collocation: IPOST will use a set of default 
values of WVNLC for the defects and event 
discrepancies if the user does not input values. 
The defaults are: 

WVNLC(nodei-fO) * magnitude of 
U(nodei+3) velocity vector. 

WVNLC(nodei+l) * magnitude of 
U(nodei+3) velocity vector. 

WVNLC(nodei+2) = magnitude of 
U(nodei+3) velocity vector. 

WVNLC(nodei+3) = magnitude of 

U(nodei+3) velocity vector, squared, divided 
by the magnitude of the position vector of 
U(nodei). 

WVNLC(nodei+4) = magnitude of 
U(nodei+ 3) velocity vector, squared, 
divided by the magnitude of the position 
vector of U(nodei). 

WVNLC(nodei+5) = magnitude of 

U(nodei+3) velocity vector, squared, divided 
by the magnitude of the position vector of 
U(nodei). 

WVNLC(nodei+6) = absolute value of 
U(nodei+6), spacecraft mass. 

IPOST can interpolate a value for WVNLC if it 
lies between two input points. 

The weighting for the subproblem non-linear 
constraints. 


WVUS REAL 500*1. DO The independent variable weighting for 

INDSVR. 
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VARIABLE 

NAME 

WVU 


STORED 

TYPE VALUE DESCRIPTION 

REAL 25*1.D0 The independent variable wei ghting for 

INDVR(iyU(i). Recommended value is WVU(i) = 
abs(U(i)). Collocation: IPOST will use a set of 
default values of WVU for the defects and event 
discrepancies if the user does not input values. 
The defaults are: 

WVU(nodei+0) = magnitude of U(nodei) position 
vector. 

WVU(nodei+l) = magnitude of U(nodei) position 
vector. 

WVU(nodei+2) = magnitude of U(nodei) position 
vector. 

WVU(nodei+3) = magnitude of U(nodei+3) 
velocity vector. 

WVU(nodei+4) = magnitude of U(nodei+3) 
velocity vector. 

WVU(nodei+5) = magnitude of U(nodei+3) 
velocity vector. 

WVU(nodei+6) = absolute value of U(nodei+6), 
spacecraft mass. 

IPOST can interpolate a value for WVU if it lies 
between two input points. 


5.2 


The following is a description of input parameters for Namelist $TRAJ. $TRAJ data 
sets correspond to as many events as the user desires. REAL refers to REAL*8, 
LOG refers to LOGICAL, INT refers to INTEGER*4, and CHAR refers to character 
variables. 


VARIABLE STORED 

NAME TYPE VALUE DESCRIPTION 

ALT(50) REAL VARIOUS Tabular atmospheric altitude data 

(km) used to calculate density and 
speed of sound at a given altitude. 
Elements in the array correspond to 
atmospheric pressure (PRES), 
molecular we ight ( WGT), and 
temperature (TPT) for IATMOS = 1. 
Table size is specified by ITABSZ. 

ALTATM(22) REAL 22*0. DO Altitude boundary of sensible atmosphere 

ALTATM (0) = Sun's altitude 

(1) = Mercury's altitude 

(2) = Venus’ altitude 

(3) = Earth’s altitude 

(4) = Mars’ altitude 

(5) = Jupiter's altitude 

(6) = Saturn's altitude 

(7) = Uranus’s altitude 

(8) = Neptune's altitude 

(9) = Pluto’s altitude 

(10) « Earth’s moon altitude 
ALTATM (11) through (22) reserved for 
other bodies’ ALTATM. 

ALTIT REAL 0. Circular altitude of park orbit (see 

MANTYP and ILNCH) 
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VARIABLE STORED 

NAME TYPE VALUE 


DESCRIPTION 


ATMS (10) REAL Atmospheric initial data. 

(Defaults are defined for Earth's atmospheric 
data.) 

1.225D0 ATMS(1) = Base density for exponential and 

tabular atmospheres. (Default is for tabular - 
kg/m .) 

1.01325D5 ATMS(2) = Base pressure for tabular 
atmosphere, (newton/m ). 

28.964D0 ATMS(3) = Base molecular weight for tab ular 
atmosphere, (kg/mole) 

3.139D-7 ATMS(4) = First order gravitational constant for 
tabular atmosphere. (1/m) 

287.D0 ATMS(5) = Specific gas constant for tabular 
atmosphere. (joules/kg/K) 

7 . 162D0 ATMS(6) = Scale height for exponential 

atmosphere, (km) 

0.D0 ATMS(7) = Base altitude for exponential 

atmosphere, (km) 

1.4.D0 ATMS(8) = Ratio of specific heat 

0.D0 ATMS(9) = T JJ). 

0.D0 ATMS(10) = T.B.D. 

BTHETA REAL 0. B-plane theta angle for escape (see 

MANTYP and ELNCH) 

BTHRST REAL l.D0,2*0.D0 Thrust alignment unit vector in the body 

system. 

C(i j,k) REAL VARIOUS The analytic elements of the planetary 

ephemeris. 

i - polynomial coefficient; third order in time; 
i = 1, 2, 3, 4 

j - planet number; 0=sun, 

3=earth,...; j * 0, 1, 2,..., 22 
k - element type: 

1 = semimajor nxia 

2 = eccentricity 

3 = inclination 

4 = longitude of ascending node 

5 = argument of periapsis 

6 = mean anomaly 


CD 

REAL 0.D0 

Drag coefficient for the vehide(can be 
overridden by $TAB input). 

CL 

REAL 0.D0 

Lift coefficient for the vehicle (can be overridden 
by $TAB input). 
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VARIABLE STORED 

NAME TYPE VALUE DESCRIPTION 

COINTP CHAR, linear’ How initial values of ordinary nodes are 

calculated from event nodes, by linear/hermite 
interpolation. 

linear' linear interpolation 
'hemite' hemite interpolation 

CONODE CHAR. '5' Defines a property of collcoation internal nodes. 

= 'float', the internal nodes float so that any 
variable length is divided into nsegph 
segments of equal time duration. 

= 'fixed', the internal nodes are fixed in time. If 
any varaible length phase increases or 
desreases in duration , then the last segment 
absorbs the entire increase or decrease. 

CRITR CHAR. TIME ’ The name of the event criterion variable, (see 

Table 3-9) 

C3 REAL 0.D0 Escape energy, (see MANTYP and ILNCH) 

DATE(i) REAL 0.D0 Calendar Date. 

(Y ear,Month,Day, Hour .Minute .Seconds ) 

DEC REAL 0.D0 Declination of launch asymptote, (see 

MANTYP and ILNCH) 

DT REAL 1.D0 Constant delta-time propagation step for the 

Multiconic, Cowell, and Encke propagators. 

DTIMR(20) REAL 20*0.D0 Reference time for activities. Initiated by 

setting DTIMR(i) = 1, Corresponding event 
criteria times are: 

TIMRF1,TIMRF2,...,TIMRF20. 

DTIMR(1) = Start of mission (first event) 
DTEMR(2) = Start of finite burn 
DTIMR(3) = Nuclear electric decay start 
DTIMR(4) = Blowdown propulsion initiation 
DTIMR(5) = Throttle level polynomial start 
DTIMR(6) = Body roll angle polynomial start 
DTIMR(7) = Body yaw angle polynomial start 
DTIMR(8) = Body pitch angle polynomial start 
DTIMR(9) * Not available 
DTIMR (10) = Not available 
DTEMR (11) through DTIMR(20) are vavailable for 
user specification. 
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VARIABLE STORED 

NAME TYPE VALUE DESCRIPTION 

DVX(3) REAL 3*0. DO Delta velocity vector to be added to the state at 

the current maneuver. Use MANTYP = 
'IMPULS'. 

EVENT(l) REAL 0.D0 Event Sequence Number for the current, phas e. 


EVENT(2) REAL 0.D0 Primary/Roving flag. 

= 0, Primary. 

= 1, Roving. 

FRAC REAL 0.0 The dry jettison mass function, evaluated 

according to input flag, JETT. For JETT: 

= 2, the dry jettison weight is equal to FRAC. 

= 3, the diy jettison weight is equal to 
WPROP (1/FRAC - 1). 

If Table FMASST is input, FRAC is overridden 
by the table lookup value. 

FRACI REAL 0.0 The initial fraction of total spacecraft mass, 

SCMASS, which is weight of propellant, 
WPROP. If Table FMASIT is input, FRACI is 
overridden by the table lookup value. 

GJ2 REAL VARIOUS J2 terms for the 

GJ2(0) = Sun’s GJ2. 

GJ2(1) = Mercury's GJ2. 

GJ2(2) = Venus's GJ2. 

GJ2(3) = Earth's GJ2. 

GJ2(4) = Mars' GJ2. 

GJ2(5) = Jupiter's GJ2. 

GJ2(6) = Saturn's GJ2. 

GJ2(7) = Uranus's GJ2. 

GJ2(8) s Neptune's GJ2. 

GJ2(9) = Pluto's GJ2. 

GJ2(10) s Earth's moon GJ2. 

12*0.D0 GJ2(11) - (22) Reserved for other bodies' GJ2. 
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VARIABLE STORED 

NAME TYPE VALUE 


DESCRIPTION 


GMU REAL VARIOUS Gravitational constants, (u), for planetary 

bodies. 

GMU(O) = Sun’s GMU. 

GMU(l) = Mercury's GMU. 

GMU(2) = Venus’s GMU. 

GMU(3) = Earth's GMU. 

GMU(4) = Mars’ GMU. 

GMU(5) = Jupiter's GMU. 

GMU(6) = Saturn's GMU. 

GMU(7) = Uranus's GMU. 

GMU(8) - Neptune's GMU. 

GMU(9) = Pluto's GMU. 

GMU(IO) = Earth's moon GMU. 

12*0.D0 GMU(ll) - (22) Reserved for other bodies' GMU. 

GO REAL 9.80665D-3 Sea level gravitational constant, (km/sec ) 


IATMOS INT. 1 Flag which specifies the type of atmospheric 

model to be used. 

= 1, user defined exponential atmospheric 
model. 

= 2, user defined tabular atmosphere data. 

IBODY INT. VARIOUS Four digit planetary code. First two digits 

correspond to the satellites of the planets; Last 
two digits correspond to the planet number. 
EBODY(O) = 0000 (Sun). 

IBODYd) = 0001 (Mercury). 

IBODY(2) = 0002 (Venus). 

IBODY(3) = 0003 (Earth). 

IBODY(4) = 0004 (Mars). 

IBODY(5) = 0005 (Jupiter). 

IBODY(6) = 0006 (Saturn). 

IBODY(7) = 0007 (Uranus). 

IBODY(8) = 0008 (Neptune). 

IBODY(9) * 0009 (Pluto). 

IBODY(10) = 0103 (Earth's moon). 

12*0.D0 IBODY(ll) - (22) Reserved for other bodies' 
identifiers. 
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VARIABLE 

NAME 

ICBODY 


STORED 

TYPE VALUE DESCRIPTION 

INT -1 Reference Body for collocation node inputs (U's) 

= - 2, use IPBODY (1) of last phase 
= - 1, use IPBODY (1) for the phase 
= i > 0, use Body i, i = 0 sun, i = 3 Earth, etc. 
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VARIABLE STORED 

NAME TYPE VALUE 


DESCRIPTION 


ICOORD INT. 1 Flag which specifies the coordinate frame in 

which S/C body orientation is defined. 

= 1, inertial ecliptic system. 

= 2, inertial planet equatorial system. 

= 3, UVW system. 

= 4, RTN system 
= 5 B-plane system 
= 6 Cone and dock system 


IDBODY INT. 1 Index into IBODY array for the planetary body 

of the state vector. 

IDFRAM(l) CHAR. 'ECLIPTIC' Frame of reference for input and output state 

vector. 

= 'ECLIPTIC', Ediptic of 1950. 

= 'EARTHEQ', Earth equatorial. 

= 'BODYEQ', Body equatorial. 

IDFRAM(2) CHAR. 'MEAN1950' Epoch of state vector. 

= 'MEAN 1950', Mean 1950 
= 'MEAN2000', Mean 2000 
= ’MEANDATE’, Mean of date. 

IDT INT. 1 Flag defining the step size for Multiconic and 

Encke propagators. 

= 1, constant step size (DT). 

* 1, calculate step size depending on maximum 
step size (DT) and scaling(STEP). 

IEPOCH CHAR. 'JULIAN' Date input type flag. 

= 'CALEND' is calendar date. 

= 'JULIAN is Julian date. 
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VARIABLE STORED 

NAME TYPE VALUE 


DESCRIPTION 


IFORCE (5) INT. 5*0 


ILNCH INT. 0 


INC REAL 0.D0 


Array of perturbing force flags. 

= 0, off. 

* 0, on. 

IFORCE(l) = number of perturbing bodies. 
IFORCE(2) = J2 force for central body (flag). 
EFORCE(3) = solar radiation pressure (flag). 
IFORCE(4) = atmospheric lift and drag (flag). 
IFORCE(5)= thrust acceleration 

= 1, generalized thrust profile (Table input) 
= 2, blowdown system 
— 3, electric (nuclear or solar) propulsion 
If IFORCE (1) is non-zero, input NPERT. 

Flag for initializing state, assumes launch from 
park orbit, (requires ALTIT = park orbit altitude 
INC = hyperbola inclination 
BTHETA = B-plane angle used if INC < O) 
= 0, no launch inputs 
= 1, input C3, RA, DEC 
= 2, input VINFX, VINFY, VINFZ 

Orbital inclination, (see MANTYP) 
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VARIABLE STORED 

NAME TYPE VALUE 

INPUTX CHAR. ’NONE’ 


IPBODY(i) INT. 0,3 


IPROP CHAR. '1STEF 


IPBODY(l) 

IPBODY(l)). 

ITABSZ INT. 22 


DESCRIPTION 

Mode of state input/update. 

= ’NONE', no state input/update. 

(reset after each $TRAJ) 

= ’CARTES', input X as cartesian state. 

= 'CONIC', input X as conic state. 

Index into IBODY array for primary and 
secondary bodies for propagator. (See IPROP) 
IPBODY(l) = Primary body index. 

IPBODY(2) = Secondary body index. 

= 0, Sun 
= 1, Mercury 
= 2, Venus 
= 3, Earth 
= 4, Mars 
= 5, Jupiter 
= 6, Saturn 
= 7, Uranus 
= 8, Neptune 
= 9, Pluto 
= 10, Earth's moon 
= 11 - 22, other user defined bodies. 

Propagator mode. 

= '1STEP', Onestep. (requires 
IPBODYd) and IPBODY(2)) 

= 'CONIC', 2-Body Conic, (requires IPBODYd)). 
= MULTIC', Stumpf Multiconic (requires 
IPBODY(l),IPBODY(2) and NPERT)). 
s 'COWELL', Cowell propagator, (requires 

= ’ENCKE’, Encke propagator, (requires 


Size of input atmospheric table parameters 
(ALT, PRES, TPT, WGT.) 
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VARIABLE STORED 

NAME TYPE VALUE 

JETT INT 0 


JULDAT REAL 0.0 
MANTYP CHAR. ‘NONE’ 


MASSFI CHAR 'NONE’ 


DESCRIPTION 

Mass fraction jettison calculation flag. 

= 0, do nothing at this event. 

= 1, jettison previously calculated and saved 
value of jettison_dry_weight. 

= 2, set WJETTM * FRAC, or table 
lookup of FMAS ST, an d save it to be 
jettisoned later (JETT = 1). 

= - 2, set WJETTM = FRAC, or table lookup of 
FMASS T, and jettison at this phase. 

= 3, set WJETTM = WPROP * (1/x - 1) and save 
it to be jettisoned later, where x = FRAC or 
table lo okup of FMASST. 

= - 3, set WJETTM = WPROP * (1/x - 1) and 
jettison at this phase., where x = FRAC or 
table lookup of FMASST. 

JETT is reset to 0 after each phase 

initialization. 

The Julian date in days. 

Type of impulsive maneuver and subsequent 

state initialization. 

= ‘NONE’, no impulsive maneuver 

= ‘LAUNCH’, launch event (requires ILNCH). 

= TMPULS’, generalized impulsive AV 
(requires DVX, DVY, DVZ) 

= ‘ORB INS’, orbit insertion (requires RAPOAP, 
RPERI, INC of desired final orbit) 

Initialize mass flag. 

= NONE', do not initialize weight this phase. 

= WPROP=', set WPROP = SCMASS * FRACI. 

= 'SCMASS = set SCMASS = WPROP/FRACI. 
MASSFI is reset to NONE' after each phase 
initialization. 


MDL INT. 1 The event cycling/phasing model. 

= 1, hit VALUE exactly. 

= 2, hit VALUE, positive derivative. 
= 3, hit VALUE, negative derivative. 
= 5, mdl=l with 360 discontinuity. 

= 6, mdl=2 with 360 discontinuity. 

= 7, mdl=3 with 360 discontinuity. 

= 8, hit if greater than VALUE. 

= 9, hit it less than VALUE. 
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DESCRIPTION 


VARIABLE 

NAME 

NAMLST 


NPERT 


NSEGPH 

NSGPHO 

NSGPWD 


NSGPWI 


STORED 
TYPE VALUE 

CHAR. ’TRAJ’ 


INT. 21*0 


INT. 0 
INT. 0 


INT. -1 


INT. -1 


Namelist flag for next input. 

= 'NONE', no more namelist, solve problem. 

= TABLE', read namelist TAB next. 

= TOP', solve problem, then read namelist TOP 
to start the next problem. 

= TRAJ', read namelist TRAJ next. 

NOTE: TRAJ and TABLE namelists can be 
read only in the first problem. 

Index into IBODY array for the perturbing 
bodies. (IFORCE(l) > 0). 

(IPROP=’MULTIC', ’COWELL’, or 'ENCKE'). 

= 0, Sam 
= 1, Mercury 
= 2, Venus 
= 3, Earth 
= 4, Mars 
= 5, Jupiter 
= 6, Saturn 
= 7, Uranus 
= 8, Neptune 
= 9, Pluto 
= 10, Earth's moon 
= 11 - 22 other user defined bodies. 

Number of collocation segments for this phase. 

Number of collocation input values of controls 
(U) for this phase. NSGPHO = 0 means to 
interpolate between events to obtain internal 
node states. 

Number of collocation input values of WVNLC 
for this phase. NSGPWD has a role similar to 
NSGPHO. The default for NSGPWD = 
NSGPHO. NSGPWD may be set = 0 for any 
phase. 

Number of collocation input values of WVU for 
this phase. NSGPWI has a role similar to 
NSGPHO. The default for NSGPWI is 
NSGPWI = NSGPHO. NSGPWI may be set = 0 
for any phase. 
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VARIABLE STORED 

NAME TYPE VALUE 


DESCRIPTION 


NT1MES INT. 0 Number of repetitions of a repeating roving 

event. 

< 0, Repeat the event NTIMES using option 2 of 

ROVET. 

= 0, Do not repeat the event. 

> 0, Repeat the event NTIMES using option 1 of 

ROVET. 

PDUMP REAL 0.D0 Fraction of remaining WPROP to be added to 

jettison mass (WJETT). 

PINC REAL 1.0D6 Print interval for the trajectory print block. 

PITCHO, REAL 3*0.D0 Coefficients for cubic time dependent motion 

PITCH1, of pitch angle. 

PITCH2, 

PITCH3 

PLANET CHAR VARIOUS Planet names associated with IBODY. 

POLEV(ij,k) REAL VARIOUS Coordinates of polar aria, Earth Equato rial 

Mean 1950. 

i - polynomial coefficient, third order in time; 
i = 1,2,3, 4 

j - planet number, 0 = sun, ...3 = Earth,... 
k - angle type; 1 = right ascension, 

2 = declination 

PRES (50) REAL EARTH Tabular pressures for atmosphere corresponding 

to ALT. (newtons/m 2 ) 

PRNC REAL -1. Output interval for PROFIL data file. 

< 0., no PROFIL output 

= 0., all events and propagation steps. 

> 0., output interval (days). 
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VARIABLE STORED 

NAME TYPE VALUE 


DESCRIPTION 


PRPDAT REAL 15*0.D0 Data for pressure regulated and blowdown 

propulsion systems. 

PRPDAT(l) = Not input. Internally set equal to 

Table input ($TAB) Vacuum thrust for 

generalized thrust propulsion. 

PRPDAT(2) = Exhaust area for generalized 
thrust propulsion. 

PRPDAT(3) = Maximum mass flow rate for 
blowdown propulsion. 

PRPDAT(4) = Propellant density for blowdown 
propulsion. 

PRPDAT(5) = Initial tank volume for blowdown 
propulsion. 

PRPDAT(6) = Ullage ratio for blowdown 
propulsion. 

PRPDAK7) = Ratio of specific heats for 
blowdown propulsion. 

PRPDAT(8) = Initial thrust for blowdown 
system. 

PRPDAT(9) = Thrust efficiency for electric 
propulsion. 

PRPD AT( 10) = Housekeeping power for electric 
propulsion. 

PRPDAT (11) = Maximum power for electric 
propulsion. 

PRPDAT (12) = decay constant for nuclear 
electric. 

PRPDAT (13) to (15) = Solar array constants 
for solar electric propulsion. 


PSCALE 

REAL 

1 . 

Propellant consumption effectiveness, 
(rate scale factor). 

RA 

REAL 

0.D0 

Right ascension of launch asymptote 
(see MANTYP and ILNCH). 

RAPOAP 

REAL 

0.D0 

Desired apoapsis radius. 
(MANTYP = ‘ORBINS’) 
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VARIABLE 

NAME 

RE 


RECT 


ROLLO, 

ROLL1, 

ROLL2, 

ROLL3 

ROVET 


STORED 

TYPE VALUE DESCRIPTION 


REAL VARIOUS Planetary radii. 

RE(0) = Sun's radius. (695950.) 

RE(1) = Mercury's radius. ( 2433.) 

RE(2) = Venus's radius. ( 6053.) 

RE(3) = Earth’s radius. ( 6378.) 

RE(4) = Mars' radius. ( 3386.) 

RE(5) = Jupiter's radius. ( 71370.) 

RE(6) = Saturn's radius. ( 60369.) 

RE(7) = Uranus’s radius. ( 24045.) 

RE(8) = Neptune's radius. ( 22716.) 

RE(9) = Pluto's radius. ( 5700.) 

RE(10) = Earth's moon radius ( 1738.) 

12*0.D0 RE(ll) - (22) Reserved for other bodies' radii. 


REAL 0.01 


Encke propagation rectification threshold, 
lAR/Rol. 


REAL 3*0.D0 


Coefficients for quadratic time dependent 
motion of roll angle. 


REAL 2*0.D0 Values for repeating roving events. 

Option 1: ROVET(l)=x, ROVET(2)=dx, where 
x = value of CRlTK on first occurrence of the 
event, and x + n*dx is the value of CRTTR on 
the n- the repetition of the event. 

Option 2: ROVET(l)=xl, ROVET(2)=x2, 

ROVET(N TIMES+ 1 )=xn , where xi = value of 
CRITR for the i-th occurrence of the event. 
Note: Thi s option is lim ited to nine repetitions. 
(See NTIMES: if NTIMES > 0, Option 1 is 
selected if NTIMES < 0, Option 2 is selected). 
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VARIABLE STORED 

NAME TYPE VALUE DESCRIPTION 


RSOI REAL VARIOUS 




12*0. DO 

SCDRAG 

REAL 

3*0.D0 

SCMASS 

REAL 

1.D6 

SCSFA 

REAL 

3*0.D0 


Radius of sphere of influence for each planetary 
body. 

RSOI(O) = Sun’s RSOI. 

RSOI(l) = Mercury’s RSOI. 

RSOI(2) = Venus’s RSOI. 

RSOI(3) = Earth’s RSOI. 

RSOI(4) = Mars' RSOI. 

RSOI(5) = Jupiter’s RSOI. 

RSOI(6) = Saturn’s RSOI. 

RSOIC7) = Uranus's RSOI. 

RSOIC8) = Neptune’s RSOI. 

RSOI(9) = Pluto’s RSOI. 

RSOI(IO) = Earth's moon RSOI. 

RSOI(ll) - (22) Reserved for SOI of other bodies. 

S/C areas normal to the three body axes used in 
calculating atmospheric effects. 

Vehicle mass instantaneous update. 


S/C areas normal to the three body axes, 
includes reflectivity coefficients for use in 
calculating solar radiation pressure. 
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VARIABLE STORED 

NAME TYPE VALUE 


DESCRIPTION 


SFC 

REAL 

1.024D17 

Solar flux constant, kg-m/s 2 . 

SPI 

REAL 

1.D6 

Engine specific impulse. 

STEP 

REAL 

1.D0 

Step size scaling factor for Multiconic and 
Encke propagators. 

THEDOT 

REAL 

VARIOUS 

Planetary rotation rate, (rad/sec) 
THEDOT(O) = Sun’s rotation rate. 


THEDOT(l) = Mercury' 8 rotation rate. 
THEDOT(2) *= Venus's rotation rate. 
THEDOT(3) = Earth's rotation rate. 
THEDOT(4) = Mars' rotation rate. 
THEDOT(5) = Jupiter’s rotation rate. 
THEDOT(6) = Saturn's rotation rate. 
THEDOT(7) = Uranus's rotation rate. 
THEDOT(8) = Neptune's rotation rate. 
THEDOT(9) = Pluto' 8 rotation rate. 
THEDOT( 10) = Earth's moon rotation rate. 


12*0. DO THEDOT(ll) = (22) Reserved for other bodies' 

rotation rates. 


THLO, 

THL1, 

THL2, 

THL3 

REAL 

3*0.D0 

Coefficients for quadratic time 
dependent setting throttle level. 

THRUST 

REAL 

1000. 

Average thrust level for impulsive maneuver. 

TOL 

REAL 

1.0D-12 

Hybrid tolerance for event cycling, 

TOL = (vd-vay(l+abs(vd)), where va = actual 
value of the CRITR var. vd = VALUE. 

TPT(50) 

REAL 

VARIOUS 

Tabular data for atmospheric molecular 
temperatures corresponding to ALT. (K) 
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VARIABLE STORED 

NAME TYPE VALUE DESCRIPTION 

TSOI REAL VARIOUS Time of sphere of influence in days for each 

planetary body, (requires IPROP = ‘1STEP’) 
TSOI(O) = Sun's TSOI. 

TSOI(l) = Mercury’s TSOI. 

TSOI(2) = Venus's TSOI. 

TSOI(3) = Earth’s TSOI. 

TSOI(4) = Mars' TSOI. 

TSOI(5) = Jupiter’s TSOI. 

TSOIC6) = Saturn’s TSOI. 

TSOIC7) = Uranus's TSOI. 

TSOIC8) = Neptune’s TSOI. 

TSOI(9) = Pluto's TSOI. 

TSOI(IO) = Earth’s moon TSOI. 

12*0.00 TSOI(ll) - (22) Reserved for other bodies. 

VALUE REAL 1.0D10 The desired value of the criterion variable 

(CRITR) at which the event is to occur. 

VINFX, REAL 3*0.D0 Hyperbolic excess velocity Qaunch 

VINFY, V-infinity) (see MANTYP and ILNCH) 

VINFZ 


WGT(50) 

REAL 

EARTH 

Tabular values for atmospheric molecular 
weights corresponding to ALT. (kg/mole) 

WJETT 

REAL 

0.D0 

Jettison mass for this event, (kg) 

WPROP 

REAL 

1.D6 

Propellant mass, (kg) 

X(i) 

REAL 

6*0 .DO 

State vector to be input/updated for 
NPUTX='CARTES’, 

X=X,Y,Z,VX,VY,VZ; for 
INPUTX=’CONIC’, X=( semi-major-axis, ecc., 
inc., long.-of-asc.-node, arg.-of-periapsis, mean- 
anomaly). 

YAWO, 

YAW1, 

REAL 

3*0.D0 

Coefficients for quadratic time dependent 
motion of yaw angle. 


YAW2, 

YAW3 
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5.3 


The following is a description of input parameters for namelist $TAB. The $TAB 
data must follow the $TRAJ event data which is applicable. Currently, only four 
parameters can be input as tables: thrust vs. time, drag coefficient vs. angle of 
attack and mach number, lift coefficient vs. angle of attack and mach number and 
S/C mass vs. C3. 

VARIABLE STORED 

NAME TYPE VALUE DESCRIPTION 

DCTRP INT. 2*1 A flag to activate the no-extrapolation 

i = 1,2 feature for argument values that exceed 

the table boundaries in x and y, respectively. 
= 0, do not extrapolate. Use the last table 
value instead. 

= 1, Extrapolate to obtain the table value. 

NAMLST CHAR ‘TRAJ 7 Namelist flag for next input. 

= *NONE\ no more namelists, end of 
input stream. 

= TOP, $TOP will follow. 

= TRAJ\ $TRAJ will follow 
= TABLE’, $TAB will follow 

TABLE ( ) REAL 0.D0 Table input array (see below). 

The table inputs for IPO ST are generalized to include: 

1) Allowable size for each table of 2500 entries. The total size of all tables is 
limited by the amount of storage allocated to IBKT. Both of these values 
can be changed by a simple program modification to satisfy user 
requirements. 

2) Generalized argument specification. The argument to be used for each 
table is specified by input and can be any computed output variable. 

3) Constant-valued, monovariant or bivariant table types. 

4) Linear interpolation capability. 

5) A no-extrapolation feature for argument values that exceed the table 
boundaries. This feature is obtained by inputting the flag IXTRP (i) to 
request no-extrapolation in x and y respectively, as follows: 

=0, Do not extrapolate. Use the last table value instead, 
si, Extrapolate to obtain the table value. 

Each table is input in namelist $TAB as the array TABLE. As a result, each 
table being input requires a separate input of namelist $TAB. 
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These tables are available in IPOST. 


CDT is the drag coefficient (will override CD in $TRAJ) 

CLT is the lift coefficient (will override CL in $TRAJ) 

FMASIT is the initial mass fraction defining WPROP as a fraction of SCMASS, 
or vice versa. FMASIT overrides FRACI. 

FMASST is the jettison mass fraction table. JETT is the jettison . FMASST 
overrides FRAC. 

SMA SST is the vehicle mass (will override SCMASS in $TRAJ) 

THRSTT is the vacuum thrust in the generalized thrust model IFORCE(5) = 1 
in $TRAJ 

The choices for independent parameters) are defined in Table 1-5. 

The elements of the array TABLE are different depending on the type of table 
being input. 

Table pointers are used by the table look-up routines for all tables (except 
constant-valued tables) to provide efficient operation of the table look-up routines. 
These pointers should always be input as 1. The number of pointers for a given 
table depends on the table type, i.e., monovariant or bivariant. For example, a 
monovariant table has one pointer which is input as TABLE (7), and a bivariant 
table has five pointers input as TABLE (10) through TABLE (14). 




The table elements for constant-valued tables are: 

TABLE(l) The Hollerith name of the table. 

TABLE(2) The table type. 

=0, Constant-valued table. 

TABLE(3) The table value. A typical constant-valued table would be input 
as follows: 

TABLE = 6HTHRSTT,0,2000000., 
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Monovariant tables are formed by a series of ordered pairs of x and f(x) which 
represent a curve that is a series of line segments. 

The size of a monovariant table is obtained by: 

n = 7 + 2 * n x 

The table elements for monovariant tables are: 

TABLE(l) The Hollerith name of the table. 

TABLE(2) The table type. 

=1, Monovariant table. 

TABLE(3) The Hollerith name of the x argument. 

TABLE(4) The number of x values in the curve, n x 

TABLE(5) The interpolation type. 

=0, Step table, i.e., no interpolation. 

=1, Linear interpolation. 

TABLEC6) The type of x values. 

=-l, Decreasing values. 

= 1, Increasing values. 

TABLE(7) The initial value of the x pointer. 

TABLE(8) The first table value (xl). A typical monovariant 
table would be input as follows 

TABLE = 3HCDT,1,6HSMACH ,n x , 1,1,1, 
xi,flxi), X2,f(x2), ...,xn x ,fixn x ), 
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Bivariant tables are formed by a family of monovariant curves, where a 
monovariant curve is input for each value of y. As a result, the arguments are x and 
y and the function is flx,y). The x arguments in each curve do not need to be the 
Mmft value; however, there must be the same number of x values in each curve. 

The size of a bivariant table is obtained by: 

n — 14 + ny * (2 * n* + 1) 

The table elements for bivariant tables are: 

TABLE(l) The Hollerith name of the table. 

TABLE(2) The table type. 

=2, Bivariant table. 

TABLE(3) The Hollerith name of the x argument. 

TABLE(4) The Hollerith name of the y argument. 

TABLE(5) The number of x values in each curve, n x 

TABLE(6) The number of y values (curves) in the family , n y 

TABLE (7) The interpolation type 

si, Linear interpolation. 

TABLEC8) The type of x values 

=-l, Decreasing x values. 

= 1, Increasing x values. 

TABLEC9) The type of y values 

=-l, Decreasing y values. 

= 1, Increasing y values. 

TABLE(IO) The initial value of the yj curve pointer (input = 1). 

TABLE(ll) The initial value of the x pointer on the yj curve (input = 1). 

TABLE(12) The initial value of the x pointer on the yj+1 curve (input = 1). 

TABLE(13) The initial value of the x pointer on the yj+2 curve (input = 1). 

TABLE(14) The initial value of the x pointer on the yj+3 curve (input = 1). 
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TABLEU5) 


Bivariant Tables (Continued) 

The first table value (yj). A typical bivariant table would be 
input as follows: 

TABLE = 6HCDT ,2,6HSMACH ,6HANGAT1, n x ,n y 

l f l,U,l, 1,1,1, 

yi, xi, Kxi,yi), X 2 , feyi), —»xn x , flxn x yi> 
y2, xi, fixi,y 2 >' X 2 , ftx2y2), -,xn x ,^xn x y2)> 

yn y , xi, Rxi, yn y ), X 2 , fix 2 , yn y ), xn x t ^ x n x ,yn y ) > 
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6.0 OUTPUT DESCRIPTION 

6.1 STANDARD OUTPUT 


During IPOST execution, a variety of output is made available to the user. A 
number of print flags, such as IPRINT, control the degree and type of output. These 
flags are discussed in Section 6.3. This section will describe the standard default 
output. 

The first output which appears during execution is an echo of the namelist inputs, 
$TOP, and all the STRAJ’s. After the namelist echo, there will be an input 
summary of what events have been chosen along with the criteria, the initial 
trajectory inputs and the subproblem setup, if any. An NPSOL summary, the 
master problem controls, constraints and optimization will also be displayed. 

For explicit optimization using decomposition, IPOST starts to solve the master 
problem by first solving each subproblem in consecutive order. The first subproblem 
is iterated until the desired target conditions are readied. A summary of results is 
displayed. After all subproblems have been solved, with their results displayed, the 
master problem can proceed to optimize the cost (objective) function. Each master 
problem control parameter is perturbed (assuming gradients are computed by finite 
differencing, and not analytically), and the sequence of subproblems is repeated, in 
order to evaluate cost and constraint parameter values. The subproblem solutions 
are not normally displayed for perturbed master problem trajectories. The objective 
gradient and constraint gradients are displayed for the spedfied master problem 
control parameters. 

The master problem now determines the appropriate master problem control 
parameter step size and direction in order to optimize cost and maintain 
constraints. For the new set of control parameters, the subproblem sequence is 
solved once again and the master problem iteration is completed when a summary 
is displayed. 

When the cost function cannot be optimized further and all constraints are satisfied, 
the master problem concludes with a convergence summary. If a PROFIL file has 
been generated, then a parameter list is displayed. 

The final optimized trajectory is now displayed. The trajectory time history display 
consists of a print block at each event (start of a phase), at user spedfied print 
intervals during a phase, and at the end of a phase (next event). 

Certain activities will trigger additional print data. For example, when propulsion 
has been activated, a print block assodated with propulsion data and a print block 
assodated with vehide orientation are displayed. Trajectory displays can also occur 
at other stages of the IPOST solution process. These occurrences are controlled by 
IPRO ($TOP) and IPRINT ($TOP). 
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IPRINT is the general control flag. As the value IPRENT increases, the volume of 
output increases. For example, IPRINT = 0 will print only cursory data for the run, 
whereas IPRINT = 5 will print massive amounts of debug printout. 

6.2 PROFIL OUTPUT 

A binary output file (PROFIL) can be written. It is activated by setting a non- 
negative value of PRNC in $TRAJ. Additional inputs required are, NPVAR and 
OUTNAM. NPVAR specifies how many variables are to be in each PROFIL block, 
and OUTNAM lists the specific variable names. For example, PRNC = 100, NPVAR 
= 2 and OUTNAM = TIMRF’, DVMAG ’ will produce mission time (days) and AV 
mag nitude in each PROFIL block starting with the current event, and then every 
100 days thereafter. Any variable in the data dictionary, corresponding to Tables 3 
- 8 and 3-9, can be used in OUTNAM. Up to 100 parameters can be included in a 
PROFIL block. 


6.3 DEBUG PR 




There are several debug print flags in EPOST. The first is the generalized debug 
flag IPRINT. IPRINT can have a value from 0-4, with the amount of debug print 
increasing as the value of IPRINT increases. With IPRINT = 0, there is no debug 
print. With IPRINT = 4, a horrendous amount of printout will occur from trajectory 
data, to subproblem sensitivity matrices, to entry and exit messages from different 
subroutines. More will be added in later revisions of this document as to what each 
level of IPRINT will supply the user. IPRINT can be set at each event, so debug can 
be performed only in problem areas. 

The debug print flag, IPRO, determines the amount of trajectory data the user sees. 
If IPRO = -1, the default, the user will only see the final targeted and optimized 
trajectory. If IPRO = 1, the original nominal trajectory will be printed, along with 
the final targeted and optimized trajectory. If IPRO = 1, all the nominal trajectories 
will be printed along with the final. Finally, if IPRO = 2, every trajectory performed 
will be printed. 

IDEB is a trial step summary debug print flag. When NPSOL is performing trial 
steps calculating the sensitivity matrices required, IDEB is defaulted to 0 and the 
user will not see these calculations. If the user would like to see these steps, IDEB 
may be set equal to one. 
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6.4 ERROR MESSAGES 

The error messages are presented in alphabetical order as follows: 

Message : BAD IEPOCH VALUE 
SourceATvpe: DYNXMI/fatal 

Unnditinn/Cnrrertive Action: Invalid value for IEPOCH. Possible values are 
’CALEND’, 'JULIAN', or 'NOINPUT. Check input data. 

Message: BAD IDFRAM VALUE 
Source/Tvoe: DYNXMI/fatal 

Condition/Corrertive Action: Invalid value for IDFRAM. Possible values are 
•ECLIPTIC' or 'EQUATOR'. Check input data. 

Message : BAD INPUTX VALUE 
SourceAFvpe: DYNXMI/fatal 

Condition/Corrertive Action: Invalid value for INPUTX. Possible values are 
'CONIC', 'CARTES', 'DELTAV, or NONE’. Check input data. 

Message : CNFUNC, UNUSABLE NOMINAL 
Source/Typei CNFUNC, warning . 

Condition/Corrertive Action: The trajectory went beyond realistic conditions, such 
as mass less than zero, and was prematurely terminated. Reconsider formulation of 
the problem. 

Message: cycxm problem 

Source/Tvp e: CYCXM/fatal 

Condition/Corrertive Action: The conditions for one of the pending events could not 
be met in 20 iterations. Check the event criteria and the integration step size 
inputs. 

Message: DINPT - ESN NOT FOUND 
Source/Tvpe: DINPT/fatal 

Condition/Corrective Action: The data for the event printed above this message as 
-ESN = XXX.XXX does not exist. 

Message: ERROR - NO PENDING EVENTS. DUMP OF IBKT 
Source/Tvne: TGOEMI/waming 

Condition/Corrective Action: No event criteria have been input. Check the value of 
each CRiTR variable. 

Message: EXCEEDED MAXIMUM NUMBER OF PENDING EVENTS. DUMP OF 
IBKT 

Source/Tvne: TGOEMI/waming 

Condition/Corrective Action: The number of events that have been specified as 
roving events exceeds the allowable maximum of 10. Reduce the number of roving 
events and rerun the job. 
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Message : *GRAD* FATAL TERMINATION 
Source/Tvpe: GRAD 2/fatal 

Condition/CoiTective Action: The indicated error condition caused the 
targeting/optimization algorithm to terminate the program execution. Check the 
problem setup and the inputs in namelist TOP. 


Message: *PAD* NOISE, INCREASE PERT 
Source/TVpe: PAD/ warning 

Condition/Corrective Action: An independent variable perturbation value, PERT or 
PERTSB, is too small, or the independent variable has no effect on the dependent 
variables when using the targeting/optimization algorithm. PAD will automatically 
adjust the pert and try the perturbed function evaluation again. If the results are 
unsatisfactory, rerun the problem with a larger input value of PERT/PERTSB. 


Message : *PAD*PATHOLOGICAL NOISE/PERTURBATION 
Source/Tvne: PAD/fatal 

Condition/Corrective Action: An independent variable perturbation causes noise 
and polarization in separate dependent variables. PAD cannot adjust the size of the 
perturbation. The problem may be pathological. Reformulate the problem. 

Message : ♦PAD* POLARIZATION, DECREASE PERT 
Source/Tvpe: PAD/warning 

Condition/Corrective Action: An independent variable perturbation value, PERT or 
PERTSB, is too large, or the independent variable has no effect on the dependent 
variables when using the targeting/optimization algorithm. PAD will automatically 
adjust the pert and tiy the perturbed function evaluation again. If the results are 
unsatisfactory, rerun the problem with a smaller input value of PERT/PERTSB. 

Message: READAT - CASE WON’T FIT 
Source/Tvne: READAT/fatal 

Cond l tiO D /CorrectlYe Act i on; The constant-valued general input data exceeds 
maxim u m number of cells. Reduce the amount of general data 


Message: READAT - FATAL INPUT ERRORS 
Source/Tvne: READAT/fatal 

Cond l tion/Correctiye Actaom The Hollerith input names shown are incorrect. Check 
all variables containing Hollerith inputs to see that the names are valid and spelled 
correctly. 


Message : RTOP - NAMELIST ERROR 
Source/Tvne: RTOP/fatal 

Condi tion/Corrective Action: A namelist error has occurred. Check input data. 

Message: RTRAJ - NAMELIST ERROR 
Source/Tvne: RTRAJ/warning 

Condition/Corrective Action: A namelist error has occurred. Check input data. 
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Message : RTRAJX, DATA BUFFER - EXCEEDED SIZE OF BUFFERS (GBKT OR 
IBKT) 


Source/Tvpe: RTRAJX/fatal , , 

Condi tinn/Cnrrertive Action: The constant-valued input data exceeds the m a xim u m 
number of cells. Reduce the amount of input data. 


Message : RTRAJX, NO. EVENT NO. WAS INPUT 
SourcaTvpe: RTRAJX/fatal 

Condition/Corrertive Action: A phase was input with no event number (EVENT), or 
it is zero. Input the event number as EVENT(l). 


Message : SERCH "name" 

Source/Tvpe: SERCH/fatal/warning 

Condition/Corrective Action: The character input variable "name" is not valid or is 
misspelled. This is a warning error only if the misspelled variable is a print 
variable. Check the character name in question. Message: Q-DYNXM - NO 
PROPAGATOR HAS BEEN CHOSEN Source/Type: DYNXM/waming 
Condition/Corrective Action: No propagator has been input. Check input data. 

Message: SETIC, IND. PHASE XXXXXX NOT FOUND 
Source/Tvpe: SETIC/fatal 

Condition/Corrective Action: Phase number XXXXXX was requested as a control 
parameter phase but does not exist. Check the phase numbers input in ENDPH and 
INDSPH. 

Message: SRADRX - BAD ORDERING, INDXI/D/LC 
Source/Tvp e: SRADRX/fatal 

Condition/Corrective Action: The specified indices of the independent and 
dependent variables are not in ascending order when using the 
targeting/optimization algorithm. Check the values of INDXI and INDXD. 

Message: SRADRX - READAT - INDXI/D/LC ARE ZERO 
Source/Tvpe: SRADRX/fatal 

Condition/Corrective Action: The specified indices of the independent and 
dependent variables are zero when using the targeting/optimization algorithm. 
Check the values INDXI and INDXD. 

Message : TGOEMI PROBLEM 
Source/Tvpe: TGOEMI/fatal 

Condition/Corrective Action: No event criteria have been input or too many pending 
events have been requested. Check the value of each CR1TR variable. 

Message: THE TIME OF PERIAPSIS WAS NOT FOUND WITHIN MAX 
ITERATIONS IN SUBROUTINE ONESTP 
Source/Tvpe: ONESTP/fatal 

Condition/Corrective Action: The time of periapsis was not found within the 
maximum number of iterations = 50. 
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Message : TSPXM - BAD IV-TO-END 
Source/Tvpe: TSPXM/warning 

Condition/Correctiv e Action: The data and the dictionary do not match. 
Message: *TTS CONTROL DIVERG IN NEWTON RAPHSN* 

Source/Tvpe: TTS/fatal 

Condition/Correctiv e Action: Select alternate initial conditions and restart. 


Message: *TTS CONVERGENCE NOT REACHED IN TARG1S 
Source/Tvpe: TTS/fatal 

Condition /Corrective Action: Select alternate initial conditions and restart. 


Message : *TTS DIFF NO. DEP. INDEP VARIABLES ♦ 

Source/Tvpe: TTS/fatal 

Condition/Correctiv e Action: Number of subproblem independent controls has to 
equal number of subproblem dependent target variables. Check input data. 

Message : *TTS FUNCTION DIVERG IN NEWTON RAPHSN* 

Source/Tvpe: TTS/fatal 

Condition/Corrective Action: Select alternate initial conditions and restart. 


Message: *TTS ITER EXCEEDED IN NEWTON RAPHSON * 
Source/Tvpe: TTS/fatal 

Condition/Corrective Action: Select alternate initial conditions and restart. 


Message: *TTS PROPAGATION MODEL NOT PROVIDED * 

Source/Tvpe: TTS/fatal 

Condition/Corrective Action: Correct input data. Allowed values are NRAPH'. 
TARG1S', NULL’, or 'NONE’. 


Message: *TTS TARG1S TARGETS NOT AT SAME PHASE* 

Source/Tvpe: TTS/fatal 

Condition/Corrective Action: Redefine inputs. All targets for one subproblem occur 
at the same phase. 


Message: *TTS WRONG NUMBER (NOT 3) INDEP VARIA* 

Source/Tvpe: TTS/fatal 

Condition/Corrective Action: Three subproblem independent control variables 
required. Check input data. 
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6.5 TPOST POST -PROCESSING 

Machine proce ssin g of EPOST generated data is in the form of IPPOST. IPPOST is 
a tool which processes an IPOST (and certain versions of POST) generated PROFIL 
datafile. It allows the user to display PROFIL data in various forms. IPPOST is an 
interactive tool which makes use of the NCAR Graphics package. There are 

approximately 55 co mmand s which provide capability ra ngin g from a simple data 
echo to generating an x-y plot of any 3 variables versus any fourth, independent, 
variable. Access to the command list, and command descriptions, can be obtained 
by entering "man" after IPPOST is brought on-line. In addition, a help menu i6 
available. 
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APPENDIX A: ACRONYMS 


BDT 

BDR 

CP 

CRAF 

DEC 

DOF 

DSM 

DT 

GEO 

GLOW 

GMAP 

HEO 

IMLEO 

IMU 

IPOST 

IPREP 

ISP 

rroo 

ius 

LEO 

LPREP 

NPOST 

NPSOT6D 

NPSOL 

PGA 

POST 

RA 

RCA 

RPERI 

RTN 

S/C 

SEP 

SOI 

SPDY 

SQP 

SVD-LUD 

TCM 

TFP 

TOF 

TPBV 

UVW 

VINF 


B dotT 
B dot R 

Computer Processor 

Comet Rendezvous Asteroid Flyby 

Declination 

Degrees of freedom 

Deep Space Maneuver 

Delta time 

Geostationary orbit 

Gross Liftoff Weight 

General Mission Analysis Program 

High Earth Orbit 

Inserted/injected Mass in Low Earth Orbit 
Inertial Measurement Unit 

Interplanetary Program to Optimize Simulated Trajectories 
Interplanetary PREProcessor 
Specific Impulse 

Interplanetary Targeting and Optimization Option 

Inertial Upper Stage 

Low Earth Orbit 

Lunar PREProcessor 

NPSOL POST 

NPSOL POST6D 

Non-linear Programming Stanford Optimization Laboratory 

Projected Gradient Algorithm 

Program to Optimize Simulated Trajectories 

Right Ascension 

Radius of closest approach 

Radius of periapsis 

Radius Tangent Normal - Coordinate system defined by the 

radius vector and the orbit nominal vector 

Spacecraft 

Solar electric propulsion 
Sphere of influence 
Seconds/day 

Sequential Quadratic Programming 

Singular value decomposition, LU decomposition 

Trajectory correction maneuver 

Time from periapsis 

Time of flight 

Two-point boundary value 

Coordinate system defined by the S/C velocity and the orbit 
nominal vector 

Velocity at infinity (hyperbolic excess velocity) 
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APPENDIX B; IPREP/LPREP USER'S GUIDE 

B-l INTRODUCTION 

IPREP (Interplanetary PREProcessor) provides rapid grid-search optimization on 
launch and arrival windows, minimum mission delta-v, or mass optimization. 
IPREP uses simple propagators (no integration) and minimal perturbing forces so 
that the trajectory calculations will be as quick as possible. 

IPREP has two modes in which it executes, depending on the current value of the 
variable 'namlst.' If namlst = inter (Namelist top sets the initial namlst value, after 
that, it is set at the end of each following namelist), then IPREP will run an 
interplanetary mission. In this case, IPREP utilizes a simple conic propagator and 
Lambert method solutions to provide the velocity requirements at each planet in an 
interplanetary trajectory. Gravity assist is also modeled as a local conic 
perturbation. A Chebychev polynomial represents low thrust segments. The grid 
search is performed on a user input time window for any body encountered during 
the mission. The best mission is calculated by the minimum cost, which is a user 
selected weighting of each encounter delta-v, arrival velocity, departure velocity, 
and initial mass. 

If namlst = lunar' then IPREP will run a lunar mission. In this case IPREP 
employs patched conics and actual targeting (utilizing a Newton-Raphson 
technique), since the earth-moon problem is more complicated than interplanetary 
conic formulations. Lunar mission conditions are generated to the same fidelity as 
interplanetary conditions. The grid search is also performed on encounter time 
(launch or arrival) and optimization is done on total mission delta-v. 

Interplanetary and lunar namelists can be chained together to run as many 
missions the user desires from one input desk. At the end of each namelist, simply 
set the namlst variable to the type of namlst desired to be read in next. When all 
the namelists are input, simply set the last value of namlst to 'none.' 

The following sections describe input and output for the three namelists in IPREP, 
analytic models and algorithms, and program architecture. Most of the detailed 
models can be found in the IPOST documentation (Reference B-l), and the mission 
analysis context of IPREP, LPREP, and IPOST can be found in the IPOST Mission 
Analysis Guide (Reference B - 2). 
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B-2.0 INPUT 



Variable 

Type 

Default 

Description 

abf (50) 

real 

-l.dO 

Aerobrake mass fraction for each 
encounter (orbit insertion) 

acc(50) 

real 

O.dO 

Acceleration due to low thrust engines 
for each leg used in launch and orbit 
insertion spirals. 

advmax(50) real 

50.d0 

Maximum delta-v allowed for s deep 
space leg arrival. 

alphaw 

real 

30. dO 

Powerplant specific mass for low 
thrust engines (kg/kw) 

altatm (22) 

real 

various 

altitude of the atmosphere for the 
planets. Used for calculating entry 
velocity during aerobraking. (km) 

bb 

real 

O.dO 

Used with dd in the low thrust 
efficiency equation = bb/(l+(dd/c)2) 

c(4,22,6) 

real 

various 

Coefficients for analytic ephemeris 
time dependent elements. The first 
array element defines the time based 
coefficient, the second the planetary 
number, and the third, the orbital 
element, a, e, - M. 

cm ax 

real 

l.d38 

Maximum value of the cost function 


allowed. Filtering is performed on all 
trajectories whose cost function is 
greater than this max value. 
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cntrl 


character 


cs real 

C8at(7,20) real 


date(6) real 


'auto' s'auto', grid is refined about the 

minimum trajectory until the 
maximum number of iterations 
(number of times the grid is to be 
refined) has been met. Only the 
minimum cost trajectory is 
printed out. 

= 'manual', grid is not refined. Data 
for each trajectory is saved and 
printed out. This will be useful for 
plotting and deter mining minimum 
trajectories 'manually'. 

= 'sort', and initial pass through 
the grid is made in the manual 
mode, the data is sorted by the 
minimum cost function, a user 
input number (nsave) is saved and 
the grid is refined about each of 
these sets of points. This is to aid 
in avoiding local minimnma 

O.dO Exhaust velocity for the lowthrust 

engines. 

O.dO Orbital elements of planetary 

satellites for calculation of satellite 
closest approaches. 
csat(14) = semi-major axis (km) 
csat(24) = eccentricity 
C8at(34) = inclination (radians) 
csat(44) -longitude of the node (rad) 
C8at(54) = argument of periapsis (rad) 
csat(6,i) = mean anomaly (rad) 
csat(74) = epoch time (Julian date) 
i is the index of the satellite given by 
isat(j,k). The number of satellites for 
each encounter is given by nsat. 

O.dO Epoch time can be input as a 

calendar date as follows: 

(1) = year, 

(2) s month, 

(3) = day, 

(4) = hours, 

(5) = minutes, 

(6) = seconds. 


dd 


real 


O.dO 


it 
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Used with bb in the low thrust 
efficiency equation 


dvmax 

real 

lO.dO 

Maximum delta-v allowed for any 
maneuver. Filtering is performed 
on any trajectory which has a 
maneuver delta-v greater than 
this value. 

ep 

real 

O.dO 

Epoch time from which launch 
and flyby times are measured 
(julian date). 

frac(50) 

real 

O.dO 

Mass fraction included in the 
mass calculation for each 
maneuver. The tankage mass is 
calculated to be this percentage of 
the propellant plus payload mass 

frame 

character 

ecliptic 

Input/output frame. Default of 
ecliptic, or can be 'equator', and the 
frame will be each planet's equator. 

gmu(0:22) 

real 

various 

Gravitational parameter for planetary 
bodies. 

| 

1 

o 

w 

real 

O.dO 

Gravitational parameter for each 
satellite input in csat. 

go 

real 

9.806d-3 

Gravitational acceleration (km/s 2 ) 

har (50) 

real 

-l.dO 

User input arrival vinf magnitude for 
low thrust encounters, if desired. 
Otherwise, it will be calculated. 

hdep (50) 

real 

-l.dO 

User input departure vinf magnitude 
for low thrust encounters. Otherwise, 
vinf will be calculated. 

hinc(50) 

real 

-l.dO 

Inclination (deg) of the orbit for 


each encounter. For launch and 
orbit insertion maneuvers, the 
inclination is input, and for 
swingby maneuvers the 
inclination is calculated. 

= -1., the inclination for the 
launch and orbit insertion 
maneuvers will be set equal to the 
declination of the launch or 
arrival asymptote, respectively. 
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ibod(50) integer 


O.dO 


idvmax(50) 

real 

50.d0 

iephem 

integer 

0 

iepoch 

character 

julian 

ilnch(50) 

integer 

1 

iorbit(50) 

integer 

1 

ipbod(50) 

integer 

0 

iprep 

integer 

2 


Encounter body for each leg. 

= 0, Sun 
= 1, Mercury 
= 2, Venus 
= 3, Earth 
= 4, Mars 
= 5, Jupiter 
= 6, Saturn 
= 7, Uranus 
= 8, Neptune 
= 9, Pluto 

Maximum delta-v for a deep-space leg 
departure. Number of satellite bodies 
input for each encounter. See isat and 
csat. 

= 0, use analytic ephemeris 
= 1, use JPL precision ephemeris. 

Flag setting epoch time input as either 
'ep’ or 'date' 

= 'julian' , input ep as julian date 
= 'calend', input date as calendar date, 
program will convert. 

Flag designating type of launch 
maneuver. 

= 1, high thrust launch 
= 2, low thrust launch 

Flag designating type of orbit insertion 
maneuver. 

= 1, high thrust orbit insertion 
= 2, low thrust orbit insertion 

Central body (primary body) for 
each leg of the mission. Can have 
the same values as ibod above. 

Print flag which determines the 
type of printout. 

= 1, short form of printout 
including only the cost function 
value and the encounter times. 

See also 'outfile'. 

= 2, long form of printout 
including detailed trajectory 
calculations. 
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isat(10,50) 

integer 

0 

Satellite encounter bodies for each leg. 
E.G.: isat (2,3) = 4. User must input 
csat(l:7,4). This satellite is the second 
moon of the third encounter body. The 
4 designates which body in the 
satellite ephemeris to look at. 

itrmax 

integer 

5 

Maximum number of grid 
refinements desired. 

manvr(50) 

character 

’none’ 

Type of maneuver to be performed 
at each encounter. 

= launch', perform launch from 
park orbit 

= 'swngby', perform delta-v flyby 
= ’orbins', perform insertion into 
park orbit. 

namlst 

character*8 

'none' 

Designates next namelist to be read. 

ncone (50) 

integer 

4 

Number of cone angles for each low 
thrust leg 

nconop (50) 

integer 

2 

= 0, all cone angles allowed 
= 2, program chooses the best ncone 
cone angles 

ncop (50) 

integer 

1 

= 0, input cs for each low thrust leg 
s 1, optimizes cs for each low thrust 
leg 

nleg 

integer 

0 

Total number of segments to be 
performed in the trajectory 
sequence. 

nmuop (50) 

integer 

1 

=0, input wmu for each low thrust leg 
si, optimize wmu for each low thrust 
leg 

npow (50) 

integer 

0 

The type of power supply for each low 
thrust leg. 

= 0, constant 
= 1, variable 

nrev(50) 

integer 

0 

The number of revolutions the S/C 
makes during a particular 
trajectory segment. 
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nsat(50) 

real 

0 

nsave 

integer 

10 

nsig(50) 

integer 

7 

outfile 

character 

'none' 

planet(22) 

plosp 

character 

real 

various 

l.dO 

ploss 

real 

O.dO 

plost 

real 

l.dO 

Prop 

character 

conic 


Number of satellite bodies input for 
each encounter. See isat and csat. 

The number of minimum 
missions from the initia l grid 
sweep which will be optimized. 
Required when cntrl = 'sort'. 

The number of significant digits 
required in calculating the optimal 
deepspace maneuver. The solution is 
very sensitive to this number and can 
be greatly affected by its value. Since 
interplanetary problems can be vastly 
different, changing this variables 
value is strongly encouraged. 

Name of the output file in which 
the trajectory information is 
stored after the run (only required 
when iprep = 1). 

Names of the planetary bodies 

Degradation level of power at failure 
in low thrust. New power equals old 
power times plosp. 

Time variation in power of low thrust 
engines, p = po e *(ploss)t 
where pO is limited power at 

Specifies a fraction of the trip where 
discontinuous power loss occurs. 

Propagator selected 


= ’conic', calculate simple Lambert 
method conic solutions between 
encounter points 

* lowthr’, use Chebytop method to 
calculate controlled low thrust 
trajectories between encounter points 
= hplane', calculate a deepspace 
maneuver along with conic 
trajectories between encounter points 
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pyload 

real 

O.dO 

rapoap(50) 

real 

O.dO 

re(22) 

real 

various 

rperi(50) 

real 

O.dO 

rpmax(50) 

real 

O.dO 

rpmin(50) 

real 

O.dO 

scmass 

real 

-l.dO 

spi(50) 

real 

l.d6 

tanks 

real 

O.dO 

thrust(50) 

real 

l.d6 

tat 

real 

90.d0 


Required payload at the end of the 
mission. The initial S/C mass will 
be 'backed out' from this required 
mass. If not input then will be 
calculated from the initial value 
for scmas8. (kg) 

Radius of apoapsis for park orbit 
before launch or after orbit 
insertion (km). 

Equatorial radius of the planetary 
bodies. 

Radius of periapsis for park orbit 
before launch or after orbit 
insertion (km). 

Maximum radius allowed for 
swingby closest approach (km). 

Minimum radius allowed for 
swingby closest approach (km). 

If input is greater than 0, this will be 
initial S/C mass 

Specific impulse for the spacecraft 
engines for each maneuver (seconds). 

Tankage mass fraction for low thrust 
engines. 

Thrust for the spacecraft engines 
for each maneuver (newtons). 

Angle the plane of the solar panels 
make with the S/C sun vector. 
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tm(3,50) real O.dO 


tmax real l.d6 


tol(5,50) real l.d38 


*tfett(2,50) 

real 

O.dO 

wmu(50) 

real 

O.dO 

wpmax 

real 

l.d38 


tm( 14 ) = minimum time of flight 
for the ith leg 

tm(2,i) = maximum time of flight 
for the ith leg 

tm(34) = time increment for the 

ith leg, which defines the grid 
size 

If tm(14) = -1, then the user has 

chosen to fix this encounter time 
by inputing the julian date of 
the encounter into tm(24) and 
zero into tm(3,i). 

Maximum total trip time in days. 
Filtering is done on missions 
whose total trip time exceeds this 
maximum value. 

Weighting used in calculation of 
the cost function for each 
segment. The cost function is 
defined more extensively in the 
analytic section for IPREP. 
tol(l ,i) = weighting on departure 
v-inf magnitude. 

tol(24) = weighting on arrival v-inf 
magnitude. 

tol(34) = weighting on departure 
mass. 

tol(44) = weighting on arrival mass. 
tol(54) = weighting on encounter 
delta-v (launch, swingby, or 
orbins) 

tol(64) = weighting on deepspace 
maneuver delta-v. 

tol(74) = weighting on time of flight for 
current leg. 

Jettison mass for deep-space bums. 
wjett(l,i) = jettison mass. 
wjett(24) s +I.d0, jettison after bum. 

= -ld.O, jettison before bum. 

Powerplant mass fraction for each of 
the lowthrust legs. 

Maximum amount of propellant 
allowed. 
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B-2.2 T.TTNAR NAMF.T.TST INPUTS 


Variable 

Type 

Default 

Description 

ep 

real 

O.dO 

Epoch time upon which all 
encounter times are based (julian 
date). 

hinc(l) 

real 

O.dO 

Indination of park orbit about 
initial body, whether primary body 
(istyle = 1,2, 3, or 4) or secondary 
body (istyle = 5). 

hinc(2) 

real 

O.dO 

Inclination of orbit about second 
encounter body, whether 
secondary body (istyle = 1,2,3, or 4) 
or primary body (istyle = 5). 

hinc(3) 

real 

O.dO 

Indination of orbit about third 
encounter body which will be the 
primary body when istyle = 3 or 4. 

iephem 

integer 

1 

= 0, use analytic ephemeris 
= l, use JPL prerision ephemeris 

ipbod 

integer 

3 

Primary body id from ibod array. 

iprint 

integer 

0 

= 0, summary print 
= 1, full print 
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istyle 


integer 


1 


itbod 

lib 


integer 


10 


integer 0 


Controls the program operation. 

= 1, launch from primary park 
orbit and target to conditions at the 
secondary (one burn) 

= 2, launch from primary park 
orbit, target to secondary, enter 
park orbit about secondary (two 
bums). 

= 3, launch from primary park 
orbit, target to secondary, enter 
park orbit about secondary, 
launch from secondary park orbit, 
target to primary, enter park orbit 
about primary (four burns). 

= 4, launch from primary park 
orbit, target to secondary is such a 
way as to return to primary with a 
zero delta-v swingby. Provides a 
free return trajectory. 

= 5, launch from secondary, target 
to primary, enter orbit about 
primary. 

= 6, start with p re cessing park orbit 
about primary, and calculate all of 
the coplanar launch opportunities 
during the user input time frame. 
Targeting of lunar closest approach 
radius and inclination. 

= 7, start with precessing park orbit 
about the primary, and calculate all 
of the coplanar launch opportunities 
to the user input libration point 
during the requested time frame. 

= 8, start with park orbit about the 
moon, processing S/S orbit about 
Earth and calculate coplanar return 
missions from the moon. 

= 9, calculate missions from moon to 
libration points. 

= 10, roundtrip with coplanar launch 
from Earth, orbit insertion at moon, 
launch from moon with coplanar 
arrival at Earth. 

Secondary body id . 

The libration requested to be targeted 
to. 
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mxitar 

integer 

100 

npi 

integer 

5 


orb(l:6) 

real 

O.dO 

ra(l:3) 

real 

O.dO 

rp(l:3) 

real 

O.dO 

’istyle' = 1 or 2 

then 


tm(l,l) 

real 

O.dO 

tm(2,l) 

real 

O.dO 

tm(3,l) 

real 

O.dO 

tm(l,2) 

real 

O.dO 

tm(2,2) 

real 

O.dO 

tm(3,2) 

real 

O.dO 

told) 

real 

l.dO 

tol(2) 

real 

l.dO 


Maximum number of iterations to 
be allowed in the Newton-Raphson 
targeting scheme. 

Number of previous iterations in 
Newton-Raphson to check against 
for divergence, either in 
increasing target error, or 
increasing control changes. 

Orbital elements at epoch when ilnch 
= 6,7. 

Same as hinc above, except applies 
to the radius of apoapsis for park 
orbits. 

Same as hinc above, except applies 
to the radius of periapsis. 


Days from epoch for initial launch 
from the primary. 

Days from epoch for final launch 
from the primary. 

Time increment to step through 
the above region. 

Minimum time of flight to the 
planet (total trip time). 

Maximum time of flight to the 
secondary SOI. 

Time increment to step through 
the above region. 

Weighting on launch delta-v from 

primary for calculating the cost 
function. 

Weighting on orbit insertion Av at 
secondary for calculating the 
cost function. 
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If 'istyle' = 3 then add 


tm(l,3) 

real 

O.dO 

tm(2,3) 

real 

O.dO 

tm(3,3) 

real 

O.dO 

tm(l,4) 

real 

O.dO 

tm(2,4) 

real 

O.dO 

tm(3,4) 

real 

O.dO 

tol(3) 

real 

l.dO 

tol(4) 

real 

l.dO 


Or if ’istyle' = 4 then add 


tm(l,3) 

real 

O.dO 

tm(2,3) 

real 

O.dO 

tm(3,3) 

real 

O.dO 


Minimum stay time at secondary 

Maximum stay time at secondary 

Time increment to step through 
the above region. 

Minimum time of flight on 

departure from secondary SOI 
to primary arrival. 

Maximum time of flight on 

departure from secondary SOI 
to primary arrival. 

Time increment to step through 
the above region. 

Weighting on launch delta-v from 
secondary for calculating the 
cost function. 

Weighting on orbit insertion Av at 

primary for calculating the cost 
function. 


Minimum time of flight on 
departure from the secondary SOI 
to primary arrival. 

Maximum time of flight on 
depanire from the secondary SOI 
to primary arrival. 

Time increment to step through 
the above region. 


If 'istyle' = 5, then the time and weighting variables are defined as 


tm(l,l) 

real 

O.dO 

Days from epoch for the initial 
departure from the secondary SOI. 

tm(2,l) 

real 

O.dO 

Days from epoch for the final 
departure from the secondary SOI. 

tm(3,l) 

real 

O.dO 

Time increment to step through 
the above region. 

tm(l,2) 

real 

O.dO 

Minimum time of flight on 
departure from the secondary SOI 
to primary arrival. 

tm(2,2) 

real 

O.dO 

Maximum time of flight on 
departure from the secondary SOI 
to primary arrival. 

tm(3,2) 

real 

O.dO 

Time increment to step through 
the above region. 

told) 

real 

l.dO 

Weighting on the launch delta-v 
from secondary for calculating the 
cost function. 

tol(2) 

real 

l.dO 

Weighting on the orbit insertion 
delta-v at the primary for 
calculating the cost function. 


If ’istyle' = 6 or 7, then the time variables are defined as 


tm(l,l) 

real 

O.dO 

Days from epoch for initial launch 
from the primary. 

tm(2,l) 

real 

O.dO 

Days from epoch for final launch 
from the primary. 

tof 

real 

O.dO 

Time of flight from LEO to arrival 
at LLO or one of the libration points 
when istyle = 6,7. 


If 'istyle' = 8 or 9, then the time variables are defined as 


tm(l,l) real 

O.dO 

Days from epoch for initial launch 
from secondary. 

tm (2,1) 


Days from epoch for final launch from 
secondary. 

tm (3,1) 


Time increment to step through the 
above region. 

tof 


Time of flight from LLO to LEO as one 
of the libration points. 
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If istyle = 10, then add to the information under istyle = 6 


tm(l,2) 

real 

O.dO 

Days from arrival at secondary for 
initial launch to primary. 

tm(2,2) 

real 

O.dO 

Days from arrival at secondary for 
final launch to primary. 

tof(2) 

reed 

O.dO 

Time of flight from LLO to LEO. 
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B-2.S SPIRAL NAMELIST INPUTS 


al 

real 

O.dO 

Semi-mqjor axis of initial orbit (km). 

&2 

real 

O.dO 

Semi-mqjor axis of final orbit (km). 

argp 

real 

O.dO 

Argument of periapsis for initial orbit. 

el 

real 

O.dO 

Eccentricity of initial orbit. 

e2 

real 

O.dO 

Eccentricity of final orbit. 

hind 

real 

O.dO 

Inclination of initial orbit. 

hinc2 

real 

O.dO 

Inclination of final orbit. 

ipri 

integer 

0 

Central body id. 

longn 

real 

O.dO 

Longitude of node for initial orbit. 

namlst 

character*8 

’none’ 

Designates next namelist to be read. 

nrev 

integer 

1 

The number of revolutions required for 
the spiral. 

nt 

integer 

0 

= 0, no tracking from previous 
solutions from one rev to the next 
= 1, track from previous solutions. 

time 

real 

O.dO 

Julian date for commencement of the 
spiral maneuver. 
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B-3.0 SAMPLE CASES 


The following sections give some sample cases for IPREP. These cases, along with 
the above description of the namelist inputs should give the user a good idea of how 
to set up their own input decks and run this program. 

There are two interplanetary sample cases: The first is a Mars rountrip mis sion 
and the second is a low thrust mission to Jupiter. The roundtrip trajectory to Mars 
has a Venus flyby on the outbound leg. The weighting on the cost function values is 
set so the total mission delta-v will be minimized, with a strong emp ha sis on the 
delta-v at the Venus swingby. The purpose is to find the minimum mission which 
also has a zero delta-v at the flyby. As can be seen from the output on the next page 
(section B-3.1), the Venus swingby delta-v was negligible, and the total mission 
delta-v was 12.0782 km/s. 

The first LUNAR sample case is a lunar free return mission. The minimum cost in 
this case will be the minimum total mission delta-v, which happens, in this case, to 
be only the launch delta-v. As can be seen in the output in section B-3.2, the S/C 
flew by the moon at an altitude of 200.0 km ( the radius of the moon is 1738 km ), 
and returned to the Earth with an altitude of 200 km. These values can be changed 
in the input by the user. 

The second LUNAR sample case is a departure from a space station orbit to LLO. 

All of the coplanar launch opportunities for the input period of time are printed out. 

The first SPIRAL sample case is an Earth-centered trajectory going from LEO to 
GEO, maintaining dose to zero eccentridty. The rev counter was set to 20. 

The second SPIRAL sample case is an Earth-centered trajectory going from LEO to 
ahighly eccentric orbit (e2 = .8) using the tracking capability of Chebytop (nt = 1). 
The rev counter for this case was also set to 20. 
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B-3.1 TPREP SAMPLE CASES 
B-3.1.1 MARS ROUNDTRIP 
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Interplanetary Trajectory Case 
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**Translunar Trajectory Case** 
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***** Trans l unar Trajectory Case***** 
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★♦★♦♦Spiral Orbit Transfer Case***** 
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spiral Orbit Transfer 


•Spiral Orbit Transfer Case***** 
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Total Time for Spiral Maneuver * 0.997209E+01 days 


B-4.0 ANALYTIC MOnFT.fi 


B-4.1 




SR TRAJECTORY SIMULATION 


Given the times of encounter of each of the planets in the trajectory sequence, a sun 
centered come approximation can be calculated from the ephemeris positions of any 
two subsequent planets and the time of flight between the encounters. With the 
positions and the time of flight known, a simple Lambert-type solution will be 
calculated to provide the sun-centered velocity at each of the encounter points. The 
Lambert algorithm can be found in reference B-3. Once the sun-centered velodties 
are calculated, then the escape velocity ( Voo ) can be obtained by subtracting away 
the planetary velocity from the heliocentric a/c velocity. 


Voo = v g/c " v planet 

The Voo's of arrival and departure at each planet can then be utilized to calculate 

launches, orbital insertions, and gravity-assist swingbys. A generalized gravity- 
assist mission is shown in figure 1. 



figure 1. Generalized Interplanetary Mission 


In this mission, a laundi maneuver is performed from planet 1. Two gravity-assist 
swingbys are performed at planets 2 and 3, which may, or may not require A V"s 
depending on the planetary geometiy. Finally, an orbital insertion maneuver is 
performed at planet 4 The heliocentric portions of the trajectory are difined by the 
ephemens positions of the planets, and the Voo’s calculated by the above equation. 
The launch, jswingby , and orbital insertion maneuvers will be discussed in the 


IPREP also has the capability of calculating low thrust trajectories between two 
encounter points. Applying the core programs used in QT2 (reference B-6) IPREP 
approbates the powered heliocentric trajectory using Chebychev polynomials. 
Output gives S/C acceleration vector magnitude and direction due to the low thrust 
engines. The thrust acceleration direction output from the low thrust routines is 
given m the cone/clock angles system shown in Figure 2 -9. The cone angle for the 
acceleration vector may be constrained from 1 to 4 constant values or allowed to 
vary continuously along with the clock angle and magnitude. 
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Since IPOST does not accept cone and dock angles for S/C orientation, IPREP 
converts the acceleration vector to magnitude, yaw, and pitch in the heliocentric, 
ecliptic system. During the periods where the cone angle is fixed, cubic polynomial 
coeffidents are calculated for magnitude, yaw and pitch. These polynomials are 
time-dependent and fit into IPOST inputs. A sample low thrust case is shown. 
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B-4.2 LAUNCH MANEUVER 


Since we know the radius of periapsis ( r p ), the radius of apoapsis ( r a ), and the 
inclination (i) of our park orbit ( these are all user inputs ), along with the v°° which 
we calculated from the Lambert solution, we can obtain the A required at the launch 
planet to reach the target planet. 


v-infinity unit vector, 
semi-major axis, 


C 3 ■ Vqq 2 
vu = unit ( Voo ) 


declination of escape asymptote, 

6 = a tan 


-(£) 


(vvu x 2 + TUy 2 ) 


<)> = atan — 


park orbit periapsis velocity, 

v " ark = (i-p + r. )) 

escape periapsis velocity, 


v n =• 


bv 0 


If the inclination of the park orbit is greater than the declination of the departure 
asymptote then no plane change is required and the Av is simply calculated 

Av = abs ( v p - vparfc ). 

If the inclination of the park orbit is less than the declination of the departure 
asymptote, then a plane change is required and the Av is calculated 


Av pc = 2v p sin 


and 


. fabsm-i"! 

J 

Av = abs ( v p - Vpark ) + Av^. 
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B-4.3 SWINGBY MANEUVER 


Since we know the incoming and outgoing v-infinity vectors to the swingby planet, 
we can calculate the required delta-v at the flyby, if one is needed. The two 
asymptotes define the required turn, which in turn defines the radius of closest 
approach. If the radius of closest approach is not within the mission limits the user 
haa d efin ed ( beneath the surface of the planet for instance ) then a delta-v will be 
required to achieve the proper turn angle with the user defined radius. Also, if the 
magnitudes of the asymptotes are different, a delta-v will be required to give the 
proper escape velocity. 


required turn angle, 


flyby radius, 


v.o(in) = I VooGn) I 

Voo(out) = I V °»(Ollt) I 

= cos . 1 fMin) v-j o um 
V v °o(in) v oo(out) ) 


(( l _ a w 


rca = |xl 




1 


V v °° (out) / 


if rca > r max then r ca = r max , or 
if r ca < r min then r^ = r m j n , 


where r max and r m j n are the user input botmds for the closest approach. If one or 
the other of these cases is true, then the needed turn angle cannot be achieved 
within the user specified bounds for closest approach. A delta-v is required to get 
the proper turn angle, 


<j) m = 2 sin 


1 + 


r ca v oo^(out) 


where f m is the turn angle the planets gravity is capable of within the bounds of 
r max r min* The delta-v required to make up the turn angle difference is then 

Av = V v oo 2 ( in) + v« 2 (out) -2v-(in)V«(out)«>8(<l> * • 


If a delta-v is not required to change the turn angle, then the delta-v calculation is 
made only for the change in magnitude, 

Av = I v^in) - v„o( 0U ^) I . 
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Also included in IPREP is a periapsis bum optimization method for minimizing the 
delta-v of the swingby. The specifics of this method can be found in reference B-4. 
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The orbit insertion maneuver is calculated using the v,* vector, the inclination of the 
desired park orbit, the radii of periapsis and apoapsis of the desired park orbit. 

v„ = IO 

rhyp = r p 

where r^yp is the hyperbolic periapsis, and r p is the desired periapsis for the park 
orbit. 

<j) = 2 sin 





p Nf rhyp ®i 
Av l = Vph -Vp 

the declination of the incoming asymptote ( A ) is greater than the inclination of 
the park orbit ( i ), then there is also a plane change required, 


Av 2 = 2 v^ sin 

Now, the total delta-v required is 

Av = Avj + Av 2 . 
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B-4.5 OPTIMUM MIDCOURSE - DESCRIPTION OF TFTF. 
PROBLEM 


In the coplanar case, the evaluation of the optimal transfer is easily accomplished 
using the solution of Lambert's problem. However, if the plane of the arrival planet 
i6 sufficiently inclined, the point at arrival is close to the peak of this inclined 
surface, and the initial and final trajectory points are at a maximum angular 
displacement with respect to the central attracting body, one can be assured that 
the single impulse solution will be non-optimal. In these and other slightly 
perturbed casesd, a two impulse solution with one of the impulses being a midcourse 
maneuver is always a superior solution. Numerous studies and implementations of 
analytic solutions have been developed to identify the placement and timing of this 
midcourse maneuver. However, it is questionable whether these analytic solutions 
are sufficiently robust to find optimal solutions when the objective function diverges 
significantly away from the one originally used in the derivation of the analytic 
solution. To eliminate this uncertainity, the authors developed a method of "stacked 
optimizers" which is not as mathematically elegant as many of the closed-form 
solutions derived in recent years but does remove the dependence upon the 
derivation of the solution. 


To employ the staked optimizer method, the problem pf the optimal midcourse 
maneuver needed to be divided up into a set of independent subproblems. Each 
subprobiem would need to converge upon a local optimum before its result would be 
fet to the upper level optiimzation. The stacked optimizer approach essentially uses 
three levels of optimizers. The initial level is a one-dimensional minimizer that 
requires a bracketed range to enclose the domain of acceptable control variables. At 
this level the control variable is the time of the midcourse maneuver and is bounded 
by the launch and arrive dates. This level of optimization simply tests values of a 
funcfconbyperturbing values of midcourse maneuver timing until a minimum value 
for the objects function, usually delta velocity is found. The function that is tested 
by the first level optimizer is itself a two dimensional optimizer. This second level 
optimizer uses the heliocentric dedination and right ascension as control variables 
Again, these variables are perturbed until a set of values are found that result in a' 
converged objective function value. These perturbed values are used as constants in 
a final thnd level optimizer. In this final level of optimization the position 
magnitude is optimized in the direction defined by the higher level angular direction 
components until a minimum objective value is identified. In all three levies, the 
objective function value is the same as that evaluated by the third level magnitude 
optamzer, it is simply carried to the top as each optimzer converges. 


U is important to note that the method used to calculate the objective function does 
not affect the optimization method used to determine where the original midcourse 
ensts, but that differences to the objective function will change where the optimal 
midcourse is placed. Essentially, this method converges on the mimimum 
regardless of how thwe objective is calculated. In this manner, each level of 
optimization is itself given a converged solution for optimal delta velocity. As each 
subproblem is satisfied with a converged solution, the next higher level recomputes 
with perturbed control variables until a global minimum is found. 
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i 


This numerical method lends itself easily to modifications of the obective function 
evaluation since the stacked optimizers will converge on the solution regardless of 
the methos used to generate the objective function. Therefore, if the mission analyst 
wished to optimize on delta velocity injection, modcourse, arrival, or any 
combination of these in a scaled or unsealed manner, he/she could expect the 
method to converge upon the true minimum for his/her problem. Additionally, the 
implementation used by the authors allows the user to constrain many of the delta 
velocities. Again, this method has no difficulty re-evaluating midcourse maneuvers 
to satisfy the constraints. 
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B-4.6 MASS MODELING 


For each delta-v which occurs in the interplanetary trajectory, there is also a delta- 
mass which is calculated by the rocket equation. 

mf = nrijft*Av/(gIf>) 

where mj is the initial mass before the bum, mf is the final mass after the bum, Av 
is the delta-v of the bum, Isp is the specific impulse of the rocket engines, and g is 
the acceleration due to gravity ( 9.806 x 10* 3 km/s 2 ). Since the change in mass is 
entirely propellant, the mass of the tankage to hold this propellant also needs to be 
calculated. Therefore 


m tank = f x m prop 
“prop = mi - mf 

and 

f = tankage mass fraction. 

Therefore the final mass for any trajectory leg bum will be 

mf = mj - m prop - m^^ - mjettison 

where mp^gon is the mass of any structure, probe, or payload jettison for the 
current leg of the mission. 

This process is repeated for each bum until a final mass is calculated for the entire 
mission, which will be the payload for the vehicle, mp a y] oa( j. This value can be 
maximi zed at arrival. Also within IPREP is the capability to input rup^]^ and 
back out what the initial mass would be by reversing the above calculations. Then 
the initial mass can be minimized when a specific payload is reouired for the 
mission. 
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B-4.7 COST MODELING 

The cost equation in IPREP for each interplanetary trajectory is as follows. 


n 


cost = 



v«(in) maS8(out) maS8 (in) 
tol 2 i + tol 3i tol 4 i 


+ 


tol 5i J 


where the toll's are the weightings on each of these cost variables for each 

trajectory leg. The defaults for each of these weights is + °° so the cost function will 
only be dependent on the variables for which the user inputs tolerances. 
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B-4.8 lprep trajectory simulation 


Since the Earth-Moon problem is much more complicated than a similar 
interplanetary problem, simple conic trajectories have insufficient precision to f ulfill 
the accuracy requirements. By allowing only the Earth to affect the trajectory until 
the s/c encounters the lunar sphere of influence (SOI), and then switching to a 
Moon-centered system, assuming it is now only moving under lunar gravitation, we 
can still use 2-body orbital mechanics. This is still an estimate, called the patched- 
conic approximation, but it is sufficient for our needs. 

The setup is much the same as in IPREP. A general Earth-Moon roundtrip 
trajectory is shown in figure B-4.2 below. 



Figure 4.2 Earth-Moon roundtrip trajectory 

In figure B-4.2 a departure is made from low-Earth orbit , the s/c is flown to lunar 
periapsis, an orbit insertion maneuver is performed to place the s/c into low-lunar 
orbit, some time later a launch maneuver is performed to return to perigee «nH 
another maneuver is performed to place the s/c back into low-Earth orbit. 

Another difference between LPREP and IPREP is that LPREP uses a Newton- 
Raphson method to target to actual conditions at the Moon and Earth. The 
development of the Newton-Raphson method can be found in reference B-l. The 
user inputs the park orbit periapsis, apoapsis and indination at both the Earth and 
Moon, and the time of flight from perigee to the lunar SOI ( or vice-versa if a return 
mission is being calculated). A point is chosen on the lunar SOI and a 2-body conic 
trqjectopr is calculated from perigee to this point utilizing methods in reference B-5. 
The trajectory state is then converted to the Moon system and the trajectory is 
again propagated by 2-body conics to perilune. The target conditions are then 
calculated, along with the errors and sensitivity matrix. The point of penetration of 
the SOI is then changed accordingly, until the proper target conditions are achieved. 

The return trajectory is solved the same way as the departure, i.e the Earth-Moon 
trajectory is calculated and then the velocity vectors are reversed. This 
gives you a Moon-Earth trajectory. 
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The free return trajectory is done much the same way as the Earth-Moon problem 
except that the time to the Moon becomes a free variable (along with the 
penetration point on the SOI) and the return perigee becomes a new target. 

The algorithm must find the correct point of penetration and time of flight to not 
only get the proper flyby conditions at the Moon but also the correct return radius. 
No delta-v is performed at the moon, thus constituting a 'free' return. 

In figures B-4.3 and B-4.4, there are two launches from a space station orbit shown. 
In figure B-4.3, the space station orbit and the trans-lunar trajectory are coplanar, 
so there is no need for a plane change maneuver at Earth launch. In other words, 5, 
the declination of the moon relative to the space station orbit at arrival is zero, or 
near zero. Figure B-4.4 shows a translunar trajectory where S is not zero. Past 
studies have shown that the cheapest lunar missions occur when this value of 8 is at 
zero. A plane change at launch is very expensive, so if the mission is limited to a 
coplanar launch, then mission opportunities will occur only when 8 is at, or near, 
zero. Figure B-4.5 is a graph of how 6, the declination of the moon relative to the 
space station orbit, changes over time. The space station orbit is processing with 
time, due to J2 effects, about 7° per day. It can be seen from figure B-4.5 that 
opportunities with coplanar launches occur approximately every 6 to 11 days. 
Trans-lunar trajectory will be flown coplanar, targeting to lunar periapsis, and the 
insertion into lunar orbit (including any necessary plane changes) will be performed. 
Lunar orbital plane changes are relatively cheap. 

Return from the moon to Earth is calculated similarly, except that the necessary 
plane change is performed at launch from the moon with a coplanar arrival into 
space station orbit at the Earth. 

The same missions can be performed to Earth-moon libration points. The libration 
point setup is shown in figure B-4.6. 

The cost of the mission is currently based on the total delta-v. The best mission will 
be the one with the minimum delta-v. There will be improved cost analysis in 
future updates. 
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Figure B - 4.6 Libration Points 
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